Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG53.6% 26.4 / 39.5 28G25G (GGT→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (6/7);  new base G (15/0);  total (21/7)
Fisher's exact test for biased strand distribution p-value = 1.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-04

TTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCT  >  NC_000913/523192‑523481
                                                                                                                                               |                                                                                                                                                   
ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTaa                                                                                                                                                 >  1:521011/1‑148 (MQ=255)
  gtgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGaa                                                                                                                                              >  2:315601/1‑149 (MQ=255)
                     tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCCCCggcgtgg                                                                                                                         >  1:520913/1‑151 (MQ=255)
                         ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGGGCCCACCcgcgtggggg                                                                                                                      >  2:239751/1‑147 (MQ=255)
                           gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCggcgtggcg                                                                                                                       >  1:697692/1‑147 (MQ=255)
                           gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCCGGGGTAAGAAATGGCGCGCCCACCcgcgtgggg                                                                                                                       >  2:70258/1‑145 (MQ=255)
                              aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAAATTGCGCGCCCACCcgcgggggggg                                                                                                                     >  2:596544/1‑142 (MQ=255)
                              aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGGGCCCACCggcgtgggggggcccg                                                                                                                >  2:5579/1‑142 (MQ=255)
                                   aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCgg                                                                                                                >  2:273015/1‑146 (MQ=255)
                                   aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCATTGTCCAGGGGTCGGCAGGGGTAAAAAATTGCGCGGCCACCGGCGGGGCGGGCCCgggg                                                                                                              >  2:669223/1‑146 (MQ=255)
                                                 tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGT‑Cccggggggg                                                                                                >  1:507921/1‑146 (MQ=255)
                                                     tAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGGCCGGGGTGT‑CCGGGCGGGGTGACtt                                                                                          >  1:14818/1‑149 (MQ=255)
                                                              aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGGGGTCGGGGTGT‑CCCGGCGGGGTGGCTTCGGGGAAacc                                                                                >  1:513686/1‑150 (MQ=255)
                                                                 ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCgg                                                                              <  2:547013/149‑1 (MQ=255)
                                                                              aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGT‑CCGGGGGGGGTGACTTTGGGCAACCCGGGAAAACCCCTGCTgg                                                               >  1:260392/1‑151 (MQ=255)
                                                                                   aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGGAGGGGTAAGAAATTGCGCGCCCCCCCGCGTGGGGGGGCCGGGGGGTGCGGGGGGGGGGGCCTTCGGGCAACCCGGGAAAAccccgccggggg                                                             >  2:7953/1‑139 (MQ=255)
                                                                                      ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                       <  1:266750/151‑1 (MQ=255)
                                                                                          cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgc                                                   <  2:753641/151‑1 (MQ=255)
                                                                                           gTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCCGGGCCGGTGTTT‑CCGGGGGGGGTGACTTTGGGGAAACCCGGAAAACCCCTGCTGGGGGGGGAggag                                                    >  1:79752/1‑147 (MQ=255)
                                                                                            tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCCGCGGGGCGTGCTCGGGGTGT‑GCGGGGGGGGTGAACTCGGGGAACCCCGGAAAACCGCTGCTGGGGGGGGAGGtg                                                    >  1:655245/1‑148 (MQ=255)
                                                                                              aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGGAAGAATTTGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGT‑CCCGGGGGGGTGAATTCGGGCAACCCGGTAAAACCCCTGCCGGGGGCGAAGGGGCTGcc                                               >  1:35106/1‑151 (MQ=255)
                                                                                                      atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgag                                        <  2:214874/150‑1 (MQ=255)
                                                                                                              gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                               >  1:219975/1‑151 (MQ=255)
                                                                                                                    cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                         <  1:65302/151‑1 (MQ=255)
                                                                                                                     ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgc                        <  1:79334/151‑1 (MQ=255)
                                                                                                                         ccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                    >  1:97426/1‑151 (MQ=255)
                                                                                                                                        cGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCAt     >  2:736945/1‑151 (MQ=255)
                                                                                                                                           cAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtct  <  2:495647/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCT  >  NC_000913/523192‑523481

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: ‑