Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,812,772 G→A 11.9% intergenic (+23/‑87) rplT → / → pheS 50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,812,7720GA11.9% 101.4 / 6.0 42intergenic (+23/‑87)rplT/pheS50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (2/3);  total (19/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.17e-01

GATTGCTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGA  >  NC_002947/2812629‑2812917
                                                                                                                                               |                                                                                                                                                 
gATTGCTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTa                                                                                                                                             <  1:773591/150‑1 (MQ=255)
     cTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCt                                                                                                                                        >  2:138121/1‑150 (MQ=255)
      tGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCcggcgg                                                                                                                                                  >  2:1426459/1‑139 (MQ=255)
      tGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCcggcgg                                                                                                                                                  <  1:1426459/139‑1 (MQ=255)
        gCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAGAGGTTAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGcc                                                                                                                                     >  1:1085928/1‑150 (MQ=255)
            gAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCGCCCGGCAGCGCTAGCACCGCCgggt                                                                                                                                 <  2:1037732/150‑2 (MQ=255)
            gAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgggt                                                                                                                                 <  2:125642/150‑2 (MQ=255)
                   gCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                          <  2:432173/150‑1 (MQ=255)
                   gCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAAc                                                                                                                          <  1:138121/150‑1 (MQ=255)
                     ttCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGc                                                                                                                                               >  2:290405/1‑127 (MQ=255)
                     ttCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGc                                                                                                                                               <  1:290405/127‑1 (MQ=255)
                     ttCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGt                                                                                                                        >  1:1344619/1‑150 (MQ=255)
                           cgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGa                                                                                                                  <  1:950571/150‑1 (MQ=255)
                              aaTCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAg                                                                                                               <  1:444822/150‑1 (MQ=255)
                                    ccGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAag                                                                                                         <  1:836251/150‑1 (MQ=255)
                                               tGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAggg                                                                                                             >  1:402173/1‑135 (MQ=255)
                                               tGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAggg                                                                                                             <  2:402173/135‑1 (MQ=255)
                                               tGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAg                                                                                              <  1:1420997/150‑1 (MQ=255)
                                               tGGCCGATCTTGCAGTGAAAGAAAAAGGGGATTTTACTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGCGCATTGAACGTCAGCGCTAGGCGAAGAGCCTTCaacg                                                                                              >  1:642591/1‑148 (MQ=255)
                                                 gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAaga                                                                                                        <  1:1421292/138‑1 (MQ=255)
                                                 gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAaga                                                                                                        >  2:1421292/1‑138 (MQ=255)
                                                 gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCtt                                                                                                   <  2:1024094/143‑1 (MQ=255)
                                                 gCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCtt                                                                                                   >  1:1024094/1‑143 (MQ=255)
                                                      tCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTcc                                                                                       <  1:1177082/150‑1 (MQ=255)
                                                      tCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgg                                                                                                                                   <  1:1285008/106‑1 (MQ=255)
                                                      tCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgg                                                                                                                                   >  2:1285008/1‑106 (MQ=255)
                                                               gAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGACGTTCTGGCTGAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATGGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                              <  2:38803/150‑1 (MQ=255)
                                                                   gAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTAt                                                                          >  2:615543/1‑150 (MQ=255)
                                                                   gAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTAt                                                                          >  2:1185258/1‑150 (MQ=255)
                                                                     aaaaGCGGCGTTTGCTGCGATTGTCGACAAAGCACTCGCCGTTCGGCCATAATGACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt                                                                        <  1:1027920/150‑1 (MQ=255)
                                                                                           gTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCt                                                  <  2:1344619/150‑1 (MQ=255)
                                                                                                      tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCt                                                                                                    >  2:1037103/1‑89 (MQ=255)
                                                                                                      tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCt                                                                                                    <  1:1037103/89‑1 (MQ=255)
                                                                                                      tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTg                                                                >  1:999636/1‑125 (MQ=255)
                                                                                                      tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTg                                                                <  2:999636/125‑1 (MQ=255)
                                                                                                          gCCGTTCTGGCTTAAGTACCCGCGACAATCCCGCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGtctc                                   >  1:852779/1‑150 (MQ=255)
                                                                                                                        aGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTAGTATATGGAGTCTGTACATGCATGACTGGCAGTCCCTTGTCTGGCAAGCCCTAGAgg                     <  2:1085928/150‑1 (MQ=255)
                                                                                                                             ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg                >  2:903/1‑150 (MQ=255)
                                                                                                                                 gACAATCACACGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAc            <  1:1412987/150‑1 (MQ=255)
                                                                                                                                    aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGCGTCCGTACAATCTTAAGCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg         >  2:749814/1‑150 (MQ=255)
                                                                                                                                    aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCAATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg         >  2:1034347/1‑150 (MQ=255)
                                                                                                                                      atgcccGGCAGCGCTAGCACCGCCGGGTATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCt       >  1:823685/4‑150 (MQ=255)
                                                                                                                                          ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                       <  1:717750/130‑1 (MQ=255)
                                                                                                                                          ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                       >  2:717750/1‑130 (MQ=255)
                                                                                                                                           cggcggTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGa  <  1:1185258/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GATTGCTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGA  >  NC_002947/2812629‑2812917

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