Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,120,786	G→A	T165I (ACT→ATT)	glcA ←	glycolate/lactate:H(+) symporter GlcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,120,786	0	G	A	98.1%	159.8 / NA	53	T165I (ACT→ATT)	glcA	glycolate/lactate:H(+) symporter GlcA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/29);  minor base C (1/0);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACCAATATAGTTAGAGGTAAAGAACTGAGTGACAGCGAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGC  >  NC_000913/3120540‑3121024                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               ggACCAATATAGTTAGAGGTAAAGAACTGAGTGACAGCGAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATAtt                                                                                                                                                                                                                                            >  1:345854/1‑251 (MQ=255)      aaTATAGTTAGAGGTAAAGAACTGAGTGACAGCGAAGCTTCCCCCACCAACCAGCGCCGCTGGCCCCGTCTCTTCCACCCCTTTCCAGCCCTCCATCCTTGCTCACAGCCCGACCGGCA‑CCAGAACCGCCAGGACCGGTAACCGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCACTATTGGCa                                                                                                                                                                                                                                        >  1:432113/1‑250 (MQ=255)       atatAGTTAGAGGTAAAGAACTGAGTGACAGCGACGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAAt                                                                                                                                                                                                                                      <  2:13363/251‑1 (MQ=255)           aGTTAGAGGTAAAGAACTGAGTGACAGCGAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGTACCAAACGCCACCGGCGCAATATTGGCAATCAg                                                                                                                                                                                                                                   <  1:165642/250‑1 (MQ=255)               tgagGTAAAGAACTGAGTGACAGCGAAGCTCTCCTCAGAAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCATCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGAACGGCGATCAGAATCGGCCTGCCCAACTCACCAATCGCCACCGACGCAATCTTGGCAATCAGacac                                                                                                                                                                                                                               >  2:111454/2‑250 (MQ=255)               agagGTAAAGAACTGAGTGACAGCGAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGAcaca                                                                                                                                                                                                                              <  1:138217/251‑1 (MQ=255)                     aaaGAACTGAGTGACAGCGAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCAAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCCACGCACCCAACGCCACCGGCGCAATATTGGAAATCAGACACAGCCccg                                                                                                                                                                                                                        >  1:275041/1‑251 (MQ=255)                                 gACAGCGAGGCTCCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTGCCCCCCTTTCCATCCGCCCTGCATTTCTACCAGCCGGAACGGCA‑GAATACCCGACAGGAACTGTAACTGACGGCCCGCCAGTGCGCCAATGAGGAACGGATCGATTCCTGTTACCTGACCGGCGTCCAGAAGCGGCCCGCCCAACGCACCAAACGGGACCGGCGCAATATTGGCATGGAGGCACAGCGCCGCTGCGTATAACg                                                                                                                                                                                                            <  2:194977/251‑1 (MQ=255)                                     gCGAAGCTTCCCCCATCACCCAGCGCCGCTGGCCACGTCTCATTCCCCCCTTTCCAGCCGTCCATCAGTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCCCCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCTCGTATAACGGttt                                                                                                                                                                                                        <  2:326249/251‑1 (MQ=255)                                       gacgCTTCCCCCAGCAACCAGCGCCGCTGGCCACTTCTCTTTCACCCCTTCCCAGCCGTCCATCATTGCTCCCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGa                                                                                                                                                                                                      <  2:258425/248‑1 (MQ=255)                                       gAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTg                                                                                                                                                                                                       >  2:69181/1‑250 (MQ=255)                                       gAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTACCGCTACCTGACCGGCGACCAGACTAGGCACGCCCAACGCACCAAACGCCACCCGCGCAATATTGGCAATCAGACACAGCCACGCCGCGAATAACGGCTTg                                                                                                                                                                                                       >  2:287632/1‑250 (MQ=255)                                         aGCTGCCCCCAGCAACCAGCCCCGCTGGCCACGTCTCTTTCACCCCTGTCCAGCCGACCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGACCGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAg                                                                                                                                                                                                    <  1:46690/251‑1 (MQ=255)                                            ttCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGccc                                                                                                                                                                                                 >  1:370560/1‑251 (MQ=255)                                                cccAGCAGCCAGCGCCCCGGGCGACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATGGCTACCAGCCGGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGc                                                                                                                                                                                             <  2:397068/251‑1 (MQ=255)                                                  cAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGccc                                                                                                                                                                                           >  1:405443/1‑251 (MQ=255)                                                      aaCCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CACGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAAGCccacc                                                                                                                                                                                        >  1:396124/1‑250 (MQ=255)                                                       aCCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCccacc                                                                                                                                                                                        >  1:249713/1‑249 (MQ=255)                                                       aCCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGAAACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCccacc                                                                                                                                                                                        >  1:283785/1‑249 (MQ=255)                                                             gccgcTGGCCACGTCTCTTTCAACCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGACTCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCACCAgcgc                                                                                                                                                                                >  1:123857/1‑251 (MQ=255)                                                                  gggCCACGTCTCTTCCCCCCCTTTCCAGCCGTCAATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATGGCGCCAATGTGGAACGGATCGATGCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCgg                                                                                                                                                                           <  2:208613/250‑1 (MQ=255)                                                                                         tCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCaa                                                                                                                                                      >  1:36904/1‑249 (MQ=255)                                                                                          ccacccGTCCATCTTTGCCCCCAGCCAGAACGGCC‑CAAGAACCGACAGGAACGGTACCTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCTCCAAagc                                                                                                                                                   <  2:405443/247‑1 (MQ=255)                                                                                           cAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAagcc                                                                                                                                                  >  1:116530/1‑251 (MQ=255)                                                                                                      catTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGGGGAACGGATCGATTCCCGTTACCTGACCGGAGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCATCGGTCATCGCCACCGGCGCACCAAAGCCAGACGCTCCtt                                                                                                                                       >  1:224116/1‑251 (MQ=255)                                                                                                       atTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACTTCCCGCCATTGCGCCAATGTGGACCGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGCCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCCCCAAAGCCAGCCGCTCCTTc                                                                                                                                      <  1:136089/251‑1 (MQ=255)                                                                                                          gctcccagccagAACGGCA‑AAAGACCCGACAGGAACGTAAACTGACTCCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAACCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAACCAGCCCCAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCa                                                                                                                                    <  2:370560/246‑1 (MQ=255)                                                                                                             accaccCAGA‑CGGCAACAAAAACCAACAGGAACGGTAACTGCCCTCCCGCCATTGCGCCAAGGTGGAACGGACCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCCACGCACCAACCGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCCCCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCaa                                                                                                                                <  2:290338/251‑1 (MQ=255)                                                                                                                            a‑caAGAACCGACAGGAACGGTAACTGACGCCCCGCCATTGCGCTAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGaa                                                                                                                 <  1:392252/251‑1 (MQ=255)                                                                                                                              aaaGAACCGACAGGACCGGTAGCCGACGTCCCGCCTTTGCGCCAATGTGGAACGGATCGATCCCCGTTCCCTGACCGGCGCCCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTAGAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACTCACCAAAGGAGaaa                                                                                                                <  2:283570/250‑1 (MQ=255)                                                                                                                                        cAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCCGCAACGCACCAAAGGAGAAACCAATCAGt                                                                                                       >  1:395117/1‑250 (MQ=255)                                                                                                                                            aaCGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCCCCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCACCAACGCACCAAAGGAGAAACCAATCAGTAacac                                                                                                  >  1:257504/1‑251 (MQ=255)                                                                                                                                               ggaaaCTGCCGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTCCCTGACCGGCGACCAGAAGCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACCCAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCTCCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTg                                                                                               <  1:100558/248‑1 (MQ=255)                                                                                                                                                     tGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACg                                                                                         >  1:368343/1‑251 (MQ=255)                                                                                                                                                        cgccccGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTGACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTg                                                                                      <  2:249713/248‑1 (MQ=255)                                                                                                                                                          tCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGAt                                                                                    >  1:368833/1‑251 (MQ=255)                                                                                                                                                            ccaccatggcgccatggtcgAACCGATCGATACCCTTTACCTGCCCGGCGACCAGAATCGGCCCGCCCAACCCCCCAACCGCACCCGGCTCAATATTGGCCATCAGACACACCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCCCACCAGTAATCGACAGCGGCGCTCCAAAGGCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCCAACGCTGAtcg                                                                                  <  2:168232/232‑1 (MQ=255)                                                                                                                                                                     cgcCAATGTGGAACGCATCGATTCCCGTTACCGGACCGGCGACCAGAATCGGCACTCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACACCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCTCCGGTAATCGCCACCGGCGCACCAAAGCTAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATg                                                                         <  1:74623/251‑1 (MQ=255)                                                                                                                                                                      gcCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGCGATgg                                                                        >  1:146537/1‑251 (MQ=255)                                                                                                                                                                          aggtgGAACGGACTGATGCCCGTATCCAGACCGGCGACAAGAATCGGCACGCCCAACGCACCAACCGCAACCGGCGCAAGATGGACAACAAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCAACCAGAAGCGCACCGGTAATCGCAACAGGCGAACAAAAGCCAGCCGCTCCTGCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGAt                                                                    <  2:133083/249‑1 (MQ=255)                                                                                                                                                                            gtgGAACTGATCGCTCCCCGTTACCTGACCGGCGACCAGAATCGGCCCGCCCAACGCACCAACCGCCACCGGCGCAATATTGGCAATCAGCCACAGCCCCGCCGCGTATACCGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATaa                                                                  <  2:354013/251‑1 (MQ=255)                                                                                                                                                                               gAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCg                                                               <  1:350198/251‑1 (MQ=255)                                                                                                                                                                                aaCGGATCGAGTCCCGTTACCGGACCGGCTACCAGAATCGCCACGCCCAACGCCCCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGa                                                              <  2:368343/251‑1 (MQ=255)                                                                                                                                                                                      tCGATTCCCGTTACCTGACCGGCGCCCAGAATCGGCACGCCCAACGCAACACACGCCACCGACGCAATATTGGCAATCAGACACAACCCCGCCGCGCATAACGGTTTGAAGCCTAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCCACGCACCAAAGGAGAAACCCATCAGTAACACCTGCAAACGCTGCTCGTCGGTGATGGAGATAAccacgatgcg                                                        >  2:459754/1‑243 (MQ=255)                                                                                                                                                                                        gATTCCCTTTACCGGACCGGCGACCAGAATCGGCACGCCCAACGCCCCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGa                                                      <  1:276603/251‑1 (MQ=255)                                                                                                                                                                                              cGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAatat                                                <  1:177089/251‑1 (MQ=255)                                                                                                                                                                                                   ccggacaggcgacaaaaaccGGCACCCACAACGCTCAAAACGCCACAGGCGCAAAATTGGAAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAACGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCCGCAACGCAACAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGCTAACCGAGCTGCGGATAATATCGAAc                                           <  2:136624/232‑1 (MQ=255)                                                                                                                                                                                                              aCCAGAACCGGCACCCCCAACGCCCCAAACACCCCCGGCTCAATATGGCCAATCAGCCACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCTGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGAGAACCGAGCTGCGGATAAGAGCTAACTGCCCGCTGGc                                <  2:283785/251‑1 (MQ=255)                                                                                                                                                                                                                 aGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAACCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAACGCCCGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAACCGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAAc                             >  1:267047/1‑251 (MQ=255)                                                                                                                                                                                                                             cccAACGCCCCAAACGCCACCGGCGCAATATTGGCAACCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAGCAGTAACACCTGCAAACGCGGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTtata                 <  1:287632/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaaCGCCCCCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACAccg      <  1:138560/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaCGCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGCTAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACAccgc     >  1:422472/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCCACCGGCGCAATATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCAACAAAGCCAGCCGCTCCTTCCATCAACGCACCACAGGAGAAACCAATCAATAACACCTGCAAACGCTGATAGTCGGTGATGGAGATAACCACGCTGCGGATAATATCGCACTGCCCGCTTGCAACGGTTAATTTATACAGGAACACctccgc  >  1:298154/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GGACCAATATAGTTAGAGGTAAAGAACTGAGTGACAGCGAAGCTTCCCCCAGCAACCAGCGCCGCTGGCCACGTCTCTTTCACCCCTTTCCAGCCGTCCATCATTGCTACCAGCCAGAACGGCA‑CAAGAACCGACAGGAACGGTAACTGACGTCCCGCCATTGCGCCAATGTGGAACGGATCGATTCCCGTTACCTGACCGGCGACCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGC  >  NC_000913/3120540‑3121024

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