Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,222	A→C	I156L (ATC→CTC)	rhtC →	threonine efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,222	0	A	C	94.6%	186.5 / ‑6.0	57	I156L (ATC→CTC)	rhtC	threonine efflux pump
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (31/22);  minor base G (3/0);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCA  >  NC_000913/4007990‑4008441                                                                                                                                                                                                                                          |                                                                                                                                                                                                                            tgaccttgacctctGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTc                                                                                                                                                                                                                           <  2:13101/225‑1 (MQ=255)                ctgGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTgg                                                                                                                                                                                                           <  1:42444/236‑1 (MQ=255)                  ggATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCg                                                                                                                                                                                                         <  1:229564/236‑1 (MQ=255)                        gTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTTTGCACCTGCGCCACAGGGCGGGCTGGCGAAAAGTGGG‑GGCAGTTTCCTGAAAGGGTTACTGACCAATCTCGCCAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTTTCGGTGGTAACGTTTGCACTACCCCGCGCTGGGGCGCTTTTTCGCTGCTCCTTGTCGAAActcttgactggctt                                                                                                                                                                                                   >  2:268712/1‑226 (MQ=255)                               gATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGGGAAAAAGGGG‑GGCAGTTTCCTGAAAGGGTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGGGTTCTCATTGTTTGTCGGGGAGAACGTTGGCACTACCGCGCGCTGGGGCCTTTTTGTGCTGGCCATTGTCGAAACGCTGGCGTGGGGTACcgtc                                                                                                                                                                                             >  2:165894/1‑235 (MQ=255)                                  gCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGGCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTTTCATTGTTTGTCGGGGATAAAGTTTGCACTACCGCGCGCTGGGGGCTTTTTTTGCTGGCCATTGTCGAAACGCTGGCGTGGTGTACCGTCGttt                                                                                                                                                                                         >  2:578748/1‑235 (MQ=255)                                                     aaaGAGGCGGTTTTTTCACCTGCGCCACAGGGCGAGCTGGCGAAAAGTGGG‑CGCCGTTTTCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTTTCGGGGATAACGTTGGCACTACCCCGCGCTTGGGGATTTTTGCGCTGTGCCTTGTCGAAACGCCGGGGTGGTTTTCCGTCGTTGTCCGCCTCTTTTGCCTgtc                                                                                                                                                                      >  2:22080/1‑234 (MQ=255)                                                            cGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAATGGG‑CGCCGGTTCCTGAAAGGGTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTTGCTCGGGGTTCTCATTGTTTGTCGGGGATAACGTTTGCACTACCGCGCGCTGGGGCATTTTTTCGCTGGTCCTTTTCGAAACGCTGGCGTGGGTTAACGTCGTTGCCAGGCTTCTTGCCCTGGCGCcaaa                                                                                                                                                                >  2:694997/1‑234 (MQ=255)                                                               tttCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcg                                                                                                                                                             <  1:165894/235‑1 (MQ=255)                                                               tttCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcg                                                                                                                                                             <  1:483826/235‑1 (MQ=255)                                                                   tGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTg                                                                                                                                                        <  1:694997/236‑1 (MQ=255)                                                                         tGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATc                                                                                                                                                  >  1:237275/1‑236 (MQ=255)                                                                          gcgcCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCa                                                                                                                                                 >  1:686172/1‑236 (MQ=255)                                                                           cgcCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCa                                                                                                                                                 >  1:449702/1‑235 (MQ=255)                                                                                  ggTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTgg                                                                                                                                         >  1:142960/1‑236 (MQ=255)                                                                                       aGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGaa                                                                                                                                     <  1:448592/235‑1 (MQ=255)                                                                                            gccaaaaATGGGCCGAAGTTCCCTAAAGGGTTTACTGCCCAATCCCGCTATTCCGAAGGCGATTATCTACTTGGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTCCCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTgg                                                                                                                                 <  2:24657/232‑1 (MQ=255)                                                                                             cGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTGGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAtt                                                                                                                              >  1:662967/1‑236 (MQ=255)                                                                                              gAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTg                                                                                                                             <  1:536218/236‑1 (MQ=255)                                                                                                              tttCCTGAAAGGTTTACTGGCCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCgggg                                                                                                              <  1:268712/236‑1 (MQ=255)                                                                                                                  cTGAAAGGTTTACTGACCAATCCCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGtt                                                                                                          <  1:473397/236‑1 (MQ=255)                                                                                                                        ggTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCCGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGc                                                                                                    >  1:691551/1‑236 (MQ=255)                                                                                                                          tttACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCg                                                                                                  >  1:593845/1‑236 (MQ=255)                                                                                                                                    aaTCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCGTTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCAGGATTTGGcat                                                                                        >  1:742396/1‑236 (MQ=255)                                                                                                                                       ctcGCTATTCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTTCCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGcattca                                                                                     <  1:484675/236‑1 (MQ=255)                                                                                                                                       ctcGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGcattca                                                                                     <  1:22080/236‑1 (MQ=255)                                                                                                                                                         aTTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCTCGGt                                                                    >  2:551305/1‑235 (MQ=255)                                                                                                                                                         aTTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGGTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGt                                                                    >  1:246410/1‑235 (MQ=255)                                                                                                                                                               tACTTTGGCTCGGTTTTCTCATTGTTTGTCGGTGATCACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCa                                                             <  2:41960/236‑1 (MQ=255)                                                                                                                                                                 cTTTGGCTGGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGa                                                           <  1:578748/236‑1 (MQ=255)                                                                                                                                                                  tttGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAATTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGAc                                                          >  1:64981/1‑236 (MQ=255)                                                                                                                                                                   gtGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGAcg                                                         <  1:219121/235‑1 (MQ=255)                                                                                                                                                                    tGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGAc                                                          >  1:281491/1‑234 (MQ=255)                                                                                                                                                                       cTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCCGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTc                                                     >  1:560593/1‑236 (MQ=255)                                                                                                                                                                           gtgtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGa                                                           >  1:634855/1‑226 (MQ=255)                                                                                                                                                                           gtgtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGa                                                           <  2:634855/226‑1 (MQ=255)                                                                                                                                                                           gtgtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCa                                                 >  2:565419/1‑236 (MQ=255)                                                                                                                                                                                ctcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTa                                            >  1:189858/1‑236 (MQ=255)                                                                                                                                                                                   aTTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGa                                                                                                                                      <  2:409398/143‑1 (MQ=255)                                                                                                                                                                                   aTTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGa                                                                                                                                      >  1:409398/1‑143 (MQ=255)                                                                                                                                                                                      gtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGGGATGCCAGACGCGTCTTCAGAGTAAGTCgg                                      >  1:660308/1‑236 (MQ=255)                                                                                                                                                                                        ttgtCGGTTATAACGTTGGCACTACGGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGAt                                    <  1:670570/236‑1 (MQ=255)                                                                                                                                                                                        ttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTTCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAACTCggtg                                    >  2:693111/1‑234 (MQ=255)                                                                                                                                                                                         tgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGTCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGAt                                    >  2:391292/1‑235 (MQ=255)                                                                                                                                                                                         tgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGGTTATTTCGCGGGGATGCCAGACGGGGCTTCATAGTAAGTCGGTTa                                   >  2:224848/1‑236 (MQ=255)                                                                                                                                                                                          gcgggTGATAACGTTGGCACTACCGCGCGCTGGGGCATGTTTGCGCTGCTCATTGTCGAAACGCTGTCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATaa                                  <  2:449702/233‑1 (MQ=255)                                                                                                                                                                                               tGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTACGAGTAATTCGGATAAGGCg                              >  2:544307/1‑235 (MQ=255)                                                                                                                                                                                                 aTAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGtt                            >  1:133306/1‑235 (MQ=255)                                                                                                                                                                                                        tGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAcgccg                     >  1:713879/1‑235 (MQ=255)                                                                                                                                                                                                         ggCACTACCGCGCGCCGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGCGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGCAGTGGTTTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCa                                                 >  1:608728/1‑206 (MQ=255)                                                                                                                                                                                                           cACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCAt                  <  1:565419/235‑1 (MQ=255)                                                                                                                                                                                                                ccGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACa            <  1:224848/236‑1 (MQ=255)                                                                                                                                                                                                                    gcgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGGGATGCCAGACGCGTCTTCAGAGTAAGTCGGCTAAGGAGCTTACGCCGCATCTGCCattgt        >  2:393552/1‑234 (MQ=255)                                                                                                                                                                                                                        tGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAtttt     >  2:186182/1‑235 (MQ=255)                                                                                                                                                                                                                          gggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCCGCGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAtttttta  >  2:198041/1‑234 (MQ=255)                                                                                                                                                                                                                          gggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa  >  1:217000/1‑236 (MQ=255)                                                                                                                                                                                                                          gggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa  >  1:463781/1‑236 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                            TGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGCGAAAAGTGGG‑CGCAGTTTCCTGAAAGGTTTACTGACCAATCTCGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCA  >  NC_000913/4007990‑4008441

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