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	98.4%	215.3 / NA	64	V326A (GTG→GCG)	galP	galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/32);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGTTCCTGACCAT  >  NC_000913/3089014‑3089506                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        agaACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATg                                                                                                                                                                                                                                                     >  1:298399/1‑251 (MQ=255)    aCAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGgct                                                                                                                                                                                                                                                  <  1:370034/251‑1 (MQ=255)                cccgcGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGc                                                                                                                                                                                                                                      <  2:100048/250‑1 (MQ=255)                        gtgtTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCgg                                                                                                                                                                                                                              >  1:229501/1‑251 (MQ=255)                        gtgtTCCGTTGCGGACTGTTGCAGGTAATGCAGCAATCCACCGGGATGAACGTCATCATGTATTACTCTCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCTAGCAAATGTTGGTGACCGTGATTGTCGGCCTGACCAACGGACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGTCGTAAACCAACGCTAACGCTGGGCCTCCTTGCGATGTCCGCTGGCATGGTCGTACTCgg                                                                                                                                                                                                                              <  2:540426/251‑1 (MQ=255)                              cTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGGTTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTTACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAat                                                                                                                                                                                                                        >  1:146000/1‑251 (MQ=255)                                gggcggacggTTGCAGGTAATGCAGCAAGTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAAGGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTGGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAatga                                                                                                                                                                                                                      <  2:66254/242‑1 (MQ=255)                                   cGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCCTGTTGACCGCTGGGGACGTAAACCAAAGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCCTGGGCGTACTCGGGACAATGATGc                                                                                                                                                                                                                   >  1:516858/1‑251 (MQ=255)                                      ccTGTTGCAGGTAAGGCAGCAATTCACCGTGAGGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCTAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ata                                                                                                                                                                                                               <  2:179743/250‑1 (MQ=255)                                         gTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCg                                                                                                                                                                                                            <  1:459185/251‑1 (MQ=255)                                               ggTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCGTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGTTATTc                                                                                                                                                                                                      <  1:398962/251‑1 (MQ=255)                                                   aTGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCCACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCCAGGGCGGACTCGGTACAATGATGC‑ATATCGGTATTCActc                                                                                                                                                                                                  >  2:486675/1‑251 (MQ=255)                                                    tGCAGCAATTCACCGGGCTGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGATATACCAACACTAACGAGCAAAATTGGGGGACCGGGATTGTCCGCCTGACCAACGGACTTACCACCTTTATCGCAATCGGCCTTGTTGAAAGCCGGGGACGTAAAACCACCCTAACGCTGGGCCTCCCTGCGGTGGGTGCTGGGCTGGGCGTACCCGGTGCAAAGATTC‑CAATCGGTAATCActct                                                                                                                                                                                                 >  2:518048/1‑251 (MQ=255)                                                     gcagcaATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACtctc                                                                                                                                                                                                <  1:374297/251‑1 (MQ=255)                                                       agcaATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCg                                                                                                                                                                                              >  1:261457/1‑251 (MQ=255)                                                          aaTTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGAACAATGATGC‑ATATCGGTATTCACTCTCCGTCg                                                                                                                                                                                           >  1:474383/1‑251 (MQ=255)                                                                             catcatGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGGCCGTGATTGTCGGCCTGACCAACGTACTTGCCAACTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAAACAACGCTAACGCTGGGCTTCCTGGCGATGGCGGCTGGCATGGGCGTAATCGGTACAATGATGCCACA‑CGGGATTCACTCTCCGTCGGCGCAGGATTtcgccatcg                                                                                                                                                                        >  2:438415/1‑251 (MQ=255)                                                                               tcatGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTtcgccctcgcc                                                                                                                                                                      >  1:320881/1‑251 (MQ=255)                                                                                catGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTtcgccatcgcca                                                                                                                                                                     <  1:470405/251‑1 (MQ=255)                                                                                catGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTtcgccatcgcc                                                                                                                                                                      >  1:389882/1‑250 (MQ=255)                                                                                             ccGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGtt                                                                                                                                                        <  1:234956/251‑1 (MQ=255)                                                                                                           aCTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTg                                                                                                                                          <  1:390995/251‑1 (MQ=255)                                                                                                            cTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGCTGTTTATTGTCGGtttttc                                                                                                                                         >  2:272178/1‑249 (MQ=255)                                                                                                              ggCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGGTTATTGTCGGTTTTGcc                                                                                                                                        >  1:140080/1‑250 (MQ=255)                                                                                                                 gggTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGccctgg                                                                                                                                    >  2:297491/1‑247 (MQ=255)                                                                                                                           aacccgaCCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTcc                                                                                                                          <  2:298399/245‑1 (MQ=255)                                                                                                                            acccTACCGAGCAAATGTTGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCCCCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCg                                                                                                                         <  1:387998/248‑1 (MQ=255)                                                                                                                               cTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTg                                                                                                                      >  1:21557/1‑251 (MQ=255)                                                                                                                                 aCCGAGAAAATGTGGGGTACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGAGGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGAt                                                                                                                    <  2:453311/251‑1 (MQ=255)                                                                                                                                    gAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTCTTCACTCTACGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGTGTGCCGGGCCGCTGATTTg                                                                                                                 >  2:40088/1‑251 (MQ=255)                                                                                                                                    gAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTg                                                                                                                 <  1:389966/251‑1 (MQ=255)                                                                                                                                      gCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTggg                                                                                                               <  1:133650/251‑1 (MQ=255)                                                                                                                                          aTGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACt                                                                                                           <  1:129238/251‑1 (MQ=255)                                                                                                                                            gtgGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACtg                                                                                                          <  1:274693/250‑1 (MQ=255)                                                                                                                                              gggggACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGc                                                                                                       >  1:294331/1‑251 (MQ=255)                                                                                                                                                      gTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCACCGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCTTGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTTTTTCGCCATCGCCATGCTGCTGATGTTTTTTGTCGTTTTTGCCATGATTGCCGGCCCCCTGATTTGGGTACTGTGCTCCGAAAt                                                                                               >  2:492335/1‑251 (MQ=255)                                                                                                                                                           tGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGc                                                                                          >  1:55172/1‑251 (MQ=255)                                                                                                                                                              cGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCg                                                                                        <  2:294331/250‑1 (MQ=255)                                                                                                                                                                   tGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATACTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGGGCTCCGAAATTCAGCCGCTGaa                                                                                   >  1:3861/1‑250 (MQ=255)                                                                                                                                                                     aCCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATTAGTGCCGGTCCGCTGATTTTGGTACTGTGCTCCGAAATTCAGCCGCTGAAAgg                                                                                >  2:132089/1‑251 (MQ=255)                                                                                                                                                                         aCGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCcgc                                                                            <  2:74900/251‑1 (MQ=255)                                                                                                                                                                                tGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGAtttt                                                                      <  1:519764/250‑1 (MQ=255)                                                                                                                                                                                      cTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATcc                                                               >  1:26282/1‑250 (MQ=255)                                                                                                                                                                                      cTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATc                                                                >  1:16168/1‑250 (MQ=255)                                                                                                                                                                                          agcgcagccGGCCTTGTTGACGTCTTGGGACGTAAACCAACGCTAAAGCTGGGCCTCCTGGTGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATTGGTATTCACTCTCCGGCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTTGGTTTTGCCAGGAGTGCCGGTCCGCTGATGTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCt                                                            <  2:32131/242‑1 (MQ=255)                                                                                                                                                                                          aTCGCAATCGGCCGTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCt                                                            <  2:229501/250‑1 (MQ=255)                                                                                                                                                                                                  cGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACGCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACt                                                    <  2:474383/250‑1 (MQ=255)                                                                                                                                                                                                        tgttgACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCCTGAGTGCCGGTCCGCTGGTTTGGGTCCTGGGCTCCGAAATTCAGCCGCTGACAGGCCGCGATTTTGGCATCACCTGCTCCCCTGcccccc                                             >  2:44531/1‑247 (MQ=255)                                                                                                                                                                                                         gttgACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTTGCATGGGCGTACTCGGTACAATTATGC‑ATAGCGGTATTCACTCGCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGATTGCCTGTCCGCTGATGTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGccacca                                             <  1:554518/250‑1 (MQ=255)                                                                                                                                                                                                           tgACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAAc                                           <  1:298066/250‑1 (MQ=255)                                                                                                                                                                                                            gcccGCTGGGGTCGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGTTGCTGGCATGGGCGTACCCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTg                                         <  2:330656/249‑1 (MQ=255)                                                                                                                                                                                                            gACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGATCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCTCCTGCTCCACTGCCACCAACTg                                         >  2:88591/1‑251 (MQ=255)                                                                                                                                                                                                               cGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGAt                                      >  1:520902/1‑251 (MQ=255)                                                                                                                                                                                                                   ggggACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCCTGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCTCTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCTCCTGCTCCACTGCCACCAACTGGATTGcc                                  <  2:207478/251‑1 (MQ=255)                                                                                                                                                                                                                      gACGTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAAc                               <  1:556375/251‑1 (MQ=255)                                                                                                                                                                                                                         gTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCCCCTGCTCCACTGCCACCAACTGGATTGCCACCATg                            >  2:45620/1‑251 (MQ=255)                                                                                                                                                                                                                         gTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATg                            <  1:438415/251‑1 (MQ=255)                                                                                                                                                                                                                         gTAAACCAACGCTAACGCTGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATg                            >  1:265852/1‑251 (MQ=255)                                                                                                                                                                                                                               cAACGCTAACGCTGGGCTTCCTGGCGATGGCGGCTGGCATGGGCGTACTCGGTACAATGAGGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCTCCAACTGGATTGCCAACATGATCGtt                      <  2:16168/251‑1 (MQ=255)                                                                                                                                                                                                                                         gctGGGCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTTTTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGt            >  1:476163/1‑251 (MQ=255)                                                                                                                                                                                                                                           ggggCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATAGCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGGGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGTTc          <  2:265852/250‑1 (MQ=255)                                                                                                                                                                                                                                            gggCTTCCTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGTTcc         >  1:339349/1‑251 (MQ=255)                                                                                                                                                                                                                                                  ccTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATTC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGTTCCTGACCa   <  2:102877/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cTGGCGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGTTCCTGACCAt  >  1:81401/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        AGAACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATGTGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACTCGGTACAATGATGC‑ATATCGGTATTCACTCTCCGTCGGCGCAGTATTTCGCCATCGCCATGCTGCTGATGTTTATTGTCGGTTTTGCCATGAGTGCCGGTCCGCTGATTTGGGTACTGTGCTCCGAAATTCAGCCGCTGAAAGGCCGCGATTTTGGCATCACCTGCTCCACTGCCACCAACTGGATTGCCAACATGATCGTTGGCGCAACGTTCCTGACCAT  >  NC_000913/3089014‑3089506

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