Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,456,653 A→T 8.3% K113* (AAG→TAG) ‡ leuS → leucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,456,6530AT 8.3% 142.7 / 3.9 48K113* (AAG→TAG) ‡leuSleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (26/18);  new base T (2/2);  total (28/20)
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.73e-01

CCGCTACCAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGAC  >  NC_002947/5456505‑5456802
                                                                                                                                                     |                                                                                                                                                     
ccGCTACCAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCaa                                                                                                                                                      >  1:270314/1‑150 (MQ=255)
ccGCTACCAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCaa                                                                                                                                                      >  1:375375/1‑150 (MQ=255)
       cAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctg                                                                                                                                               >  1:939154/1‑150 (MQ=255)
         gcgcATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctggg                                                                                                                                             >  2:146419/1‑150 (MQ=255)
              tgacaggCAACAACAGCCTACAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctgggcctgg                                                                                                                                        >  2:103276/6‑150 (MQ=255)
                   ggCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATc                                                                                                                                   <  2:1036393/150‑1 (MQ=255)
                      aaGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGAc                                                                                                                                >  2:837741/1‑150 (MQ=255)
                      aaGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGACGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAAAGCCTGGGCCTGGCCATCGAc                                                                                                                                >  2:384355/1‑150 (MQ=255)
                               cTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACgg                                                                                                                       >  1:677212/1‑149 (MQ=255)
                               cTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGt                                                                                                                       <  2:947034/150‑1 (MQ=255)
                                  cAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGaa                                                                                                                    <  2:375375/150‑1 (MQ=255)
                                       gATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTcac                                                                                                               >  2:1076143/1‑150 (MQ=255)
                                            gCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCt                                                                                                          >  1:343466/1‑150 (MQ=255)
                                              tGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGt                                                                                                                       <  2:586327/135‑1 (MQ=255)
                                              tGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGt                                                                                                                       >  1:586327/1‑135 (MQ=255)
                                                ggACGCTTTCGGCATGCCCGCGGAAAACGCGGCGGTGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACGGGGCACATGAAGTCACCACCTGcca                                                                                                      >  2:335458/1‑148 (MQ=255)
                                                 gACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACAGCGACTACAGGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCaaa                                                                                                     <  1:334295/150‑1 (MQ=255)
                                                      tttCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctgggc                                                                                                                                            >  2:848539/1‑106 (MQ=255)
                                                      tttCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctgggc                                                                                                                                            <  1:848539/106‑1 (MQ=255)
                                                                 ccGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGtt                                                                                     <  1:146419/150‑1 (MQ=255)
                                                                                cgATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                      <  2:1003044/150‑1 (MQ=255)
                                                                                 gATGAAAAACAACGCCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTCGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGtt                                                                     >  2:194207/1‑150 (MQ=255)
                                                                                    gAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCAc                                                                  >  1:754429/1‑150 (MQ=255)
                                                                                    gAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGCTCAc                                                                  <  1:837741/150‑1 (MQ=255)
                                                                                      aaaaCAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCAccc                                                                >  1:1029035/1‑150 (MQ=255)
                                                                                        aacaacGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGc                                                              >  1:656435/1‑150 (MQ=255)
                                                                                                  ccccGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGaaa                                                    >  2:597560/1‑150 (MQ=255)
                                                                                                   cccGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGaaaa                                                   >  2:347769/1‑150 (MQ=255)
                                                                                                     cGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCTATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAg                                                 <  1:245421/150‑1 (MQ=255)
                                                                                                     cGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCTATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAg                                                 >  2:245421/1‑150 (MQ=255)
                                                                                                          aaGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGcatc                                            <  1:347769/150‑1 (MQ=255)
                                                                                                                  gTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAcc                                                                                        >  2:59453/1‑98 (MQ=255)
                                                                                                                  gTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAcc                                                                                        <  1:59453/98‑1 (MQ=255)
                                                                                                                     cGAAAACATCGACTACATGAGGTGCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCaa                                 <  1:335458/150‑1 (MQ=255)
                                                                                                                     cGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGTCTGTTCGAAAAAGGCATCATCTACCGCaa                                 <  1:103276/150‑1 (MQ=255)
                                                                                                                       aaaaCATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGa                               >  2:601725/1‑150 (MQ=255)
                                                                                                                            ccagtcgatggccaggCCCAGGCTCTTGAGCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCa                                                                              <  2:953797‑M1/82‑1 (MQ=255)
                                                                                                                            ccagtcgatggccaggCCCAGGCTCTTGAGCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCa                                                                              >  1:953797‑M1/17‑98 (MQ=255)
                                                                                                                            aTCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGt                          >  1:975530/1‑150 (MQ=255)
                                                                                                                               gACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTAcc                       <  2:1029035/150‑1 (MQ=255)
                                                                                                                                 cTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGt                     <  1:194207/150‑1 (MQ=255)
                                                                                                                                       gAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTg               >  1:999152/1‑150 (MQ=255)
                                                                                                                                       gAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTg               >  1:656579/1‑150 (MQ=255)
                                                                                                                                       ccaggCCCAGGCTCTTGAGCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAcc                                                                                        <  2:1043954‑M1/72‑1 (MQ=255)
                                                                                                                                       ccaggCCCAGGCTCTTGAGCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAcc                                                                                        >  1:1043954‑M1/6‑77 (MQ=255)
                                                                                                                                           aCCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGAc           >  1:1111790/1‑150 (MQ=255)
                                                                                                                                             CCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCCCCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                      >  2:1032266‑M1/1‑89 (MQ=255)
                                                                                                                                             CCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                      <  1:1032266‑M1/89‑1 (MQ=255)
                                                                                                                                                GCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGc      >  2:270314‑M1/1‑150 (MQ=255)
                                                                                                                                                  cTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCcggcg     <  2:939154/150‑1 (MQ=255)
                                                                                                                                                     aaGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGAc  <  2:1111790/150‑1 (MQ=255)
                                                                                                                                                     aaGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGAc  <  1:601725/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCTACCAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGAC  >  NC_002947/5456505‑5456802

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: