Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,279,023	T→C	F182S (TTC→TCC)	ghxP →	guanine/hypoxanthine transporter GhxP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,023	0	T	C	98.6%	246.0 / NA	70	F182S (TTC→TCC)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (35/34);  minor base G (0/1);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGC  >  NC_000913/4278777‑4279268                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      tGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGAGACGTGTTCCTGATGGGTGTGCTGTTTCCGGCAATTTCTGCCACGGGTATAAGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTACCGGTTTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCg                                                                                                                                                                                                                                                   >  2:39050/1‑251 (MQ=255)  gCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCgg                                                                                                                                                                                                                                                  >  1:69137/1‑251 (MQ=255)   cTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGt                                                                                                                                                                                                                                                 >  1:615547/1‑251 (MQ=255)     ggggCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTg                                                                                                                                                                                                                                                <  2:361451/250‑1 (MQ=255)      gggCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAt                                                                                                                                                                                                                                              <  1:38990/251‑1 (MQ=255)            catgttAGCGTTCCGGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTGCTGCCACGGGGATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGTTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGCGTCGGTCTGGTGATTAAAACCCCGCTTGATAGTCTGCCCTTTGCGCTGTGTGATTTCGCGGCCTCCCCGGTGAttctcac                                                                                                                                                                                                                                        <  2:423251/247‑6 (MQ=255)                ttAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGt                                                                                                                                                                                                                                  >  1:716751/1‑251 (MQ=255)                    cGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGt                                                                                                                                                                                                                              <  1:616735/251‑1 (MQ=255)                           gTCGGGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatac                                                                                                                                                                                                                         <  1:31958‑M1/249‑5 (MQ=255)                              gcgcTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAATCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaac                                                                                                                                                                                                                    <  2:536133‑M1/251‑10 (MQ=255)                                   gggTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaa                                                                                                                                                                                                               <  2:216638‑M1/251‑15 (MQ=255)                                      tGCCGTGTTCCTGATGGGTGTGCTGTTTACGTTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGTGATTGTTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGGCTGGTGATTAAAAACCCGCTTGATGGTCGGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGTTGAT‑‑TGTCACTGGTAg                                                                                                                                                                                                                                <  1:37641/231‑1 (MQ=255)                                      tGCCGTGTACCTGGTGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTCTAAGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGAACGGGGATTGGAATAGGTCTGTTCCTGCTGCTCATTGCCGATAACGGTGTCGGTCTGGTGATTAAACACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCAGGTGAT‑‑TGTCACTGGTAg                                                                                                                                                                                                                                >  2:37641/1‑231 (MQ=255)                                      ggcCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacag                                                                                                                                                                                                            <  2:69137‑M1/250‑18 (MQ=255)                                       gCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTACCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGA‑‑ATGTCACTGGTAGGt                                                                                                                                                                                                                              >  2:420974/1‑232 (MQ=255)                                       gCCGTGTTCCTGATGGGTGTGCTGTTGACGGTAATTTCTGCCTCGGGTATCCGTAGCTGGATTTTGCGCAACTGGCCTCACGGTGTGGCGCACGGCCCGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGt                                                                                                                                                                                                                              <  1:420974/232‑1 (MQ=255)                                        ccGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacagga                                                                                                                                                                                                          <  1:731739‑M1/251‑20 (MQ=255)                                                  gATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaac                                                                                                                                                                                                  <  1:144839‑M1/249‑28 (MQ=255)                                                   aTGGGTGTGCTTTTTACGTTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCattcctaacgacaacaggaaaaaaaaccat                                                                                                                                                                                               <  1:568994‑M1/251‑31 (MQ=255)                                                   aTGGGTGTGCTTTTTACGGTAATGTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCGCGGTGTGGCGCACGGCACGGGGAGTGGTATCGGTCTGTTCCTGCTGCTCAGTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGTTGATTTCGCGACCTCCCCGGGGA‑‑GTGTCACTGGTAGGTCatacctaacgtcaacaggaaaaaaaaccat                                                                                                                                                                                               <  2:363212‑M1/251‑31 (MQ=255)                                                     gggTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatac                                                                                                                                                                                             <  2:237064‑M1/251‑33 (MQ=255)                                                      ggTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAGACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacc                                                                                                                                                                                            >  2:69351‑M1/1‑218 (MQ=255)                                                              gattacGGTAATTTCGGCCACGGGTATCCGTAGCTGGAGATTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTGGATGGTCTGCCCGTTGCGCGGGGGGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCataccgaacgacaacaggaaaaaaaaccatacccaaaaggt                                                                                                                                                                                    <  2:733100‑M1/249‑42 (MQ=255)                                                              gTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggt                                                                                                                                                                                    <  2:459175‑M1/251‑42 (MQ=255)                                                                tgacgGTAATTGCTCCCACGGGTAGCCGTAGCTGGATTTTGCGCAGCTGGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCGGCTCAGTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGGTGCGCTGGGTGATGTCGCGACCTGCCCTGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccagacccaaaaggta                                                                                                                                                                                   <  2:226260‑M1/248‑43 (MQ=255)                                                                  aCGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtg                                                                                                                                                                                <  1:321492‑M1/251‑46 (MQ=255)                                                                   cGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCCCGGCACGGGGCTTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtg                                                                                                                                                                                >  1:303050‑M1/1‑205 (MQ=255)                                                                    ggTAATTTCTGCCACGGGTATCCGTAGCTGGAGTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGTTTGGTATCGGTCTGTTCCTGCTGCTCATGGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATGTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatgcctaacgacaacagtaaaaaaatccatacccaaaaggtagtgaa                                                                                                                                                                              <  2:354885‑M1/251‑48 (MQ=255)                                                                           tCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGGTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgttt                                                                                                                                                                        <  1:322551‑M1/250‑54 (MQ=255)                                                                                   gggTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaac                                                                                                                                                               >  2:218338‑M1/1‑189 (MQ=255)                                                                                      tATCCGTAGCTGGATTTTGCGCAACTTGACTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaacggtagtgaatcgtttgcttttaactgg                                                                                                                                                            >  1:603347‑M1/1‑186 (MQ=255)                                                                                       aTCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccataccaaaaaggtagtgaagcgtttgcttttaactggg                                                                                                                                                           >  2:262745‑M1/1‑185 (MQ=255)                                                                                             aGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgt                                                                                                                                                     >  1:129989‑M1/1‑179 (MQ=255)                                                                                                 ggATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaagggagtgaatcgtttgcttttaactggggattgttca                                                                                                                                                  <  2:623580‑M1/250‑76 (MQ=255)                                                                                                    ttttGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtg                                                                                                                                              >  1:354464‑M1/1‑172 (MQ=255)                                                                                                      ttGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtgg                                                                                                                                             <  2:603347‑M1/250‑81 (MQ=255)                                                                                                       tGCGCAACTTGCCTCACGGTGTGGCGCTCGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggc                                                                                                                                            <  2:293600‑M1/250‑82 (MQ=255)                                                                                                        gcgcAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcg                                                                                                                                           <  2:716751‑M1/250‑83 (MQ=255)                                                                                                            aaCTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcg                                                                                                                                      >  1:450518‑M1/1‑164 (MQ=255)                                                                                                            aaCTTGCCTCACGGTGGGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcg                                                                                                                                      <  2:306033‑M1/251‑88 (MQ=255)                                                                                                                       cGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACTGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacatgaaaaaaaaccatacccaaaaggtagtgaagcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaat                                                                                                                           <  1:649845‑M1/251‑99 (MQ=255)                                                                                                                             ggCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgtta                                                                                                                       >  2:126059‑M1/1‑147 (MQ=255)                                                                                                                             ggCGCACGGCACGGGGATTGGGATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaag                                                                                                                     <  2:172620‑M1/251‑105 (MQ=255)                                                                                                                              gcgcACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGA‑‑GTGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaaga                                                                                                                    <  1:366794‑M1/251‑106 (MQ=255)                                                                                                                                gcacggcacgGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAGACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagact                                                                                                                  <  1:69351‑M1/251‑108 (MQ=255)                                                                                                                                       acggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggag                                                                                                           <  1:594938‑M1/251‑115 (MQ=255)                                                                                                                                       acggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTACCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaaaaggaaaaacaaccatacccaaaaggtagtgaatagtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttacgactgtggg                                                                                                             >  1:46099‑M1/1‑137 (MQ=255)                                                                                                                                           ggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggaggg                                                                                                         >  1:115598‑M1/1‑133 (MQ=255)                                                                                                                                            gATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggaggg                                                                                                         >  2:539230‑M1/1‑132 (MQ=255)                                                                                                                                                ggTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatct                                                                                                  <  1:250242‑M1/251‑124 (MQ=255)                                                                                                                                                               ctgctgCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCAGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaccaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcc                                                                                   >  1:662863‑M1/1‑113 (MQ=255)                                                                                                                                                                    gctCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctgg                                                                              >  1:309586‑M1/1‑108 (MQ=255)                                                                                                                                                                       cATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGAGATTTCGCGCCCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccacaaggtagtgaatcgtttgcttttaactggggattgttcaggtggagttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcagggaccgcctctgggat                                                                           >  2:68622‑M1/1‑105 (MQ=255)                                                                                                                                                                          tGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggcgaatctgtggcaggaaccgcctctggtatagg                                                                        >  1:116257‑M1/1‑102 (MQ=255)                                                                                                                                                                                                      aaaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGAAATCACCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggcagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattattttctct                                            >  1:105019‑M1/1‑74 (MQ=255)                                                                                                                                                                                                       aaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtg                                                                                                                                              >  2:214552‑M1/1‑73 (MQ=255)                                                                                                                                                                                                       aaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtg                                                                                                                                              <  1:214552‑M1/152‑80 (MQ=255)                                                                                                                                                                                                            ccGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctct                                      >  1:304040‑M1/1‑68 (MQ=255)                                                                                                                                                                                                              gCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttc                                    >  1:154192‑M1/1‑66 (MQ=255)                                                                                                                                                                                                                  gATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTAataccaaacgaaaacagaaaaaaaaaccatacccaaaaggtagtgaatcttttgtttttaacgggggattgttcgggtggcgttgcggtgttttgaatgttaagattgggggggggaaattctggggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcg                                <  1:539230‑M1/251‑190 (MQ=255)                                                                                                                                                                                                                   aTGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgc                               >  1:249562‑M1/1‑61 (MQ=255)                                                                                                                                                                                                                        cTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCCATGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcccctggtataggggggggcgaagatagccttatttcctctccctcttcttcgctgtac                          >  2:664891‑M1/1‑56 (MQ=255)                                                                                                                                                                                                                        cTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaacc                                                                                      <  1:487528‑M1/191‑136 (MQ=255)                                                                                                                                                                                                                        cTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCACGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaacaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaacc                                                                                      >  2:487528‑M1/1‑56 (MQ=255)                                                                                                                                                                                                                            ccGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTAatacctaacgacaacaggacaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcagggaccgcctctggtataggggggggcggagatagcattatttcctctccctcttcttcgctgtacc                         >  2:731550‑M1/1‑52 (MQ=255)                                                                                                                                                                                                                                tGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTAACTGGTCGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataac                  >  1:51685‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                   gcTGGGAGATATCGCGACCACCACGGTGAT‑‑TGTCCCTGGTAGGTCatccctaacgccaacaggacacaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtataggggggggcgaagaacgcattatttcctctccctctcttcgctgtgccacataaccggc               >  2:89010‑M1/1‑45 (MQ=255)                                                                                                                                                                                                                                    cTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagctaaccggc              >  2:562678‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                                              tCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaacgcccctggtatagggggtggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcacccttaatt    >  1:477547‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                                                gcgACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcaccattaatttt  >  1:410076‑M1/1‑32 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      TGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGC  >  NC_000913/4278777‑4279268

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