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

GGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCA  >  NC_000913/66287‑66772                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 ggggAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccg                                                                                                                                                                                                                                              >  1:185733‑M1/1‑246 (MQ=255)   ggAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTAAGCGGTCATGACGATGTAAACGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccg                                                                                                                                                                                                                                              >  1:238031‑M1/1‑244 (MQ=255)    gAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagac                                                                                                                                                                                                                                          >  1:60971‑M1/1‑243 (MQ=255)     aaTGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacc                                                                                                                                                                                                                                         <  1:192646‑M1/251‑10 (MQ=255)        gCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgc                                                                                                                                                                                                                                      <  2:185733‑M1/251‑13 (MQ=255)                                ggCGTGTCGGAGGAGGGCTTTTTCGGTCCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATAGCGTCAGCGGGCATGACGCTGTAATCGACGCCGGAAGGTTTGAGCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCTTCGcccgatacccgcacgcgtcggcctgcaggcatggag                                                                                                                                                                                                              <  2:315063‑M1/251‑37 (MQ=255)                                       cGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTTTAATACCGGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                                                                                                                                                                                       <  1:127448‑M1/251‑44 (MQ=255)                                             aggGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGAGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagag                                                                                                                                                                                                   >  1:42990‑M1/1‑202 (MQ=255)                                                     ttCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgat                                                                                                                                                                                          >  1:59475‑M1/1‑194 (MQ=255)                                                        gTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaa                                                                                                                                                                                      <  1:81836‑M1/251‑61 (MQ=255)                                                          accctcaacCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCTGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTGTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcatgtatggataaacagtagagagttgcgataaaaa                                                                                                                                                                                     <  2:308514‑M1/246‑62 (MQ=255)                                                           cctgcaaccACCTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGAGGTCAGGACGCTGTAAGCGACGCCGGAAGGTGTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagatccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagc                                                                                                                                                                                   <  2:270865‑M1/247‑64 (MQ=255)                                                            cTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcg                                                                                                                                                                                  >  2:58563‑M1/1‑187 (MQ=255)                                                                 aCCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcagg                                                                                                                                                                             <  2:60971‑M1/251‑70 (MQ=255)                                                                      ttCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtaaggtagg                                                                                                                                                                         >  2:293806‑M1/1‑177 (MQ=255)                                                                         aCCGGTTGCGATGGTAACCGCGACCATAGCGTCAGCGGGCATGACGCTGTAATCGACGCCGGAAGGTTTGATCGCAAAGCCGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGGTTCTAATACCCGGCGcccgagacccgcacgcgtcgaccggcaggtatggagaaacagtatagagttgcgataaaaagcgtcagttaggagcc                                                                                                                                                                     <  2:247033‑M1/251‑78 (MQ=255)                                                                                     gCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTGTGAGCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTGCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagatagttgcgataaaaagcgtcaggtagtatccgctaatcttat                                                                                                                                                          <  2:43898‑M1/250‑89 (MQ=255)                                                                                                    atCGTCAGCGGTCATGTCGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatg                                                                                                                                          <  2:164138‑M1/251‑105 (MQ=255)                                                                                                    atCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagcgagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatg                                                                                                                                          >  1:286885‑M1/1‑147 (MQ=255)                                                                                                           gCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagc                                                                                                                                   <  1:362482‑M1/251‑112 (MQ=255)                                                                                                             ggTCAGGACGCTGTAATCGACGCCGGAAGGTTTGAGCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGTCGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaa                                                                                                                                 <  1:239972‑M1/251‑114 (MQ=255)                                                                                                                cATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCGGACGTTGCCCCATGTGAGCGGGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagt                                                                                                                              <  2:272341‑M1/251‑117 (MQ=255)                                                                                                                           aaTCGACGCCGGCAGGTTTGATCACAAAGACGCCGCGCTCGCGCTCAAAGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagcaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtg                                                                                                                   >  2:213075‑M1/1‑124 (MQ=255)                                                                                                                              cGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctcatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgca                                                                                                                 >  2:184796‑M1/1‑121 (MQ=255)                                                                                                                                   ccGGAAGGTTTGATCACAAAGACGCCGCTCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGATGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagaaccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaggcgtcaggtaggatccgctaatcttatggataaaactgctatggcatcgcaacgggtgacgccgtgcaaataat                                                                                                           >  2:213896‑M1/1‑116 (MQ=255)                                                                                                                                          gTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtg                                                                                                    <  1:269635‑M1/251‑143 (MQ=255)                                                                                                                                               aTCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggactt                                                                                               <  1:94767‑M1/251‑148 (MQ=255)                                                                                                                                                            cgcgCTCGCGATCAACGGCGCTGACGTTGCCACATGTGAGCGTGACCAGGTTGTGTTTTCGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagccccgcaagcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgatt                                                                                  >  1:18481‑M1/1‑91 (MQ=255)                                                                                                                                                                  cgcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattataga                                                                             <  2:96472‑M1/250‑166 (MQ=255)                                                                                                                                                                             gcgcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttac                                                                 <  1:199244‑M1/251‑178 (MQ=255)                                                                                                                                                                              cgcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacg                                                                <  1:5506‑M1/251‑179 (MQ=255)                                                                                                                                                                               gcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgc                                                               >  2:315682‑M1/1‑72 (MQ=255)                                                                                                                                                                                     gTTGCCCAATGTGAGCGTGACCAGGTTGGGTTTTGGCAGCGCCAGGTTGGCTTATAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgagttttt                                                         >  1:86378‑M1/1‑66 (MQ=255)                                                                                                                                                                                                                      tGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaat                         <  1:108134‑M1/250‑218 (MQ=255)                                                                                                                                                                                                                              ccAGGTTGGCTGCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaatgcttttaat                <  1:142156‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                                         tCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgacgccgtgcaaataatcaatgtggacttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaatgcttttaataagcggggtt      <  1:3322‑M1/250‑237 (MQ=255)                                                                                                                                                                                                                                            aaTACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggcgaaaaagtagagagttgcgataaaaagcgtcaggtaggatccgctaatcttatggataaaaatgctatggcatagcaaagtgtgccgccgtgcaaataatcaatgtggccttttctgccgtgattatagacacttttgttacgcgtttttgtcatggctttggtcccgctttgttacagaatgcttttaataagcggggttaccg  >  1:344577‑M1/1‑11 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 GGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCA  >  NC_000913/66287‑66772

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: ‑