Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,580	0	T	G	45.2%	53.3 / 15.6	43	I235L (ATC→CTC)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (16/7);  new base G (0/19);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 2.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTT  >  NC_000913/907432‑907729                                                                                                                                                       |                                                                                                                                                       aGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGggggatg                                                                                                                                                       <  1:372208/150‑1 (MQ=255)     gTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGggggatggggat                                                                                                                                                  >  1:737117/1‑151 (MQ=255)               gAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAgcgg                                                                                                                                       <  2:2380/151‑1 (MQ=255)                aGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAgcggc                                                                                                                                      <  2:539741/151‑1 (MQ=255)                    cACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCtt                                                                                                                                  <  1:935059/151‑1 (MQ=255)                                  tcatcCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGcc                                                                                                                     >  2:331258/1‑150 (MQ=255)                                   catcCGCCAGCAGATCGTGTTTCATCCCCTG‑AATCAACGCCCGGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCa                                                                                                                    <  2:177872/150‑1 (MQ=255)                                   catcCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCACCGCCCTGGCGGTGGGGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCCGCGGGTGGGGAGGGGGATGGGGATCAC‑CCAGCGGCGCTTACCGCTCACCGCCAt                                                                                                                   <  2:154626/151‑1 (MQ=255)                                      ccgccgGCAGAGCGTTTTTCAGCCCCTG‑AATCAACGCCCGGGGGGTGGGGGGGGGGACGAAATTAACCACATGGAAAACCCA‑CCCCGAAACCCGGCGGGGGGGGGGGGGAAGGGGGATAAA‑CCAGCGGCGCTTACCCCCCACCGCCATaaa                                                                                                                <  2:226606/151‑1 (MQ=255)                                                      tttCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                >  2:503699/1‑151 (MQ=255)                                                      tttCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                >  1:498730/1‑151 (MQ=255)                                                      ttccAGCCCCTG‑AACCACCGCCCTGGGGGTGGGGGGGGGGAGGAAAATAATCACATTAAAAACCCA‑CCCCAAAACCGGGCGGGGGGGGGGGGGGAGGGTGATCAC‑CAGGGGGGGCTTACCGCTACCCGCCATAAAATGTTCAAACTGTTg                                                                                                 <  2:661376/147‑1 (MQ=255)                                                       tttcaccccTGAAATAAACCCCCTGGGGGGGGGGGGGGGTGGGAAAGTAATCACATTGAAAACCCC‑CCCCAAAACCCGGGGGGTGGGGGGGGGGAGGGGGATCAC‑CCAGGGGCGTTTACCGCTCACCGCCATAAAATGTTCAAACt                                                                                                     <  1:603073/141‑1 (MQ=255)                                                       ttCAGCCCCGA‑AAACACCCCCCTGGCGGTGGGGGGCGGTACGAAAGTAACCACTTTGAAAACCCA‑CCCCGAAACCCAGCGGGGGGGGGGGGGGAGGGGGATCAC‑CCAGCGGCTTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                <  1:888003/150‑1 (MQ=255)                                                          aGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                            >  1:781666/1‑151 (MQ=255)                                                            ccctTG‑AATCACCGCCTGGGGGGGGGTGGGCGGTAGGAAAATAATCACTTGGAAAACCCC‑CCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAC‑CGAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                            <  1:861779/149‑1 (MQ=255)                                                            ccccTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                           >  1:109912/1‑150 (MQ=255)                                                             cccTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAAc                                                                                         >  1:578055/1‑151 (MQ=255)                                                              ccTG‑AATCAACGCCCGGGCGGGGGGGGGGGGTCCGAAATTAATCAAATTGAAAACCCA‑CCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCAA‑CCGGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                        <  1:454823/151‑1 (MQ=255)                                                               ccgaaaTCAACGCCCGGGGGGTGGGGGGCGGTACGGAAGTAACCACATTGAAAACCCC‑CACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAA‑CCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                        <  1:291935/147‑1 (MQ=255)                                                                         cGCCCGGGGGGTGGTGGGGGGTGGGAAAGTAATCACATT‑AAAAACCACCCCCAAAACCCAGGGGGGGGGGGGGGGAAGGGGGACCAA‑CGAGGGGCGTTTCCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                               <  2:702929/150‑1 (MQ=255)                                                                              tGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                         >  1:374862/1‑151 (MQ=255)                                                                              ggggggtggtGGGGGGTAGGGAATTAATCACATGGAAAAACCA‑CACCGAAACCCAGCGGGTGGGGAGGGGGAGGGGGATCAC‑GCAGCGGCGTTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGt                                                                           <  2:928614/145‑1 (MQ=255)                                                                                   gtggtgGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                    >  1:499506/1‑151 (MQ=255)                                                                                     ggtggGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCggcggc                                                                   >  2:238760/1‑150 (MQ=255)                                                                                        gggCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                               >  1:547795/1‑151 (MQ=255)                                                                                             gTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGa                                                              >  2:285839/1‑147 (MQ=255)                                                                                              tACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                         >  2:583192/1‑151 (MQ=255)                                                                                                     gTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                  >  1:849823/1‑151 (MQ=255)                                                                                                      tAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                 >  1:172068/1‑151 (MQ=255)                                                                                                           aaaTTGAAAACCCA‑CAAAAAAACCCCGGGGGGGGGGGGGGGGAGGGGGACCAA‑CCAGGGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                             <  1:543901/148‑1 (MQ=255)                                                                                                            caTTGAAAACCCA‑CCCCAAAACCCAGGGGGTGGGGGGGGGGAGGGGGATAAA‑CCAGGGGCGCTCCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTgg                                            <  2:412223/150‑1 (MQ=255)                                                                                                                 agaACCCC‑CCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGAATCAC‑CCGGCGGCGCTTACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                         <  1:62361/148‑1 (MQ=255)                                                                                                                   aacccA‑CACCAAACCCCAGGGGGGGGGGGGGGGGGGGGGGATCAC‑CCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatga                                    <  1:755079/148‑1 (MQ=255)                                                                                                                     cccA‑CACAGAAACCCGGCGGGTGGGGGGGGGGTGGGGGACCAC‑CCAGCGGCCCTCCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATc                                  <  1:371575/150‑1 (MQ=255)                                                                                                                        ctccccGAAACCCAGGGGGGGGGGGGGGGGATGGGGATCCA‑CCATCGCCGTTACCCGCTCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                   <  1:430004/144‑1 (MQ=255)                                                                                                                          ccccaaaacccgggggggggggggggggAGGGGGATCAC‑CCGGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                               <  1:583340/136‑1 (MQ=255)                                                                                                                             caaaaCCCAGGGGGGGGGGGGGGGGAGGGGGATAAC‑CCGGCGCCGTTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaa                            <  2:174605/148‑1 (MQ=255)                                                                                                                             cGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                           >  1:399709/1‑151 (MQ=255)                                                                                                                                   ccgggggggggggggggggAGGGGGCACAC‑CGGGGGGCGTTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                     <  2:379639/146‑1 (MQ=255)                                                                                                                                   ccggCGGGGGGGGAGGGGGAGGGGGATCAA‑CGAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                         <  1:220598/144‑1 (MQ=255)                                                                                                                                      gCGGGGGGGGAGGGGGATGGGGAAAAA‑CCAGCGCCGCTTGCCGCGCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAg                   <  1:829219/150‑1 (MQ=255)                                                                                                                                               gggggggTGGGGGATCAC‑CCGGCGCCGTTCCCCTCCCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc         <  2:745125/149‑1 (MQ=255)                                                                                                                                               gAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc         >  2:570896/1‑151 (MQ=255)                                                                                                                                                 gggggA‑GGGGAACAAGCCGGCGGCTTTTCCCGCCAGCCGCCATAAAATGTACAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcag          <  2:970273/148‑1 (MQ=255)                                                                                                                                                      atggggatCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt  >  2:273656/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTG‑AATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCA‑CACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAA‑CCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTT  >  NC_000913/907432‑907729

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