Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,740 C→G 13.8% intergenic (+29/+70) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7400CG13.8% 138.2 / 11.2 50intergenic (+29/+70)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (24/19);  new base G (4/3);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

CACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTC  >  NC_002947/1243598‑1243887
                                                                                                                                              |                                                                                                                                                     
cACGAAGAGCACGGCTAAGTGTGCCCTGCCTACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGc                                                                                                                                                  <  3:98507/148‑1 (MQ=255)
  cGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCgg                                                                                                                                               >  8:304054/1‑149 (MQ=255)
          aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGttttttca                                                                                                                                       <  4:256357/149‑3 (MQ=39)
             gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGc                                                                                                                                    >  8:197194/1‑149 (MQ=255)
                             cAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGcc                                                                                                                   >  3:298294/1‑148 (MQ=255)
                                   gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGcc                                                                                                                       <  8:35924/138‑1 (MQ=255)
                                     aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                      <  6:264822/123‑1 (MQ=255)
                                     aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                      >  5:264822/1‑123 (MQ=255)
                                     aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCgg                                                                                                          <  3:17802/149‑1 (MQ=255)
                                             gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                                  >  4:285251/1‑149 (MQ=255)
                                                   ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                                            >  1:158431/1‑149 (MQ=255)
                                                    gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaa                                                                                             >  6:255544/1‑147 (MQ=255)
                                                     cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaagag                                                                                          <  5:255544/149‑4 (MQ=255)
                                                       tGAACGGTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTTTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCt                                                                                        <  7:299941/149‑1 (MQ=255)
                                                        gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                       >  2:50103/1‑149 (MQ=255)
                                                         aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                       <  7:304054/148‑1 (MQ=255)
                                                          aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                       >  8:62287/1‑147 (MQ=255)
                                                           cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                    <  4:298294/149‑1 (MQ=255)
                                                            gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                    >  1:276112/1‑148 (MQ=255)
                                                              ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTTCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                    >  5:246538/1‑146 (MQ=255)
                                                                    tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                  >  6:219384/1‑142 (MQ=255)
                                                                    tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                  <  5:219384/142‑1 (MQ=255)
                                                                      aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                          <  3:285251/148‑1 (MQ=255)
                                                                      aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCg                                                                         <  7:197194/149‑1 (MQ=255)
                                                                             tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACGGGCCTTCGTGGCCGTTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGaaa                                                                   <  1:50103/148‑1 (MQ=255)
                                                                                       cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                                        >  8:129113/1‑149 (MQ=255)
                                                                                        tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGACCGCGAATGCCTTTAGCAGATCAATCGTTGAAATAGCGCCAgccg                                                       >  7:285056/1‑149 (MQ=255)
                                                                                         gAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgc                                                      >  2:297274/1‑149 (MQ=255)
                                                                                            cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                                      <  2:60717/68‑1 (MQ=17)
                                                                                            cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                                      >  1:60717/1‑68 (MQ=17)
                                                                                             gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                                           >  7:253533/1‑124 (MQ=255)
                                                                                             gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                                           <  8:253533/124‑1 (MQ=255)
                                                                                             gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCATAAGCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc                                                     >  1:108291/1‑148 (MQ=255)
                                                                                                  cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                                             >  8:75362/1‑149 (MQ=255)
                                                                                                    cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                           <  2:276112/149‑1 (MQ=255)
                                                                                                               tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                 <  7:62287/148‑1 (MQ=255)
                                                                                                                aaTTGCTGCGTGCATGAAAAAAGCGGCCTTCGTGGCGGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGGAGCTCAATGGTTGAAATCGCGGCCGCCGCCCATTTGTTTCCAGCGGTTGAc                                >  7:53452/1‑148 (MQ=255)
                                                                                                                  ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCGGATCAATCGtt                                                                       >  1:27697/1‑107 (MQ=255)
                                                                                                                  ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                       <  2:27697/107‑1 (MQ=255)
                                                                                                                       gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCg                         >  2:215914/1‑148 (MQ=255)
                                                                                                                        cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGc                        <  2:158431/148‑1 (MQ=255)
                                                                                                                        cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa                       <  5:56570/149‑1 (MQ=255)
                                                                                                                        cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCAATGCCATTAGCAGATAAATAGTTTAAATACCGCCAGCCGCCCATTTCATTACAGCGGTTGACGATGccccc                       >  6:56570/1‑146 (MQ=255)
                                                                                                                           gCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                    >  3:157922/1‑149 (MQ=38)
                                                                                                                             ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                    >  2:67314/4‑147 (MQ=255)
                                                                                                                               gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                           >  4:195099/1‑122 (MQ=255)
                                                                                                                               aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                    >  4:324684/2‑145 (MQ=255)
                                                                                                                               aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                    <  3:324684/144‑1 (MQ=255)
                                                                                                                                 aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                 >  3:10366/1‑130 (MQ=255)
                                                                                                                                 aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                 <  4:10366/130‑1 (MQ=255)
                                                                                                                                              cGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc  <  7:75362/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTC  >  NC_002947/1243598‑1243887

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