Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_007779	66,533	3540 bp→CCCGAGACCCGCACGCGTCGACCTGCA		[araD]–araB	[araD], araA, araB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_007779	66533	70072	3540	39 [0]	[0] 30	[araD]–araB	[araD],araA,araB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_007779	= 66532	0 (0.000)	26 (0.770) +CCCGAGACCCGCACGCGTCGACCTGCA	21/146	0.6	100%	coding (19/696 nt)	araD	L‑ribulose‑5‑phosphate 4‑epimerase
	?	NC_007779	70073 =	0 (0.000)					intergenic (‑25/‑314)	araB/araC	L‑ribulokinase/DNA‑binding transcriptional dual regulator

GTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTG  >  NC_007779/69981‑70172                                                                                             |                                                                                                    ttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGgag                                                                                             >  3:19566‑M2/93‑101 (MQ=255)     cccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAAc                                                                                         <  3:1236473‑M2/13‑1 (MQ=255)        atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGt                                                                                      <  3:501450‑M2/16‑1 (MQ=255)        atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGt                                                                                      <  1:866331‑M2/16‑1 (MQ=255)         tgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                                     <  1:940358‑M2/17‑1 (MQ=255)            gagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                                  >  3:751502‑M2/82‑101 (MQ=255)              gcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagagag                                                                                <  3:1492010‑M2/22‑1 (MQ=255)                 tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTg                                                                             >  1:853185‑M2/77‑101 (MQ=255)                      aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATa                                                                        <  3:716565‑M2/30‑1 (MQ=255)                        gttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaa                                                                      <  3:1325298‑M2/32‑1 (MQ=255)                          tgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaaa                                                                    >  1:359513‑M2/68‑101 (MQ=255)                             gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                                 <  3:760702‑M2/37‑1 (MQ=255)                             gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                                 >  3:1049216‑M2/65‑101 (MQ=255)                             gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                                 <  1:1038000‑M2/37‑1 (MQ=255)                                ttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGGAACAGTAGAGAGTTGCGATAAAAAGCGTCa                                                              >  1:626750‑M2/62‑101 (MQ=255)                                    cagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                          <  1:1027495‑M2/44‑1 (MQ=255)                                    cagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                          <  1:742311‑M2/44‑1 (MQ=255)                                       cgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCaggaagga                                                       >  1:854086‑M2/55‑101 (MQ=255)                                       cgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGa                                                       <  3:1237789‑M2/47‑1 (MQ=255)                                       cgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGa                                                       >  1:536308‑M2/55‑101 (MQ=255)                                               tggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                               <  3:135384‑M2/55‑1 (MQ=255)                                                 gcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                             <  1:37983‑M2/57‑1 (MQ=255)                                                      taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATg                                        <  3:98954‑M2/62‑1 (MQ=255)                                                         tacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGAt                                     <  3:1599697‑M2/65‑1 (MQ=255)                                                           cccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaa                                   >  3:444568‑M2/35‑101 (MQ=255)                                                              ggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                                >  3:1210241‑M2/32‑101 (MQ=255)                                                                  cccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCt                            <  1:492610‑M2/74‑1 (MQ=255)                                                                      agacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATgg                        <  1:884112‑M2/78‑1 (MQ=255)                                                                          ccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATa                    >  3:335032‑M2/20‑101 (MQ=255)                                                                                            aGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTg  >  3:1492154‑M2/2‑101 (MQ=255)                                                                                             |                                                                                                    GTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTG  >  NC_007779/69981‑70172

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑