Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,387,896 T→C 10.6% *198W (TGA→TGG)  PP_1209 ← cold‑shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,387,8960TC10.6% 121.4 / 5.4 47*198W (TGA→TGG) PP_1209cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (2/3);  total (23/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 = 8.98e-01

AAGAGGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAG  >  NC_002947/1387749‑1388043
                                                                                                                                                   |                                                                                                                                                   
aaGAGGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTc                                                                                                                                                    <  6:38365/149‑1 (MQ=255)
   agGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAAc                                                                                                                                                 <  8:219482/149‑1 (MQ=255)
     gCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCTCAGGTTTGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGACCAAcgg                                                                                                                                               <  2:246581/149‑1 (MQ=255)
        ccGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAAcggcgg                                                                                                                                            <  4:315435/149‑1 (MQ=255)
          gCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGACCAAcggctgac                                                                                                                                          <  6:94579/149‑5 (MQ=255)
              cctcATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGACGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACt                                                                                                                                      >  7:60879/1‑149 (MQ=255)
               ctcATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGGCGGTCAGCCGTTGGTCAACGGCGGTTTACt                                                                                                                                      <  1:30169/148‑1 (MQ=255)
               ctcATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTg                                                                                                                                     >  6:353051/1‑149 (MQ=255)
                tcATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTg                                                                                                                                     >  1:283531/1‑148 (MQ=255)
                tcATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTg                                                                                                                                     >  3:91525/1‑148 (MQ=255)
                 cATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCCCTGGTTAGCCGATGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCCGTTTAATGc                                                                                                                                    >  4:57951/1‑148 (MQ=255)
                      ttGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCAc                                                                                                                              <  6:145444/149‑1 (MQ=255)
                      ttGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCAc                                                                                                                              <  4:91525/149‑1 (MQ=255)
                         aGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTc                                                                                                                           >  8:75407/1‑149 (MQ=255)
                         aGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTc                                                                                                                           >  1:162019/1‑149 (MQ=255)
                              ttCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCgg                                                                                                                      >  8:351799/1‑149 (MQ=255)
                              ttCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCg                                                                                                                       <  8:216197/148‑1 (MQ=255)
                                     ctgctATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAg                                                                                                               <  3:57951/149‑1 (MQ=255)
                                      tgctATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGACCAAcggc                                                                                                                                              >  3:173225/1‑117 (MQ=255)
                                      tgctATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGACCAAcggc                                                                                                                                              <  4:173225/117‑1 (MQ=255)
                                         tATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGt                                                                                                                                                     <  8:59541/107‑1 (MQ=255)
                                         tATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGt                                                                                                                                                     >  7:59541/1‑107 (MQ=255)
                                           tCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGccc                                                                                                           <  2:283531/147‑1 (MQ=255)
                                           tCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACCTCTTCGGCCTGCAGGCCCtt                                                                                                         >  8:345111/1‑149 (MQ=255)
                                                  cAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCgg                                                                                                  <  8:60879/149‑1 (MQ=255)
                                                   aaGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGt                                                                                                 <  7:351799/149‑1 (MQ=255)
                                                             ccACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCg                                                                                        >  6:116884/1‑148 (MQ=255)
                                                                 aGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAAc                                                                                   <  5:353051/149‑1 (MQ=255)
                                                                   gTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACt                                                                                  >  7:51245/1‑148 (MQ=255)
                                                                        atgCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCAcgc                                                                            >  6:271074/1‑149 (MQ=255)
                                                                          gctgcCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACgcgc                                                                          >  7:2439/1‑149 (MQ=255)
                                                                              ccTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTg                                                                        <  5:271074/147‑1 (MQ=255)
                                                                                   aGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTc                                                                  <  8:51245/148‑1 (MQ=255)
                                                                                    ggTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTcg                                                                                                     <  6:335933/112‑1 (MQ=255)
                                                                                    ggTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTcg                                                                                                     >  5:335933/1‑112 (MQ=255)
                                                                                        cagcaCGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCa                                                             <  6:229595/148‑1 (MQ=255)
                                                                                              ggCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCg                                                      <  2:162019/149‑1 (MQ=255)
                                                                                                 ccTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCATTCGACCAGAATGCGGTg                                                   >  8:144185/1‑149 (MQ=255)
                                                                                                        ggCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCa                                            <  2:292882/149‑1 (MQ=255)
                                                                                                                     cGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCTCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGccc                               <  7:345111/149‑1 (MQ=255)
                                                                                                                      gTTGATAGCGTCGGCCGGTCAGCCGTTGGTNAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGccc                               >  1:29310/1‑148 (MQ=255)
                                                                                                                          aTAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGaaa                           <  7:110040/148‑1 (MQ=255)
                                                                                                                           tAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATg                                   >  5:364114/1‑139 (MQ=255)
                                                                                                                           tAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATg                                   <  6:364114/139‑1 (MQ=255)
                                                                                                                             gCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTg                        >  6:349028/1‑148 (MQ=255)
                                                                                                                               gTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCAc                     >  3:339457/1‑149 (MQ=255)
                                                                                                                                          cgCCGTTGACCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCAt          >  4:179405/2‑149 (MQ=255)
                                                                                                                                          aGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCAt          >  2:264272/1‑149 (MQ=255)
                                                                                                                                                  gTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAg  >  5:286621/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGAGGCCCCGCCTCCTCATTATTGAGTTGTTCCAGGCTGCTATCGCCGGCAAGCCAGCCCCCACAGGTATGATGCTGCCTTCAGGTGCAGCACGGCCCTTGTGGGCGCTGGTTTGCCGTTGATAGCGTCGGCCGGTCAGCCGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAG  >  NC_002947/1387749‑1388043

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