Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,089,159	(G)5→6	coding (876/1395 nt)	galP →	galactose:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,089,154	1	.	G	95.7%	160.4 / ‑0.3	47	coding (871/1395 nt)	galP	galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (25/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑T‑GGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACT  >  NC_000913/3089017‑3089284                                                                                                                                            |                                                                                                                                   aCAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑T‑gg                                                                                                                                  >  1:247240/1‑140 (MQ=255) aCAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑T‑gg                                                                                                                                  >  1:388202/1‑140 (MQ=255)   aGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGGGCAAATG‑T‑gggg                                                                                                                                >  2:320367/1‑140 (MQ=255)               cgcgcgGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTg                                                                                                                     <  2:247240/140‑1 (MQ=255)                                               aaTGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATc                                                                                     >  2:318598/1‑140 (MQ=255)                                                 tGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGc                                                                                   >  1:228020/1‑140 (MQ=255)                                                 tGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGc                                                                                   >  1:234622/1‑140 (MQ=255)                                                 tGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGc                                                                                   >  2:444692/1‑140 (MQ=255)                                                  gcagcaATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCa                                                                                  <  1:318598/140‑1 (MQ=255)                                                    agcaATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAAt                                                                                <  2:85816/140‑1 (MQ=255)                                                    agcaATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAAt                                                                                <  2:234622/140‑1 (MQ=255)                                                       aaTTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCgg                                                                             >  1:147011/1‑140 (MQ=255)                                                         ttCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGcc                                                                           <  2:388202/140‑1 (MQ=255)                                                               gggATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAAt                                                                                <  1:74309/129‑1 (MQ=255)                                                               gggATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAAt                                                                                >  2:74309/1‑129 (MQ=255)                                                                               tGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGGGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTaa                                                     >  1:310683/1‑140 (MQ=255)                                                                                    tACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaa                                                >  1:418535/1‑140 (MQ=255)                                                                                    tACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaa                                                >  1:144662/1‑140 (MQ=255)                                                                                       gcgcCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaacgc                                             <  2:36564/140‑1 (MQ=255)                                                                                        cgcCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaacgct                                            >  1:357859/1‑140 (MQ=255)                                                                                        cgcCGAAAATCTTCGAACTGGCGGGTTATA‑CAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaacgcta                                           >  2:138515/1‑140 (MQ=255)                                                                                           cGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGTTTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaacgctaac                                         >  1:499842/1‑140 (MQ=255)                                                                                           cGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCaacgctaac                                         >  1:308216/1‑140 (MQ=255)                                                                                                aTCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTgg                                    <  2:357859/140‑1 (MQ=255)                                                                                                 tCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTggg                                   <  2:329978/140‑1 (MQ=255)                                                                                                 tCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTggg                                   <  2:144662/140‑1 (MQ=255)                                                                                                   ttCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCt                                 >  2:229950/1‑140 (MQ=255)                                                                                                      gAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTcc                              <  2:147011/140‑1 (MQ=255)                                                                                                           ggCGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTg                         <  2:228020/140‑1 (MQ=255)                                                                                                             cGGGTTATACCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGAt                       <  1:436120/140‑1 (MQ=255)                                                                                                                    taCCAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGgctgct                >  1:161474/1‑140 (MQ=255)                                                                                                                      ccAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTgg              >  2:101968/1‑140 (MQ=255)                                                                                                                      ccAACACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTgg              >  1:433726/1‑140 (MQ=255)                                                                                                                        aaCACTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCa            <  1:138515/140‑1 (MQ=255)                                                                                                                          cacTACCGAGCAAATGCG‑GGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATg          >  1:47485/1‑140 (MQ=255)                                                                                                                          cacTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGgctgc                 >  1:427656/1‑133 (MQ=255)                                                                                                                          cacTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGgctgc                 <  2:427656/133‑1 (MQ=255)                                                                                                                          cacTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATg          >  2:69275/1‑140 (MQ=255)                                                                                                                          cacTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATg          >  2:203217/1‑140 (MQ=255)                                                                                                                            cTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATgg         >  2:369760/1‑139 (MQ=255)                                                                                                                            cTACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATgg         <  1:369760/139‑1 (MQ=255)                                                                                                                             tACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGc       <  1:444692/140‑1 (MQ=255)                                                                                                                             tACCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGc       <  2:418535/140‑1 (MQ=255)                                                                                                                              aCCGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCg      >  1:237902/1‑140 (MQ=255)                                                                                                                               ccGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGt     >  2:52717/1‑140 (MQ=255)                                                                                                                                cGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTa    <  1:52717/140‑1 (MQ=255)                                                                                                                                cGAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTa    >  2:62049/1‑140 (MQ=255)                                                                                                                                 gAGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATgatccact   >  2:206304/1‑133 (MQ=255)                                                                                                                                  aGCAAATGTT‑GGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACt  <  1:69275/140‑1 (MQ=255)                                                                                                                                  aGCAAATG‑TGGGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACt  <  1:101968/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                   ACAGCAACTTCCGCCGCGCGGTGTTCCTTGGCGTACTGTTGCAGGTAATGCAGCAATTCACCGGGATGAACGTCATCATGTATTACGCGCCGAAAATCTTCGAACTGGCGGGTTATACCAACACTACCGAGCAAATG‑T‑GGGGGACCGTGATTGTCGGCCTGACCAACGTACTTGCCACCTTTATCGCAATCGGCCTTGTTGACCGCTGGGGACGTAAACCAACGCTAACGCTGGGCTTCCTGGTGATGGCTGCTGGCATGGGCGTACT  >  NC_000913/3089017‑3089284

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