Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,3850AG63.3% 27.4 / 22.9 31N70N (AAT→AAC) ybcVDUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base A (4/7);  new base G (8/11);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTG  >  NC_000913/579160‑579633
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
acGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAAttt                                                                                                                                                                                                                                            >  1:311784/1‑240 (MQ=255)
                 ttATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCttn                                                                                                                                                                                                                                                   <  1:37953/216‑2 (MQ=255)
                            cGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTGTCGCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTGCTAGAGGATGCGTATTTCTTGAGAATTTAACATTTAGAACCTTTTTAAGTCCTTTTATTAACACGgt                                                                                                                                                                                                     <  1:123348/251‑1 (MQ=255)
                             gTGGAAATATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTAGCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACAttn                                                                                                                                                                                                                                     <  1:29926/218‑2 (MQ=255)
                             gTGGAAAGACTTGTGAAGTGTTCTGAATGCTCTCAGTACCTAGAAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTGGTAGATATTTGTAATCCATCGGAAACCTCCTGCTTTAGCAAGATACTCCCTGCATTGCTGAAATGTGCTTTATCTTGATTTCAACCGATCATAGGACGTATCTCTCAGATGCGTCTTTCTTGCGAATTTAACATTTACACCCTTTTTCAGTCCTTTTATTAACACGgtg                                                                                                                                                                                                    >  1:275552/1‑251 (MQ=25)
                                gAAAGATTTGGGTAGTTGTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAATCTCCTGCTTTAGCAAGAGTTTCCCTTTATTGCTTAAATGTTATTTCTCTTGATTTCAACCCATCATAGGACGTTTCTATAAGATTCGTATTTCTTGAGAATTTAACATTTTCAACCTTTTTAAGTCCTTTTATTAACACTGTGTTa                                                                                                                                                                                                 <  1:302492/251‑1 (MQ=37)
                                                                  aaaTAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGCTATTTGTAATCCATCGGAAAACTCCTGATTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTgg                                                                                                                                                               >  1:105676/1‑251 (MQ=255)
                                                                                   cAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACTTTGCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTAc                                                                                                                                                                                                                                   <  1:213901/166‑1 (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)
                                                                                                                                    gAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGCAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACCAACTTTTTAAGTCCTTTTATTAACACGgtg                                                                                                                                                                                                    >  1:221046/1‑148 (MQ=255)
                                                                                                                                         cTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTACAAGATGCGTGCTTCTTGAGCATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACCCGATGTGAATATTATCTGTGGCTAGATAGTCAATATACTGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCt                                                                                        >  1:299182/1‑251 (MQ=255)
                                                                                                                                         cTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTATCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATNAACNNNnngt                                                                                                                                                                                                   >  2:199981/1‑144 (MQ=255)
                                                                                                                                           ccTGCTTTAGCAAGATTTTCCCTGTATTGATGCAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTATATAAGATGCGTGTTTCTTGAGAATTTAACATTTACACCCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTATAAATCttt                                                                                      >  1:166731/1‑251 (MQ=255)
                                                                                                                                                   aGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTAGCATAGGACGTTTCTAGAAGATGCGTGTGTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAacac                                                                                                                                                                                                        <  1:44280/129‑1 (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)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
ACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTG  >  NC_000913/579160‑579633

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