Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,089,260	T→C	V326A (GTG→GCG)	galP →	galactose:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,089,260	0	T	C	96.9%	220.5 / NA	65	V326A (GTG→GCG)	galP	galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/36);  minor base A (2/0);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGC  >  NC_000913/3089014‑3089489                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      agaACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATg                                                                                                                                                                                                                                   >  1:535961/1‑251 (MQ=255)    accgCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACGTTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGTGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTTGGGACGTAAACCAACGCTAACGCTGGTCTTCCTGGCGATGgct                                                                                                                                                                                                                                <  1:449001/248‑1 (MQ=255)        cAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGctgctg                                                                                                                                                                                                                            <  2:482511/251‑1 (MQ=255)        cAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGctgctg                                                                                                                                                                                                                            <  1:3427/251‑1 (MQ=255)         aaCTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGCGCAAATGTGGGGGACCGTGATTGTCGGCCTGAACAACGTACTTTCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGGTGGGCTTCCTTGCGTTGGCTGCTgg                                                                                                                                                                                                                           >  2:118685/1‑251 (MQ=255)           cTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGGGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTAATTGCCACCTTTATCGCAATCGGCCTTGTTTACCGCTGGGGACGTAAACCAACGCTAACGACGGGCTTCCTGGAGATGGCTGCTGGcc                                                                                                                                                                                                                         >  2:60837/1‑250 (MQ=255)              ccgccgCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCAGCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCTATGGCTGCTGGCATgg                                                                                                                                                                                                                      <  2:550854/251‑1 (MQ=255)               cgccgcGTGGTGTGCCGTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATGCTGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTAACGAGCAAATGTGGTGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAAGCGTGCTTGTTGACCGCTGGGGACGTAATCCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATggg                                                                                                                                                                                                                     <  1:568001/251‑1 (MQ=255)               cgccgcGCGGTGTGCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATggg                                                                                                                                                                                                                     <  2:395574/251‑1 (MQ=255)                gccgcGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACGGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAAGCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGc                                                                                                                                                                                                                    <  2:575573/251‑1 (MQ=255)                     gcgGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACt                                                                                                                                                                                                               >  2:145269/1‑251 (MQ=255)                      cgGTGTTCCTTGGCGTGCTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACTTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTc                                                                                                                                                                                                              <  1:520297/251‑1 (MQ=255)                            tCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTAc                                                                                                                                                                                                         <  1:54986/250‑1 (MQ=255)                                   cGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATATTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGc                                                                                                                                                                                                 >  1:506374/1‑251 (MQ=255)                                    gtgcTGTTGCAGGGAATGCAGCAATTCACCGGTATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGTCCTTGTTGACCTCTGGGGACGTAAACCAACGCTAACGCTGGGCTGCCTGGCGATGGCTGCGGGCAAGGGCGTACTCGGTACAATTATGCa                                                                                                                                                                                                <  2:545084/248‑1 (MQ=255)                                            gCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGGGATTGTCGGCCTTACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGt                                                                                                                                                                                         >  2:458301/1‑250 (MQ=255)                                               ggTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTc                                                                                                                                                                                     <  2:16410/251‑1 (MQ=255)                                                   aTGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCActc                                                                                                                                                                                 <  1:218099/251‑1 (MQ=255)                                                     gcagcaATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCActct                                                                                                                                                                                <  1:288625/250‑1 (MQ=255)                                                         aaatactccGGGATGAACGTCAGCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAAGGTGGGGGACCGTGATTGTCGGCCTGCCCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACTCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACGCGGTACAATGATGCATATCGGTATTCACTCTCCGTc                                                                                                                                                                           <  1:312408/244‑1 (MQ=255)                                                            ttCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCCACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGACATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCCGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCCCTCCCCGTCGgc                                                                                                                                                                        >  2:218388/1‑251 (MQ=255)                                                             tCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCTGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCGTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCGGGCATGGGCGTGCTCGGTACAATGATGCATGTCGGTATTCACTCTCCGTCGgc                                                                                                                                                                        <  2:539402/250‑1 (MQ=255)                                                              cACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGgcgc                                                                                                                                                                      <  1:187112/251‑1 (MQ=255)                                                              cACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGgcg                                                                                                                                                                       <  1:378599/250‑1 (MQ=255)                                                                         aCGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTtcgcc                                                                                                                                                           <  1:463390/251‑1 (MQ=255)                                                                             catcatGTATTACGCGCCGAAAATCGTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTtcgccatc                                                                                                                                                        <  1:111053/250‑1 (MQ=255)                                                                                   gTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCAtgc                                                                                                                                                 <  1:145269/251‑1 (MQ=255)                                                                                   gTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTGGTTGACCGCTGGGGACGTAAAGCAACGCTAACGCTGGGCTTCCTGGCGATGGCTTCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATCTCGCCATCGCCAtg                                                                                                                                                  <  1:40215/250‑1 (MQ=255)                                                                                   gTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGAACGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCAtgc                                                                                                                                                 >  1:121759/1‑251 (MQ=255)                                                                                        aCGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATgctgct                                                                                                                                             <  2:29802/250‑1 (MQ=255)                                                                                        aCGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATgctgct                                                                                                                                             <  1:126761/250‑1 (MQ=255)                                                                                                       tCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCCCCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGt                                                                                                                                 <  2:506374/247‑1 (MQ=255)                                                                                                       tCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCg                                                                                                                               <  1:520263/249‑1 (MQ=255)                                                                                                        cGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTTATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGt                                                                                                                             <  1:4895/250‑1 (MQ=255)                                                                                                        cGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTAGCCATCGCCATGCTGCTGATGGTTATTGTCGGtt                                                                                                                            >  1:318503/1‑251 (MQ=255)                                                                                                                cGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGACACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTACTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCAGTCGGAGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATg                                                                                                                    >  1:65133/1‑251 (MQ=255)                                                                                                                  ggTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATg                                                                                                                    >  1:333237/1‑249 (MQ=255)                                                                                                                   gTTACACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGt                                                                                                                 >  1:79646/1‑251 (MQ=255)                                                                                                                      ataCCAACACTACCGAGCACATGTGGGGGCCAGTGATTGTCGGCATGACCACAGTACTTGACACCTTTATCGCAATCGGACTTGTTGGCCGATGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTCCTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATATCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGCGTGcc                                                                                                              >  2:347942/1‑251 (MQ=255)                                                                                                                                     aGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTgg                                                                                               >  1:560186/1‑251 (MQ=255)                                                                                                                                                  gACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCGGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGGATTTCGCCATCGCCATGCTGCTGATGGTTATTGGCGGTTTTGCCATGAGTGCCGGTCCGCTGAGTTGGGTACTGTGCTCCg                                                                                  <  2:243668/251‑1 (MQ=255)                                                                                                                                                      gTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGGACTGTGCTCCGAAAt                                                                              >  2:371389/1‑251 (MQ=255)                                                                                                                                                         aTTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCATCTTGCCCGCTGGGGACGTAACCCAACGATAACGCTGGGCTTCCTGCCGATGACTGATGGCATGGGCGCACTCGGTAACATGTTGCATATCGGTATTCACTCTCCGTCGGCGTTGTATTTCGACATCGCCATGCTGCTGATGTTTATTGTCGGTTTTACCATGAATGCCGGTCCGCTGATTTGGGTCCTGTGCTCCGAAATTc                                                                            >  1:292738/1‑250 (MQ=255)                                                                                                                                                              cGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGc                                                                      >  1:229072/1‑251 (MQ=255)                                                                                                                                                               ggCCTGACCACCGTACTTGCCGCCTGTATCGCAATCGGCCTTGTTGACCGCTGGTGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCt                                                                     <  1:135533/251‑1 (MQ=255)                                                                                                                                                                gCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTTTTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGCTTTGGGTACTGTGCTCCGAAATTCAGCCGCTg                                                                    >  2:205874/1‑251 (MQ=255)                                                                                                                                                                    gACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGaa                                                                  <  1:10207/249‑1 (MQ=255)                                                                                                                                                                      ccAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGc                                                              >  2:265990/1‑251 (MQ=255)                                                                                                                                                                        aaCGTACTTGCCACCTTGATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCcg                                                            <  2:430779/251‑1 (MQ=255)                                                                                                                                                                         aCGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCcgc                                                           >  2:471236/1‑251 (MQ=255)                                                                                                                                                                          cGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCcgcg                                                          >  1:429558/1‑251 (MQ=255)                                                                                                                                                                             aCTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGAtt                                                       <  1:339267/251‑1 (MQ=255)                                                                                                                                                                                     ccTTTATCGCAATCGGCCTTGTTGACCGCTGGGGAAGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGAACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTCTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGCTTTGGGTACTGTGCTCCGAAATTCAGCAGCTGAAAGGCCGCGATTTTGGCATc                                               >  1:404022/1‑251 (MQ=255)                                                                                                                                                                                         tATCGCAATCGGCCTTGTTGGCCGCTGGGGACGTAAACCAAAGATAACGCTGGGCTTCCTGGCGAAGGATGATGGCATGGACGAACTCGGTACAATGGTGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCCTCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCAGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCt                                           >  2:565580/1‑251 (MQ=255)                                                                                                                                                                                            cGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCt                                        <  1:422547/251‑1 (MQ=255)                                                                                                                                                                                              cAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTAGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTcc                                      >  1:555849/1‑251 (MQ=255)                                                                                                                                                                                                      cTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCCCTGccac                              >  2:104162/1‑251 (MQ=255)                                                                                                                                                                                                               cGCTGGGGACGTAAACCAACGCAAACGCTGGGCGTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATGCACTCTCCGTCTGCGCAGTATTTCGCCATCGCCATGCTGCTGAGGTTTAGTGGCGGTTTTGCCATGAGTGCCTGTCCGCTGATTTGGGTACTTGGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGAt                     <  1:347942/251‑1 (MQ=255)                                                                                                                                                                                                                  tGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCACCTGGATTGc                  >  1:316980/1‑251 (MQ=255)                                                                                                                                                                                                                   ggggACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGCTTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGcc                 >  2:247166/1‑251 (MQ=255)                                                                                                                                                                                                                   ggggACGTAAACCAACGCTAACACTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGc                  <  1:310464/250‑1 (MQ=255)                                                                                                                                                                                                                      gACGTAAACCAACGCTAACGCTGGGCTTCCTGCAGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAAc              >  1:118513/1‑251 (MQ=255)                                                                                                                                                                                                                        cGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCTGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCTATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACAt            <  1:241066/251‑1 (MQ=255)                                                                                                                                                                                                                                 acgctaacgctGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTgg   <  2:535961/251‑1 (MQ=255)                                                                                                                                                                                                                                  cgctaacgctGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTCTTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGCTTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGgc  >  1:391652/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      AGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGCATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGC  >  NC_000913/3089014‑3089489

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