Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913523,3350TG41.9% 68.0 / 12.0 44G25G (GGT→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (7/18);  new base G (18/0);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 8.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523188‑523483
                                                                                                                                                   |                                                                                                                                                      
ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGgtg                                                                                                                                                       <  2:509995/149‑1 (MQ=255)
     tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGaa                                                                                                                                                 >  1:708829/1‑150 (MQ=255)
            cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcgc                                                                                                                                          <  2:736013/150‑1 (MQ=255)
                      ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggc                                                                                                                                <  2:545807/150‑1 (MQ=255)
                           tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCCGGGGTCGGCCGGGGGAAGAATTGGCGCGCCCACCcgcg                                                                                                                               >  1:518060/1‑146 (MQ=255)
                           tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtgg                                                                                                                            <  1:394742/149‑1 (MQ=255)
                            gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGTTCGGCAGGGGTAAGAAATTGCGGCCCCACCcgcgtgggg                                                                                                                          >  2:94206/1‑148 (MQ=255)
                            gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGTAAGAAATTGCGCGCCCACCcgcgggggg                                                                                                                          >  2:422493/1‑148 (MQ=255)
                                  aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTc                                                                                                                     <  1:423537/149‑1 (MQ=255)
                                           atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGGAGGGGTAAGAAATTGCGCGCCCCCCCGCGTGGGGGGGTCGgggtggtccg                                                                                                          >  2:679754/1‑146 (MQ=255)
                                             atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcggg                                                                                                         <  1:370788/150‑1 (MQ=255)
                                                aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcggg                                                                                                     <  1:667101/151‑1 (MQ=255)
                                                     tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGGAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTcggggggg                                                                                                     >  1:641005/1‑146 (MQ=255)
                                                             tgggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCCCCCGCGTGGGGGGCTCGGGGTGTCCGGGCGGGGT‑GACtttcgg                                                                                         >  2:311291/2‑145 (MQ=255)
                                                                 cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCCGGGTGTCCCGGCGGGGG‑GACTTCGGGCAAccc                                                                                   >  1:327074/1‑151 (MQ=255)
                                                                 cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCTCGGTGGGGCCGGGGGGGGT‑GGCTTTGGGCAAacc                                                                                   >  2:473111/1‑151 (MQ=255)
                                                                   cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGCAGGGGTAAGAAATGGGGCGCCCCCCGGCGTGGGGTGCTCGGTGTGTCCGGGGGGGGT‑GACGTCGGGCAAccc                                                                                   >  2:133169/1‑149 (MQ=255)
                                                                     ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCAGGGGTCGGGAGGGGGTAGAATTTGCGCGCCCACCGGCGGGGCGGGGTCGGGGGGTCCCGGCGGGGG‑GACCTCCGGCAAACCgg                                                                                 >  2:465600/1‑149 (MQ=255)
                                                                      tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTTGGCCGGGGTAAGAAATGGGGCGCCCACCCGCGTGGGGGGCTCGGGGTGGCCGGGGGGGGT‑GGATTCGGGCAAacc                                                                                   >  2:635907/1‑146 (MQ=255)
                                                                      tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTTAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCTGGGTGTCCCGGCGGGGT‑GGATTCGGGCAACCCggga                                                                               >  2:197234/1‑148 (MQ=255)
                                                                         atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCCGGGGTAAGAATTTGGGCGCCCACCCGCGTGGCGGGCTCGGGGGGTCCCGGGGGGGT‑GGCTTCGGGGCACCCCGGaaa                                                                             >  2:478967/1‑149 (MQ=255)
                                                                         atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAAt                                                                            <  1:413637/150‑1 (MQ=255)
                                                                         atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAAt                                                                            <  2:42167/150‑1 (MQ=255)
                                                                                gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCg‑ctgct                                                                   <  1:559488/151‑1 (MQ=255)
                                                                                 aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTTCGGGCGGGGTGGCCTTGGGCCACCCCGGAAAATccccgc                                                                      >  2:219937/1‑145 (MQ=255)
                                                                                 aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGCGGCGTGCTCGGGGTGGCCGGGGGGGGT‑GAACTCGGGGAAACCGGGAAAACCg‑c                                                                       >  2:509160/1‑146 (MQ=255)
                                                                                    aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTgggg                                                               <  2:224673/151‑1 (MQ=255)
                                                                                      cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAAGGGTCGGCGGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGGTCGGTGGGGCCGGGGGGGGT‑GGATTCTGGCAACCCGGGAAATTCG‑CTCCTggggg                                                              >  1:705472/1‑150 (MQ=255)
                                                                                         cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGaa                                                          <  1:569539/151‑1 (MQ=255)
                                                                                            cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGGGTGGCTGGGGTGTCCCGGCGGGGG‑GACTTCCGGGAACCCGGGGAATCCC‑CCGCTGGGGGCgccggg                                                       >  1:613763/1‑146 (MQ=255)
                                                                                              cgTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCCGCGTGGCGGGCCCGGTGTGTCCGGGGGGGGT‑GGCCTCGGGCAACCCCGGAAAACCG‑CCGCTGGGGGCGGAg                                                         >  2:61051/1‑147 (MQ=255)
                                                                                              cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGTGTGTTCGGGGGGGGT‑GGCTTCGGGCAACCCCGGAAATCCC‑CTGCCGGGGGCGGAgggg                                                      >  2:303561/1‑148 (MQ=255)
                                                                                                    ggggaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGTGCCGACCGGCGTGGCGTGCTCGGTGTGTCCGGACGGGAT‑GCCTTCGGGCAGCCCGGTAAATCCG‑CTGCTGGGGAAGAAGGTTCTACCCg                                               <  2:192076/148‑1 (MQ=255)
                                                                                                    ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTCGGGCGGGGT‑GACTTTGGGCAACCCCGGAAAACCG‑CTGCTGGGGGGGAAGGGGCTGccc                                                >  1:589909/1‑150 (MQ=255)
                                                                                                         tatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGCGTAAGAATTGGCGCCCCCACCGGCGTGGCGTGCTCGGTGTGGCCGGGCGGGGT‑GACTTCGGGCAACCCCGCAAAACCG‑CTGGTGGGGGGGAAGGGGCCGCCCGGCg                                            >  2:384122/1‑149 (MQ=255)
                                                                                                          atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCgag                                          <  2:312354/150‑1 (MQ=255)
                                                                                                                cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGa                                     <  1:568731/149‑1 (MQ=255)
                                                                                                                cAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGGAAGAATTGGGGCGGCCACCCGCGTGGCGTGCGCGGTGTGTCCGGGGGGGGG‑TACTTTGGGGAACCCCGGAAATCCGCCTCCTGGGGG‑GGAGGTGCTGTCCCGCGGGGCgg                                      >  2:285285/1‑148 (MQ=255)
                                                                                                                  gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                                 <  1:50633/151‑1 (MQ=255)
                                                                                                                      ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                             >  2:557184/1‑151 (MQ=255)
                                                                                                                        cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                           >  1:308134/1‑151 (MQ=255)
                                                                                                                           tCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                        >  1:156446/1‑151 (MQ=255)
                                                                                                                            cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                       <  1:362501/151‑1 (MQ=255)
                                                                                                                                 tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                   <  2:768900/150‑1 (MQ=255)
                                                                                                                                    ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTg               <  1:59193/151‑1 (MQ=255)
                                                                                                                                                  ctgtAAGAATTCGCGCGCCCACCGGCGTGTCGTGCTCGGGGAGTCCGAGCGGGAT‑GACTTCCAGAAACCCCGCAAATCCG‑CTGCTGGGGGTGCAGGTGATGCCCGGCGAGACGGACCTTGCCCTGCCCGGCCCGCCGCCGATCCTtctct  >  2:656688/2‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT  >  NC_000913/523188‑523483

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