Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,456,659 Δ1 bp 8.9% coding (343/2607 nt) leuS → leucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,456,6590C. 8.9% 129.1 / 8.1 45coding (343/2607 nt)leuSleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base . (2/2);  total (24/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 = 3.30e-01

CAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGAC  >  NC_002947/5456512‑5456802
                                                                                                                                                    |                                                                                                                                               
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                                                                                                                       >  1:586327/1‑135 (MQ=255)
                                       tGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGt                                                                                                                       <  2:586327/135‑1 (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)
                                                                                                                                                    |                                                                                                                                               
CAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGAC  >  NC_002947/5456512‑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: