Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,4570AG65.5% 28.9 / 19.1 29I46I (ATT→ATC) ybcVDUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base G (9/10);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATT  >  NC_000913/579226‑579703
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aaaTAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGCTATTTGTAATCCATCGGAAAACTCCTGATTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTgg                                                                                                                                                                                                                                     >  1:105676/1‑251 (MQ=255)
                           ntAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTTCTGTCATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACAt                                                                                                                                                                                                          >  1:197852/2‑251 (MQ=255)
                                    tttATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGTCGTTGTGAc                                                                                                                                                                                                    <  1:248978/248‑1 (MQ=255)
                                    tttATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATGAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGAc                                                                                                                                                                                                    <  1:248970/248‑1 (MQ=255)
                                               ggATATTTGTAATCCATCGGAAAACTNCTGCTTNNGCNAGATTTTCCCTGTATTGCTNAANNNTNNNTNCTCTTNATNNCAACCTATCNTAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTNTTAAGTCCTNNTATTAACACGGTGTTATCGTTNTCTNNCNNANTGTNaa                                                                                                                                                                                                                                                 >  2:73894/1‑192 (MQ=255)
                                               ggATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCGAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGt                                                                                                                                                                                            <  1:73894/245‑1 (MQ=255)
                                                          acccATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTGTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTtgtg                                                                                                                                                                                                      <  1:207084/222‑1 (MQ=255)
                                                                 ngAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGACTATTATCTGTGGCTAGATAGTAAATATAAtgtg                                                                                                                                                                                                                >  1:32291/2‑207 (MQ=255)
                                                                       cTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTACAAGATGCGTGCTTCTTGAGCATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACCCGATGTGAATATTATCTGTGGCTAGATAGTCAATATACTGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCt                                                                                                                                                              >  1:299182/1‑251 (MQ=255)
                                                                         ccTGCTTTAGCAAGATTTTCCCTGTATTGATGCAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTATATAAGATGCGTGTTTCTTGAGAATTTAACATTTACACCCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTATAAATCttt                                                                                                                                                            >  1:166731/1‑251 (MQ=255)
                                                                                             ccTGTATTGCTGAAGTGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGAc                                                                                                                                                                                                    >  1:135289/1‑191 (MQ=255)
                                                                                                gTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTTCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATAtttn                                                                                                                                       <  1:33376/249‑2 (MQ=255)
                                                                                                gTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTNc                                                                                                                                       <  1:255525/249‑1 (MQ=255)
                                                                                                                        tGATTTTAACCGATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTTACATTTACAACCTTTGTAAGGTCTTTTATTAACACGGTGATATCGATTTCTAACACGATGTGAATATTATCTGTGGCGAGATAGTAAATATAAtgtg                                                                                                                                                                                                                <  1:174405/152‑1 (MQ=255)
                                                                                                                              cAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTTCAACCTTTTTAAGGCCTTTTATTAACACGATGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTGTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAc                                                                                                                        <  1:321663/234‑1 (MQ=255)
                                                                                                                                 ccTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTANCATNTACANCCNTTTGAAGTTCTTTTANTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGNNGNNNGNTAGTAAat                                                                                                                                                                                                                        <  2:222160/135‑1 (MQ=255)
                                                                                                                                 ccTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTACGTTCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGCTAGTAAat                                                                                                                                                                                                                        >  1:222160/1‑135 (MQ=255)
                                                                                                                                   tATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTCCAACCTCTTTAAGTTGTTTTATTAACACGGCGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATACTGTGAGACGTTGTGACGTTTTAGATCAGAATAAAACAATTCACAGTCTAAATCtttt                                                                                                                                                           >  1:154095/1‑194 (MQ=255)
                                                                                                                                      cATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCttt                                                                                                                                    <  1:111016/214‑1 (MQ=255)
                                                                                                                                                  cTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTTCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATc                                                                                                                                                               <  1:47500/175‑1 (MQ=255)
                                                                                                                                                             gTATTTCTTGAGAATTTAACATTTACAACCTTTGTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCGAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGGTTAAATCTTTTCGCACTTGATCGAATATTTCTGTAAAAAGGGCAGCCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTg                                                                        <  1:286083/251‑1 (MQ=255)
                                                                                                                                                                                        aCNNNNNTAANTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGc                                                                                                                                                        <  2:197852/144‑1 (MQ=255)
                                                                                                                                                                                                  gTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTACAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACctc                                   >  1:17826/1‑251 (MQ=255)
                                                                                                                                                                                                          attaACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAATCCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTgtgtn                           <  1:199981/251‑2 (MQ=255)
                                                                                                                                                                                                               cacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGCTAGTACATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGAACGAATATTTCTTTAAAAATCTCAACCTGAGCCATTGGTAAAACCTTCCATGCGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTttgtt                      >  1:226759/1‑251 (MQ=255)
                                                                                                                                                                                                               cacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGCTAGTAAATCTACTGTGAGAAGTTGTGACGTTTTAGTTCAGAATAAAACACTTCACAGTCTAAATCTTTTCGCCCTTGATCGAATATTTCTTTCAAAATGGCAACCTGa                                                                                                                    >  1:152462/1‑157 (MQ=255)
                                                                                                                                                                                                                          cNTTNTCTAANACGANGTGAATATTATCTGTNGCTAGATAGTAAATATAATNNGNGACGTNGTGAc                                                                                                                                                                                                    <  2:135289/66‑1 (MQ=255)
                                                                                                                                                                                                                               tCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTc                                                                                                  >  1:144267/1‑159 (MQ=255)
                                                                                                                                                                                                                                   acacGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATAtt  >  1:27400/1‑251 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
AAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATT  >  NC_000913/579226‑579703

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