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

GACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAG  >  NC_007779/66437‑66629                                                                                                |                                                                                                  gACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccga                                                                                              >  3:683033‑M1/1‑96 (MQ=255)  aCGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgag                                                                                             >  1:696833‑M1/1‑95 (MQ=255)     ccgcGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacc                                                                                          <  3:1147362‑M1/101‑10 (MQ=255)       gcgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccg                                                                                        >  3:58234‑M1/1‑90 (MQ=255)              cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                                 <  1:827311‑M1/101‑19 (MQ=255)              cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                                 <  1:81340‑M1/101‑19 (MQ=255)               gATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtc                                                                                <  1:1022192‑M1/101‑20 (MQ=255)               gATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtc                                                                                <  1:893269‑M1/101‑20 (MQ=255)                   aaCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacc                                                                            <  3:79110‑M1/101‑24 (MQ=255)                    aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacct                                                                           <  3:1197157‑M1/101‑25 (MQ=255)                     cGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctg                                                                          >  1:97869‑M1/1‑76 (MQ=255)                      ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgc                                                                         <  1:1263977‑M1/101‑27 (MQ=255)                      ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgc                                                                         <  3:1264204‑M1/101‑27 (MQ=255)                                ttGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggag                                                               >  3:19566‑M1/1‑65 (MQ=255)                                    cccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaac                                                           <  3:1236473‑M1/101‑41 (MQ=255)                                       aTGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                                        <  3:501450‑M1/101‑44 (MQ=255)                                       aTGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                                        <  1:866331‑M1/101‑44 (MQ=255)                                        tgtgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagta                                                       <  1:940358‑M1/101‑45 (MQ=255)                                           gAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagag                                                    >  3:751502‑M1/1‑54 (MQ=255)                                             gCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagag                                                  <  3:1492010‑M1/101‑50 (MQ=255)                                                tGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttg                                               >  1:853185‑M1/1‑49 (MQ=255)                                                     aGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgata                                          <  3:716565‑M1/101‑58 (MQ=255)                                                       gTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaa                                        <  3:1325298‑M1/101‑60 (MQ=255)                                                         tgtgTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaa                                      >  1:359513‑M1/1‑40 (MQ=255)                                                            gtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcg                                   >  3:1049216‑M1/1‑37 (MQ=255)                                                            gtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcg                                   <  1:1038000‑M1/101‑65 (MQ=255)                                                            gtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcg                                   <  3:760702‑M1/101‑65 (MQ=255)                                                               ttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggaggaacagtagagagttgcgataaaaagcgtca                                >  1:626750‑M1/1‑34 (MQ=255)                                                                   cAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                            <  1:1027495‑M1/101‑72 (MQ=255)                                                                   cAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                            <  1:742311‑M1/101‑72 (MQ=255)                                                                      cgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga                         <  3:1237789‑M1/101‑75 (MQ=255)                                                                      cgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga                         >  1:536308‑M1/1‑27 (MQ=255)                                                                      cgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggaagga                         >  1:854086‑M1/1‑27 (MQ=255)                                                                              tGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaa                 <  3:135384‑M1/101‑83 (MQ=255)                                                                                gCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatc               <  1:37983‑M1/101‑85 (MQ=255)                                                                                     tAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatg          <  3:98954‑M1/101‑90 (MQ=255)                                                                                        tACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggat       <  3:1599697‑M1/101‑93 (MQ=255)                                                                                          cccggcgcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataa     >  3:444568‑M1/4‑7 (MQ=255)                                                                                             ggCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaa  >  3:1210241‑M1/1‑4 (MQ=255)                                                                                                |                                                                                                  GACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAG  >  NC_007779/66437‑66629

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