Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,817,106 T→G 8.2% intergenic (+26/+97) PP_t44 → / ← PP_5526 tRNA‑Pro/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,817,1060TG 8.2% 131.1 / 3.5 49intergenic (+26/+97)PP_t44/PP_5526tRNA‑Pro/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base G (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

GTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGCTGTTCCTCGCCTA  >  NC_002947/2816966‑2817252
                                                                                                                                            |                                                                                                                                                  
gTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGa                                                                                                                                            <  7:240737/149‑1 (MQ=255)
  tctcGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGa                                                                                                                                            <  7:139522/147‑1 (MQ=255)
  tctcGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGa                                                                                                                                            >  8:139522/1‑147 (MQ=255)
      ggCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGa                                                                                                                                       <  1:41713/148‑1 (MQ=255)
      ggCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGAt                                                                                                                                      <  1:160594/149‑1 (MQ=255)
        cAAGATTGAGAACATGCTTCAAAGCACAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGAttt                                                                                                                                    >  1:182644/1‑149 (MQ=255)
          aGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGAttttt                                                                                                                                  <  8:232235/149‑1 (MQ=255)
               gagaACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTAt                                                                                                                              <  2:208893/148‑1 (MQ=255)
                         ttCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGtt                                                                                                                   >  1:14327/1‑149 (MQ=255)
                          tCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGttt                                                                                                                  >  3:199746/1‑149 (MQ=255)
                          tCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGttt                                                                                                                  >  1:264911/1‑149 (MQ=255)
                            aaaGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGtttt                                                                                                                 <  2:264911/148‑1 (MQ=255)
                                cccAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGGCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGtc                                                                                                            <  8:170909/149‑1 (MQ=255)
                                     tCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTctcgct                                                                                                       >  2:211242/1‑149 (MQ=255)
                                        gggCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTctcgctcg                                                                                                     >  6:173674/1‑148 (MQ=255)
                                             tAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAtttt                                                                                               <  1:211242/149‑1 (MQ=255)
                                             tAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAtttt                                                                                               <  5:173674/149‑1 (MQ=255)
                                                     cccGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccact                                                                                       >  5:119736/1‑149 (MQ=255)
                                                        ggTAGCGCACTTGCATGTGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccactttt                                                                                    <  7:231514/149‑1 (MQ=255)
                                                         gTAGCGCACTTGCATGGGGTGCTAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccacttttt                                                                                   <  4:5183/149‑1 (MQ=255)
                                                              gcACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccactttttgtgc                                                                               >  2:27558/1‑148 (MQ=255)
                                                                  ttGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGATCGAATTTTGTGCCACTTTTTGTGACACTg                                                                          >  7:189627/1‑149 (MQ=255)
                                                                              cAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGcc                                                                                                                           >  4:141792/1‑88 (MQ=255)
                                                                              cAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGcc                                                                                                                           <  3:141792/88‑1 (MQ=255)
                                                                                              cGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCt                                              >  1:133065/1‑149 (MQ=255)
                                                                                               gAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGAGCCATCTg                                             >  3:198694/1‑149 (MQ=255)
                                                                                                 aTCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGa                                           <  2:182644/149‑1 (MQ=255)
                                                                                                  tCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGAGCCATCTGGAg                                          <  4:198694/149‑1 (MQ=255)
                                                                                                     cTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAtttt                                                                                               <  1:46250/93‑1 (MQ=255)
                                                                                                     cTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAtttt                                                                                               >  2:46250/1‑93 (MQ=255)
                                                                                                            ccGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTc                                 <  6:272292/148‑1 (MQ=255)
                                                                                                             cGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCgg                               <  2:14327/149‑1 (MQ=255)
                                                                                                                ccAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCgg                               >  3:148064/1‑146 (MQ=255)
                                                                                                                ccAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGgag                            <  1:27558/149‑1 (MQ=255)
                                                                                                                        cataaaaaaaTCCAGTCACCGTTAAGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGga                             >  5:222176/7‑140 (MQ=38)
                                                                                                                        cataaaaaaaTCCAGTCACCGTTAAGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGga                             <  6:222176/134‑1 (MQ=38)
                                                                                                                         ccTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCgg                               >  2:253782/1‑137 (MQ=255)
                                                                                                                         ccTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCgg                               <  1:253782/137‑1 (MQ=255)
                                                                                                                           tAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGgttg                                     <  4:7463/129‑1 (MQ=255)
                                                                                                                           tAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGgttg                                     >  3:7463/1‑129 (MQ=255)
                                                                                                                           tAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGt                 >  6:156448/1‑149 (MQ=255)
                                                                                                                            aaaaaaTCCAGTCACCGTTAAGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGga                             <  7:239669/134‑1 (MQ=38)
                                                                                                                            aaaaaaTCCAGTCACCGTTAAGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGga                             >  8:239669/3‑136 (MQ=38)
                                                                                                                             gAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTg                <  8:189627/148‑1 (MQ=255)
                                                                                                                             gAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGc               >  4:241033/1‑149 (MQ=255)
                                                                                                                               aaaTCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGCTg             <  4:199746/149‑1 (MQ=255)
                                                                                                                                 aTCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGc               <  7:30603/145‑1 (MQ=255)
                                                                                                                                 aTCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGc               >  8:30603/1‑145 (MQ=255)
                                                                                                                                          cTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGCTGTTCCTCGCCTa  <  2:133065/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTTCGTGGCTCCAGATGGCTCATCTGGAGGGTTGTTTCGGGAGATACAGCTCGTGCTGTTCCTCGCCTA  >  NC_002947/2816966‑2817252

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