Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,736,789	0	G	C	45.5%	62.2 / 13.2	44	L168V (CTG→GTG)	ydhP	putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (20/4);  new base C (1/19);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.09e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCAGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGTTTTGGCACCACGCTTGCG  >  NC_000913/1736644‑1736937                                                                                                                                                   |                                                                                                                                                     gCACCAGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGcccc                                                                                                                                                  <  1:471658/150‑1 (MQ=255)  cACCAGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCccgc                                                                                                                                                >  2:10598/1‑151 (MQ=255)   aCCAGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCccgcc                                                                                                                                               >  1:877532/1‑151 (MQ=255)    ccAGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCccgc                                                                                                                                                <  1:84434/149‑1 (MQ=255)      aGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGtt                                                                                                                                            >  1:118330/1‑151 (MQ=255)          cccAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGcc                                                                                                                                         >  1:703492/1‑150 (MQ=255)                 ccGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATg                                                                                                                                  >  2:360557/1‑150 (MQ=255)                                              gggctcATAAATCCCGCCATTTTTTTTTTCTCTTAGGGCCCTTCCCCTGCCCCCCTTTAAGGTAATAAGAAAAAAACACTTCCCATGA‑AATTCCCTCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAAcc                                                                                                    <  1:736846/146‑1 (MQ=255)                                                   cATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTc                                                                                               >  2:304349/1‑151 (MQ=255)                                                        aTACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTcgccgc                                                                                          >  1:370925/1‑151 (MQ=255)                                                           cccccAGCTTTTTTTTCACTCAAGTCCGTCCCCCTCCCCCCCTTTTTGGT‑ATGAAAAAAAACAACCTTCCATTGAAAACCCCTCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCa                                                                                                   <  1:522706/136‑1 (MQ=255)                                                                  ctctTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGc                                                                                  >  2:847172/1‑149 (MQ=255)                                                                   ttttttttCACTCCAGGTCGTCCCCCTGCCCCCCTTTTAGGTAATGAAAAAAACAGACTTACCATTA‑AAATCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGcc                                                                                 <  1:882548/147‑1 (MQ=255)                                                                           cACTTCGGGTCGTCCCCCTCCCCCCCTTTTAGGTAATGAAAAAAACAGACTTCCAATTA‑AAACCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATAtt                                                                        <  1:380793/150‑1 (MQ=255)                                                                               tCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGATGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCACCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGcc                                                                    >  1:26396/1‑150 (MQ=255)                                                                                  ggTCGTCCCCCTGCCCCCCTTTTGGTTAATGAAAAAAACAGATTTCCCATTA‑AAACCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAg                                                                  <  1:129237/149‑1 (MQ=255)                                                                                  ggTCGTCCCCCTCCCCCCCTTTTAGTTAAGAAGAAAAACAAACTTCCCTTAA‑AAACCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGc                                                                     <  1:322686/146‑1 (MQ=255)                                                                                  ggTCGTCCCCCTCCCCCCCTTTTAGGTAATGAAAAAAACAGACTTACCATTA‑AAACCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGc                                                                     <  1:805217/146‑1 (MQ=255)                                                                                    tCGTGCCCCGCCCCCCCCTTGAGTAAATGAGAAAACAACACTCCCAATGA‑AAAACCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGt                                                               <  2:400286/150‑1 (MQ=255)                                                                                      ggcccccTCCACCCCTTTTAGGTAAAGAAAAAAACAGACTCACCATAG‑AACCCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTa                                                             <  2:455121/147‑1 (MQ=255)                                                                                        gCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAAcc                                                          >  2:618790/1‑151 (MQ=255)                                                                                         ccccgcccccccttttGGGAAAGGAAAAAAACAGCCTTCCCATGG‑AAATCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAAcc                                                          <  1:394428/137‑1 (MQ=255)                                                                                          cccggccccccttttAGTAAATGAGAAAAACAGCTTTACCATTA‑AAACCCCCCCCCCCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAAcccc                                                        <  1:823890/139‑1 (MQ=255)                                                                                               aaccaccTTTAGGAAATGAGAAGAACAGACTTACCATAG‑AAATCACTACCAGCCCCGCAGTTGCCAGAAATGACAGCCGCCAGCCGATGGCTTCATCCAACCAGGTCGCTGCCGCCACGCCACCGATATTTGCCAGGGTTAACCCCATaa                                                    >  2:271603/2‑150 (MQ=255)                                                                                                  caccTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACAt                                                >  2:377576/1‑151 (MQ=255)                                                                                                      tttAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTg                                            >  2:212726/1‑151 (MQ=255)                                                                                                       ttAGTTAAGGAAAAAAACAGACTTCCCATTA‑AAACCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTg                                            <  1:708894/150‑1 (MQ=255)                                                                                                       ttAGATAAAGAAAAGAACGGATTTA‑CATTGAAAACCCCTCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAg                                              <  2:64511/148‑1 (MQ=255)                                                                                                        tAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCAGATGGTTTCACCCAACCAGGGCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTg                                            >  2:16274/1‑149 (MQ=255)                                                                                                          ggTAATGAGAAGAAAAATCTTCCAATAG‑AAATCCCCCCCCCCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGt                                             <  1:433692/146‑1 (MQ=255)                                                                                                             aagaaaaagaaCACATTTCCCATGA‑AAATCCCCCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACt                                      <  2:263854/144‑1 (MQ=255)                                                                                                              aTGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTgc                                    >  1:918056/1‑151 (MQ=255)                                                                                                               tGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATACACATAGTGGCAACTgc                                    >  2:152724/1‑150 (MQ=255)                                                                                                                   aaaaaaaGATTTCCCATTA‑AAACCCCCCCCACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCt                                 <  1:109002/143‑1 (MQ=255)                                                                                                                    agaaCAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCACCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGc                              >  1:717515/1‑151 (MQ=255)                                                                                                                            cTTACCATTG‑AAATCACTCCCAGCCCCTCCGTTGCCAGAAATGACATCCGCCAGCCTATGGTTTCACCCAACCAGGTCGCCGCCGGCCCGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTAt                      >  2:797917/1‑151 (MQ=255)                                                                                                                              tACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGt                    >  2:357100/1‑151 (MQ=255)                                                                                                                                acattaaaacccccccccACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGtt                   <  2:600277/136‑1 (MQ=255)                                                                                                                                    tgaaaaTCCCTCCCGGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGttt                        <  1:23461/139‑1 (MQ=255)                                                                                                                                        acccccacccACCCCCGCCGTGGCCAGAAAGGACAGCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGTTTTg                <  2:634131/139‑1 (MQ=255)                                                                                                                                           cccccccACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGTTTTGGcaccac         <  1:343747/146‑1 (MQ=255)                                                                                                                                           cccccccACCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGTTTTGGCACCACg        <  2:609498/147‑1 (MQ=255)                                                                                                                                           cACTACCACCGCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACTCCACAAACATACTGGCAACTGCGCTGGCCTGTTTATGTTTTGGccccac         >  2:155824/1‑150 (MQ=255)                                                                                                                                                 cAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGTTTTGGCACCACGCTTGCg  >  1:127540/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCACCAGCTCCCAGTACCGTCGTCAGCAATGCAGACAGCACCTGCGGACGCATTAATACCGCCAGCTCTTTTTTCACTTCAGGTCGTGCCCCTGCACCACCTTTAGGTAATGAGAAGAACAGACTTACCATTG‑AAATCACTCCCAGCCCCGCCGTTGCCAGAAATGACATCCGCCAGCCGATGGTTTCACCCAACCAGGTCGCCGCCGGCACGCCACCGATATTTGCCAGGGTTAACCCCATAAACATAGTGGCAACTGCGCTGGCCTGTTTATGTTTTGGCACCACGCTTGCG  >  NC_000913/1736644‑1736937

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