Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	121,507	0	A	T	41.5%	55.7 / 14.2	41	L15L (CTT→CTA)	aroP	aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (9/15);  new base T (17/0);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 3.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTG  >  NC_000913/121362‑121657                                                                                                                                                  |                                                                                                                                                         gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGGTTTTTAGGcc                                                                                                                                                     >  2:293777/1‑151 (MQ=255)       tCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCccgcgt                                                                                                                                               >  2:664801/1‑149 (MQ=255)        cAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCtt                                                                                                                                               <  2:186188/150‑1 (MQ=255)                                ggCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCt                                                                                                                    <  2:112113/151‑1 (MQ=255)                                       cccAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTTGCTGCTCGCCG‑GGCTG‑TTGATCTTcc                                                                                                                 >  1:638367/1‑147 (MQ=255)                                         cAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTTGCTGCCCGCCG‑GGCTG‑TTGGCCTTCcat                                                                                                               >  1:50899/1‑147 (MQ=255)                                          aGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTAGCTGCTCCCCC‑TGCTG‑TTGGCCTTcc                                                                                                                 >  2:323935/1‑144 (MQ=255)                                                aTCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTc                                                                                                    <  2:455740/151‑1 (MQ=255)                                                  cccTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCg                                                                                                   <  1:706904/150‑1 (MQ=255)                                                   ccTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTg                                                                                                 <  1:323959/151‑1 (MQ=255)                                                     tGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGGCGCGCCTTTGCTGCTCCCCG‑TGCTG‑TTGTCCTTCCCTCATGAAAACTc                                                                                                    >  2:480891/1‑146 (MQ=255)                                                     tGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCC‑GGCTGTTTGCCCTTCC‑TCATGaaaac                                                                                                      >  2:483538/1‑142 (MQ=255)                                                     tGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGCTTTTAAAGGCCCGCTTTTAGTGGTCTGCG‑TGCTG‑TTGGCCTTTCATTATGAAAACTCGTGCg                                                                                               >  1:291572/1‑151 (MQ=255)                                                            gCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCggtggt                                                                                          <  1:409059/149‑1 (MQ=255)                                                            gCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCggtggt                                                                                          <  1:237203/149‑1 (MQ=255)                                                              ggACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGGCGCGCTTTAGCTGCTCGCCG‑GGCTG‑TTGACCTTCCCTCATGGAACCTCGTGCGGGGggtggt                                                                                      >  2:114119/1‑149 (MQ=255)                                                                      ttACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCCTCATGAAACCTCGTGCGGGGGTTGttttttt                                                                                 >  2:659690/1‑148 (MQ=255)                                                                               cgcTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCG‑TGCTG‑TTGACCTTCCCTCATGAAACCCCGTGCGGGGGTTGTTTTTTTTTTCTACGCAg                                                                      >  2:367482/1‑150 (MQ=255)                                                                               cgcTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGt                                                                     <  2:790431/151‑1 (MQ=255)                                                                                  tACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTTATCTACGCAGTg                                                                    >  2:652982/1‑149 (MQ=255)                                                                                       aGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAAACTCGTGCGGTGGGTGTTTTTTTGATCTACGCAGTGATGCgt                                                              >  1:719598/1‑150 (MQ=255)                                                                                             aaCCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTTAGCCGCGCCTTCGGTGCGCCCCG‑TGCTG‑TTGTCCCTCCCTCCTTAAACCTCCTGCGGGGGTTGGTTTTTTTATTTACGCAGTGAGGGgggtgt                                                          >  2:72409/1‑148 (MQ=255)                                                                                                   gTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCG‑GGCTG‑TTGACCTTCCCTCATGAAAACTCGTGCGGTGGTTGTTTTTTTTATTTACGCAGGGATGCGTGGGTAAGttt                                                    >  2:630821/1‑148 (MQ=255)                                                                                                   gTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCa                                                 <  2:324883/151‑1 (MQ=255)                                                                                                      ccTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAAtt                                              <  1:114138/151‑1 (MQ=255)                                                                                                        tATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATTAAAACTCGTGCGGGGGGTGTTTTTTTTATCTACGCAGTGATGCGTGTGTAAGTTTTCAAtt                                              >  1:26833/1‑149 (MQ=255)                                                                                                            gcgcCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCCCGCCTTAGCTGCTCCCCG‑TGCTG‑TTGGCCTTTCATCCTGAAACCTCGGGCGGGGGGTGTTTTTTTTTTCTTCGCAGTGGTGCGGGTGTTAGTTTGTCATTTCGttt                                        >  2:434524/1‑151 (MQ=255)                                                                                                             cgcCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTttg                                       <  1:378498/151‑1 (MQ=255)                                                                                                                cACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCCTTAGCTGCTCGCCG‑TGCTG‑TTTGCCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTTATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTTCGTttgtt                                     >  2:106870/1‑150 (MQ=255)                                                                                                                    cAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAAACTCGTGCGGTGGGTGTTTTTTTTATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTTTAtt                                >  2:369505/1‑151 (MQ=255)                                                                                                                     aGCGCGATAAGCTGAATATGGCGGTTTTTAAGGGCGCGCGTTAGCTGCTCGCCG‑TGGTG‑TTGACCCTCCATCATGAAAACTCGGGCGGTGGGTGTTTTTTTGGTCTAAGCAGGGATGCGTGTGGAAGGTTGCAATTCCGGTTGTTGTAtt                                >  2:175869/1‑150 (MQ=255)                                                                                                                      gcgcGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTAGCCGCTCGCCG‑TGCTG‑GTGTCCTTTCATCATTAAAACTCGGGCGGGGGGTGTTTTTTTTATTTTCGCAGTGATGCGTGTGGAAGTTTTTAATTCCGTTTGTTGtttt                                >  2:402778/1‑149 (MQ=255)                                                                                                                            tAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTATCTGCTCCCCGTTGCGT‑TTTCCCTTCC‑TCATTAAACCCCGGGCGGGGGGTGTTTTTTTTTTTTAAGCAGGGGTGCGGGTGTTAGTTTTCAAATTCGTTTTTTGTATTTATTTGtt                        >  2:755632/1‑151 (MQ=255)                                                                                                                                   aaTATGGCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGGCCTTCCCTCATGAAACCTCGTGCGGTGGGTGTTTTTTTTTTCTTCGCCGTGGTGCGTGGGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTTCAt                   >  2:532830/1‑149 (MQ=255)                                                                                                                                    atatGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCa                 >  1:351768/1‑150 (MQ=255)                                                                                                                                      atGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCaaa               <  2:398062/150‑1 (MQ=255)                                                                                                                                      atGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCaaa               <  1:503759/150‑1 (MQ=255)                                                                                                                                         gCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAaaga             >  1:221537/1‑149 (MQ=255)                                                                                                                                          cGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGt          <  1:35367/151‑1 (MQ=255)                                                                                                                                                 tAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAAtt   >  2:212755/1‑151 (MQ=255)                                                                                                                                                 tAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAAt    <  1:653031/150‑1 (MQ=255)                                                                                                                                                  aaGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTg  <  1:122178/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                         GCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCG‑TGCTG‑TTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATGCGTGTGTAAGTTTGCAATTCCGTTTGTTGTATTAATTTGTTTACATCAAAGAAGTTTGAATTG  >  NC_000913/121362‑121657

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