Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,705,336 A→T 8.0% S13T (TCC→ACC)  PP_5007 ← polyhydroxyalkanoate granule‑associated protein GA2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,705,3360AT 8.0% 136.3 / 2.8 50S13T (TCC→ACC) PP_5007polyhydroxyalkanoate granule‑associated protein GA2
Reads supporting (aligned to +/- strand):  ref base A (27/19);  new base T (2/2);  total (29/21)
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.22e-01

TCGGCCTTCTCGCCATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGACGGAATACCCAGGCGATT  >  NC_002947/5705209‑5705481
                                                                                                                               |                                                                                                                                                 
tcGGCCTTCTCGCCATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGgtgt                                                                                                                               >  1:270880/1‑148 (MQ=255)
  ggCCTTCTCGCCATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGttct                                                                                                                            >  7:28846/1‑149 (MQ=255)
      ttCTCGCCATCCTTCACCAGCGAGTCGAACAACTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGAACTTGCGGGAGTATTTCACGAACCCGCCGACCCAGGAGATGCCTTCTTTTTCGGTGTtattct                                                                                                                         >  8:168076/1‑148 (MQ=255)
      ctctcCCCATCCTTCACCACCGACTCGAACAGCTTCCGGCCGTCCTTGTCGATCTTCGAATATTTACCCAGCCCCGCCATCCATTTCTTGCGGGAGTATTTCTCGATCCCGCCGACCCATGAGCTGCCTTCTTTTTCGGTGTtcttct                                                                                                                         <  7:168076/147‑1 (MQ=255)
           gCCATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTgc                                                                                                                      <  5:82933/146‑1 (MQ=255)
            ccATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTgccagc                                                                                                                  >  1:239646/1‑149 (MQ=255)
            ccATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGACAGATCTTGCGGGAGTATTTCTCGATCACGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTATTCTTgccagc                                                                                                                  >  7:219936/1‑149 (MQ=255)
                  ttCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCGTGTTTTTCGGTGTTCTTCTTGCCAGCCATCCt                                                                                                            >  8:25475/1‑149 (MQ=255)
                    cACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCt                                                                                                            <  7:72185/147‑1 (MQ=255)
                    cACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTg                                                                                                           <  8:28846/148‑1 (MQ=255)
                      ccAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGct                                                                                                         >  6:63930/1‑148 (MQ=255)
                       cAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTtcttct                                                                                                                         <  6:236596/131‑1 (MQ=255)
                       cAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTtcttct                                                                                                                         >  5:236596/1‑131 (MQ=255)
                         ccGAGTCGAACAGCTTCCGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTcc                                                                                                     <  8:230315/148‑1 (MQ=255)
                             gTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTAt                                                                                                 <  2:239646/149‑1 (MQ=255)
                                      cTTCGGGCCGTCCTGGTCGATATTAGAATATATACCCAGCCACGACAGCCAGATATTTCGGGAGTATTTCTCGATCCCGCAGACACAGGAGCTGACTTCTTTTTCGGTGTTATTCATGCCAGCCATCCTGATCTCCTTATGGttttttc                                                                                        >  2:237092/1‑145 (MQ=255)
                                       ttCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGTGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTgcg                                                                                       >  3:244659/1‑149 (MQ=255)
                                       ttCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGTGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTgcg                                                                                       >  3:248628/1‑149 (MQ=255)
                                       ttCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTcttctt                                                                                                                        >  6:147169/1‑116 (MQ=255)
                                       ttCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTcttctt                                                                                                                        <  5:147169/116‑1 (MQ=255)
                                          gggCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGa                                                                                    <  5:63930/149‑1 (MQ=255)
                                                tCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCt                                                                              >  1:42611/1‑149 (MQ=255)
                                                 ccTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCt                                                                              >  7:243140/1‑148 (MQ=255)
                                                     gTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAagc                                                                         <  8:243140/149‑1 (MQ=255)
                                                      tCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAagca                                                                        <  2:10321/149‑1 (MQ=255)
                                                        gATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGTGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAagcag                                                                       <  4:244659/148‑1 (MQ=255)
                                                        gATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGTGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAagcag                                                                       <  4:248628/148‑1 (MQ=255)
                                                         aTCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAagcag                                                                       <  8:169548/147‑1 (MQ=255)
                                                         aTCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCGTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAagcag                                                                       >  7:169548/1‑147 (MQ=255)
                                                                  tAGATACCCCGCCCCGCCAGCCAGCTCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACt                                                            <  8:219936/149‑1 (MQ=255)
                                                                  tAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACt                                                            >  3:41755/1‑149 (MQ=255)
                                                                     aTACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGt                                                          >  3:57009/1‑148 (MQ=255)
                                                                        cccAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCaa                                                      <  4:57009/149‑1 (MQ=255)
                                                                                     gccagATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACaa                                         <  2:270880/149‑1 (MQ=255)
                                                                                         gATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGcc                                     >  1:10543/1‑149 (MQ=255)
                                                                                             ttGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCaa                                                              <  4:113361/120‑1 (MQ=255)
                                                                                             ttGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCaa                                                              >  3:113361/1‑120 (MQ=255)
                                                                                                 gggAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGt                              >  3:261272/1‑148 (MQ=255)
                                                                                                 gggAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGt                              >  3:257177/1‑148 (MQ=255)
                                                                                                        tttCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGtt                        >  3:30892/1‑147 (MQ=255)
                                                                                                        tttCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGtt                        >  3:222482/1‑147 (MQ=255)
                                                                                                           ctcGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGAc                   >  4:225584/1‑149 (MQ=255)
                                                                                                           ctcGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGAc                   >  7:62502/1‑149 (MQ=255)
                                                                                                              gATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGACgg                 >  2:8045/1‑148 (MQ=255)
                                                                                                                      cgACCCAGTAGCTGCATTCTTTTTCGGTTTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGATCAAGCAGCTCATGCAACTGTTCAAGCTTGATAGACAACGCCTCAACGTCATGTTTAGACGGAATACCCAg        <  1:129208/149‑1 (MQ=255)
                                                                                                                        aCCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGACGGAATACCCAGGc      >  5:119118/1‑149 (MQ=255)
                                                                                                                            aGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGACGGAATACCCAGGCGAtt  <  1:237092/149‑1 (MQ=255)
                                                                                                                             ggAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGACGGAATACCCAGGCGAtt  >  1:136849/1‑148 (MQ=255)
                                                                                                                              gAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCaa                                                              >  6:85982/1‑87 (MQ=255)
                                                                                                                              gAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCaa                                                              <  5:85982/87‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                 
TCGGCCTTCTCGCCATCCTTCACCAGCGAGTCGAACAGCTTCGGGCCGTCCTGGTCGATCTTCGAATAGATACCCAGCCCCGCCAGCCAGATCTTGCGGGAGTATTTCTCGATCCCGCCGACCCAGGAGCTGCCTTCTTTTTCGGTGTTCTTCTTGCCAGCCATCCTGCTCTCCTTATGGTTTGTGCGCGACGCGCTCAAGCAGCTCATGCAACTGTTCAAGCTTGATCGACAACGCCTCAACGTCATGTTTAGACGGAATACCCAGGCGATT  >  NC_002947/5705209‑5705481

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