Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 741,504 A→C 8.0% intergenic (+101/‑399) PP_t12 → / → PP_0635 tRNA‑Thr/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947741,5040AC 8.0% 127.6 / 4.4 50intergenic (+101/‑399)PP_t12/PP_0635tRNA‑Thr/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base A (25/21);  new base C (2/2);  total (27/23)
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.00e-01

CGCACTCGTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCG  >  NC_002947/741355‑741616
                                                                                                                                                     |                                                                                                                
cgcACTCGTAACGCGAAGGTCGCAGGTTCGCTTCCTGTCTCCGGCACCAGAATCCGTTCGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTaa                                                                                                                  <  1:351839/150‑1 (MQ=255)
           cgcgAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGc                                                                                                       <  2:44489/150‑1 (MQ=255)
               aaGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGGGTTGCGCGTAAGCGCAACCCGCTTgg                                                                                                   >  1:186781/1‑150 (MQ=255)
                aGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGt                                                                                                  >  1:815678/1‑150 (MQ=255)
                        ggTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGtggtgg                                                                                          >  2:664472/1‑150 (MQ=255)
                          ttCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGcc                                                                                        >  2:320289/1‑150 (MQ=255)
                                tCCTGTCTCCGGCACCAGAATCAGTTTCTCCAGGACGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAGCCCGCTTGGTTGTGGTGGCCAAGTgg                                                                                  <  1:605089/150‑1 (MQ=25)
                                 ccTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTggg                                                                                 <  1:1104725/150‑1 (MQ=255)
                                  cTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGTCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                                >  2:1135938/1‑150 (MQ=255)
                                  cTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                                <  1:753229/150‑1 (MQ=255)
                                         cGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGGGTTGCGCTTACGCGCAACgcgcca                                                                                                     <  1:942059/122‑3 (MQ=38)
                                         cGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGGGTTGCGCTTACGCGCAACgcgcca                                                                                                     >  2:942059/1‑120 (MQ=38)
                                         cGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                                         >  2:553555/1‑150 (MQ=255)
                                          ggCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                                        <  2:496758/150‑1 (MQ=255)
                                             aCCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCaa                                                                                                            >  1:1198543/1‑111 (MQ=255)
                                             aCCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCTTACGCGCaa                                                                                                            <  2:1198543/111‑1 (MQ=255)
                                             aCCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                                     <  2:216550/150‑1 (MQ=255)
                                               cAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                                   >  1:319216/1‑150 (MQ=255)
                                               cAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                                   >  1:741291/1‑150 (MQ=255)
                                                 gAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                                 <  2:741291/150‑1 (MQ=255)
                                                 gAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                                 <  1:1004300/150‑1 (MQ=255)
                                                 gAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                                 <  2:126682/150‑1 (MQ=255)
                                                 gAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                                 >  1:126682/1‑150 (MQ=255)
                                                      aGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGTCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                                            <  1:1135938/150‑1 (MQ=255)
                                                       gTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                           >  1:730955/1‑150 (MQ=255)
                                                         tGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                                         <  2:1076689/150‑1 (MQ=255)
                                                           gTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                                       <  2:319216/150‑1 (MQ=255)
                                                           gTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                                       >  1:320308/1‑150 (MQ=255)
                                                                tGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAccc                                                  <  1:664472/150‑1 (MQ=255)
                                                                  gCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCgg                                                >  1:705052/1‑150 (MQ=255)
                                                                   cGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTTCACACCCGGc                                               >  1:630604/1‑150 (MQ=255)
                                                                           cGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                                       <  1:824982/150‑1 (MQ=255)
                                                                                ccTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTTCACACCCGGCAAGGGGAAGTGtt                                  <  2:630604/150‑1 (MQ=255)
                                                                                 cTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAAccc                                                                                                         >  2:1193687/1‑78 (MQ=255)
                                                                                 cTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAAccc                                                                                                         <  1:1193687/78‑1 (MQ=255)
                                                                                 cTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAAccc                                                                                                         >  2:591135/1‑78 (MQ=255)
                                                                                 cTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAAccc                                                                                                         <  1:591135/78‑1 (MQ=255)
                                                                                      aCCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCgg                            >  1:513924/1‑150 (MQ=255)
                                                                                       ccAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa                           >  2:1198638/1‑150 (MQ=255)
                                                                                        cAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                          <  1:1197185/150‑1 (MQ=255)
                                                                                          aTAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGc                        <  1:553555/150‑1 (MQ=255)
                                                                                             gAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAg                     >  1:818348/1‑150 (MQ=255)
                                                                                                gCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt                  >  1:186743/1‑150 (MQ=39)
                                                                                                 cccGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                                          <  1:493047/93‑1 (MQ=255)
                                                                                                 cccGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                                          >  2:493047/1‑93 (MQ=255)
                                                                                                   cGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTc               >  1:959930/1‑150 (MQ=38)
                                                                                                      cTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTcgcg            >  1:433938/1‑150 (MQ=34)
                                                                                                         gtgCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGt         <  2:815678/150‑1 (MQ=32)
                                                                                                              ggCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACcg    >  1:879729/1‑150 (MQ=32)
                                                                                                                cTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACcgcg  >  2:660907/1‑150 (MQ=11)
                                                                                                                cTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACcgcg  >  1:693601/1‑150 (MQ=11)
                                                                                                                                                     |                                                                                                                
CGCACTCGTAACGCGAAGGTCGCAGGTTCGATTCCTGTCTCCGGCACCAGAATCAGTTGGTCTATGGCGATTTCTCGCTACCTGGCACCAATAGAAGCCCGCCTTGTGCGGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAAGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCG  >  NC_002947/741355‑741616

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