Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,742 T→A 13.4% intergenic (+31/+68) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7420TA13.4% 128.6 / 10.6 52intergenic (+31/+68)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (24/21);  new base A (4/3);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

CACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243598‑1243888
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                       >  1:60717/1‑68 (MQ=17)
                                                                                            cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                                       <  2:60717/68‑1 (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                     <  3:324684/144‑1 (MQ=255)
                                                                                                                               aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                     >  4:324684/2‑145 (MQ=255)
                                                                                                                                 aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                  <  4:10366/130‑1 (MQ=255)
                                                                                                                                 aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                  >  3:10366/1‑130 (MQ=255)
                                                                                                                                              cGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc   <  7:75362/148‑1 (MQ=255)
                                                                                                                                               gTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg  <  7:129113/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243598‑1243888

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: