Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,265,834 G→C 11.9% A720A (GCC→GCG)  quiP ← acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,265,8340GC11.9% 115.2 / 8.5 42A720A (GCC→GCG) quiPacyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base C (3/2);  total (22/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 = 5.25e-01

GCCGCTGCTGGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAGGCTGCGGGC  >  NC_002947/1265696‑1265983
                                                                                                                                          |                                                                                                                                                      
gCCGCTGCTGGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCTTGCCGCAGGCGGCGa                                                                                                                                             <  1:592004/150‑1 (MQ=255)
gCCGCTGCTGGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGa                                                                                                                                             >  2:196448/1‑150 (MQ=255)
     tgcttGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCATGCTAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTCCGGCAATCGGCGAACGGc                                                                                                                                        <  1:153793/150‑1 (MQ=255)
               tATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATc                                                                                                                              >  2:252890/1‑150 (MQ=255)
                  cTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCAcc                                                                                                                           <  1:196448/150‑1 (MQ=255)
                      aaCGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAgg                                                                                                                       <  1:842608/150‑1 (MQ=255)
                       aCGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGc                                                                                                                      <  2:830610/150‑1 (MQ=255)
                          tttCGGCCTGGCCAAAATCGACAATCATCCTTGCCGAGGCTGGCAGATGCGTGTCTAAGTCTCTGCCCCAGGCGAATGGTTTGAGGGCCAGGGTAGTGAAGTCCCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCAt                                                                                                                   <  1:221798/149‑1 (MQ=255)
                           ttCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATg                                                                                                                  <  1:107269/150‑1 (MQ=255)
                            tCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATgg                                                                                                                 <  2:197524/150‑1 (MQ=255)
                               gCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCa                                                                                                                    <  1:277409/144‑1 (MQ=255)
                               gCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCa                                                                                                                    >  2:277409/1‑144 (MQ=255)
                               accTAGCCGCAGTCAACGAACATCGATTCCGAGCCGGGCGGATTGGGGTCTAAGGGGTGGGGCCAGGGGAAGGGTGTGAGGGCCAGGGTAGTGAAGGCGCCGCCTGCGGCAAGCGGCGAACGGCTGGTCACATCACCCAGGCCATGGTAg                                                                                                              >  2:141763/2‑150 (MQ=255)
                                             gacgaTCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCTAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCa                                                                                                >  1:602832/1‑150 (MQ=255)
                                             gacgaTCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCa                                                                                                >  2:797961/1‑150 (MQ=255)
                                             gacgaTCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGCGCCGGCCa                                                                                                >  1:315541/1‑150 (MQ=255)
                                              acgaTCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAg                                                                                               >  1:427069/1‑150 (MQ=255)
                                                gaTCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCTTGCCGCAGGCGGCGactt                                                                                                                                          <  2:748285/105‑4 (MQ=255)
                                                gaTCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCTTGCCGCAGGCGGCGactt                                                                                                                                          >  1:748285/1‑102 (MQ=255)
                                                    atcCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTAc                                                                                         >  1:841154/1‑150 (MQ=255)
                                                    atcCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAAAGGCTGATCGCCTGCCCCATGCCATGGTAGGCCGGTTCGGGCCAGCGGTAc                                                                                         >  2:310047/1‑150 (MQ=255)
                                                         tGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTg                                                                                   <  1:252890/150‑1 (MQ=255)
                                                           ccGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCa                                                                                 <  1:23899/150‑1 (MQ=255)
                                                              aGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCt                                                                              <  1:797961/150‑1 (MQ=255)
                                                                                  aaGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCgt                                                          >  1:414145/1‑150 (MQ=255)
                                                                                           ccccAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGcc                                                 >  2:505030/1‑150 (MQ=255)
                                                                                             ccAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTAc                                                                                         >  1:113165/1‑109 (MQ=255)
                                                                                             ccAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTAc                                                                                         <  2:113165/109‑1 (MQ=255)
                                                                                                      tGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGTCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                                      <  1:643845/150‑1 (MQ=255)
                                                                                                      tGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                                      <  1:505030/150‑1 (MQ=255)
                                                                                                           tgAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATAc                                 >  2:577435/1‑150 (MQ=255)
                                                                                                            gAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCACAGCTACTGGGTGTCG‑TCTTTGTACTGCCAGGTTCGTCGGTCGGCCCCCAGCTGTGCGATAcc                                >  2:662883/1‑150 (MQ=255)
                                                                                                                gCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCgg                                                                                                   <  1:664591/80‑1 (MQ=255)
                                                                                                                gCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCgg                                                                                                   >  2:664591/1‑80 (MQ=255)
                                                                                                                gCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGcc                            >  1:372591/1‑150 (MQ=255)
                                                                                                                   aGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTcc                         <  2:414145/150‑1 (MQ=255)
                                                                                                                         gTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTgcg                   >  1:708821/1‑150 (MQ=255)
                                                                                                                             gttCGCCGCTTGCCGCAGGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGcgcccg               >  1:268462/3‑150 (MQ=255)
                                                                                                                                 gccgctTGCCGCAGGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAg           >  1:272380/1‑150 (MQ=255)
                                                                                                                                     ccTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                                      >  1:646551/1‑119 (MQ=255)
                                                                                                                                     ccTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                                      <  2:646551/119‑1 (MQ=255)
                                                                                                                                        gCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAGGCTGCgg    <  1:310047/150‑1 (MQ=255)
                                                                                                                                          ggCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAGGCTGCGGGc  >  1:563214/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                      
GCCGCTGCTGGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACT‑GGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAGGCTGCGGGC  >  NC_002947/1265696‑1265983

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