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

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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 66532	0 (0.000)	32 (1.040) +CCCGAGACCCGCACGCGTCGACCTGCA	31/428	0.2	100%	coding (19/696 nt)	araD	L‑ribulose‑5‑phosphate 4‑epimerase AraD
	?	NC_000913	70073 =	0 (0.000)					intergenic (‑25/‑314)	araB/araC	ribulokinase/DNA‑binding transcriptional dual regulator AraC

CTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACG  >  NC_000913/69831‑70317                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     ggcgtgtcggaggagggctttttcggtccttcaaccacttcaccggtttcgatgctaaccacgaccatagcgtcagcgggcatgacgctgtaatcgacgccggaaggtttgagcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccttcgcccgatacccgcacgcgtcggcctgcaGGCATGgag                                                                                                                                                                                                                                              <  2:315063‑M2/9‑1 (MQ=255)        cggaggagggctttttcgtaccttcaaccacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttttaataccgggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGt                                                                                                                                                                                                                                       <  1:127448‑M2/16‑1 (MQ=255)              agggctttttcgtaccttcaaccacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggagctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                                                                                                                                                                                   >  1:42990‑M2/230‑249 (MQ=255)                      ttcgtaccttcaaccacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGAt                                                                                                                                                                                                                          >  1:59475‑M2/222‑250 (MQ=255)                         gtaccttcaaccacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaa                                                                                                                                                                                                                      <  1:81836‑M2/33‑1 (MQ=255)                           accctcaaccacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgctgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgtgttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaTGTATGGATAAACAGTAGAGAGTTGCGATaaaaa                                                                                                                                                                                                                     <  2:308514‑M2/34‑1 (MQ=255)                            cctgcaaccacctcaccggtttcgatgctaaccacgaccatatcgtcagaggtcaggacgctgtaagcgacgccggaaggtgtgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagatccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGc                                                                                                                                                                                                                   <  2:270865‑M2/36‑1 (MQ=255)                             cttcaaccacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                                                                                                                                                                                  >  2:58563‑M2/215‑251 (MQ=255)                                  accacttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAgg                                                                                                                                                                                                             <  2:60971‑M2/42‑1 (MQ=255)                                       ttcaccggtttcgatgctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTAAGGTAgg                                                                                                                                                                                                         >  2:293806‑M2/205‑250 (MQ=255)                                          accggttgcgatggtaaccgcgaccatagcgtcagcgggcatgacgctgtaatcgacgccggaaggtttgatcgcaaagccgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttgggttctaatacccggcgcccgagacccgcacgcgtcgaccggcaGGTATGGAGAAACAGTATAGAGTTGCGATAAAAAGCGTCAGTTaggagcc                                                                                                                                                                                                     <  2:247033‑M2/50‑4 (MQ=255)                                                      gctaaccacgaccatatcgtcagcggtcatgacgctgtaatcgacgccggaaggtgtgagcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggctgctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGATAGTTGCGATAAAAAGCGTCAGGTAGTATCCGCTAATCTTAt                                                                                                                                                                                          <  2:43898‑M2/61‑1 (MQ=255)                                                                     atcgtcagcggtcatgtcgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATg                                                                                                                                                                          <  2:164138‑M2/77‑1 (MQ=255)                                                                     atcgtcagcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGCGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATg                                                                                                                                                                          >  1:286885‑M2/175‑251 (MQ=255)                                                                            gcggtcatgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGc                                                                                                                                                                   <  1:362482‑M2/84‑1 (MQ=255)                                                                              ggtcaggacgctgtaatcgacgccggaaggtttgagcacaaagacgccgcgctcgcgatcaacggcgctgtcgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaa                                                                                                                                                                 <  1:239972‑M2/86‑1 (MQ=255)                                                                                 catgacgctgtaatcgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgcggacgttgccccatgtgagcgggaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgt                                                                                                                                                              <  2:272341‑M2/89‑1 (MQ=255)                                                                                            aatcgacgccggcaggtttgatcacaaagacgccgcgctcgcgctcaaaggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGCAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTg                                                                                                                                                   >  2:213075‑M2/152‑251 (MQ=255)                                                                                               cgacgccggaaggtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTCATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCa                                                                                                                                                 >  2:184796‑M2/149‑250 (MQ=255)                                                                                                    ccggaaggtttgatcacaaagacgccgctctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggatgtgttttggcagcgccaggttggcttctaatacccggcgcccgagaaccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAGGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAACTGCTATGGCATCGCAACGGGTGACGCCGTGCAaataat                                                                                                                                           >  2:213896‑M2/144‑251 (MQ=255)                                                                                                           gtttgatcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAAtgtg                                                                                                                                    <  1:269635‑M2/115‑1 (MQ=255)                                                                                                                atcacaaagacgccgcgctcgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACtt                                                                                                                               <  1:94767‑M2/120‑1 (MQ=255)                                                                                                                             cgcgctcgcgatcaacggcgctgacgttgccacatgtgagcgtgaccaggttgtgttttcgcagcgccaggttggcttctaatacccggcgcccgagccccgcaagcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGAtt                                                                                                                  >  1:18481‑M2/119‑251 (MQ=255)                                                                                                                                   cgcgatcaacggcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGa                                                                                                             <  2:96472‑M2/138‑1 (MQ=255)                                                                                                                                              gcgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAc                                                                                                 <  1:199244‑M2/150‑1 (MQ=255)                                                                                                                                               cgctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAcg                                                                                                <  1:5506‑M2/151‑1 (MQ=255)                                                                                                                                                gctgacgttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAcgc                                                                                               >  2:315682‑M2/100‑251 (MQ=255)                                                                                                                                                      gttgcccaatgtgagcgtgaccaggttgggttttggcagcgccaggttggcttataatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGAGttttt                                                                                         >  1:86378‑M2/94‑251 (MQ=255)                                                                                                                                                                                       tggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAAt                                                         <  1:108134‑M2/190‑1 (MQ=255)                                                                                                                                                                                               ccaggttggctgctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTtaat                                                <  1:142156‑M2/199‑1 (MQ=255)                                                                                                                                                                                                          tctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGtt                                      <  1:3322‑M2/209‑1 (MQ=255)                                                                                                                                                                                                             aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGCGAAAAAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGCCGCCGTGCAAATAATCAATGTGGCCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg                                  >  1:344577‑M2/39‑251 (MQ=255)                                                                                                                                                                                                                          cgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCg                      >  2:97064‑M2/26‑250 (MQ=255)                                                                                                                                                                                                                                             cctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACg  <  1:186470‑M2/245‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACG  >  NC_000913/69831‑70317

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