Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,265,830 C→T 10.5% G722S (GGC→AGC)  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,8300CT10.5% 98.1 / 3.4 38G722S (GGC→AGC) quiPacyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base C (18/16);  new base T (2/2);  total (20/18)
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.57e-01

GCCGCTGCTGGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAG  >  NC_002947/1265696‑1265974
                                                                                                                                      |                                                                                                                                                 
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)
                                                         tGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTg                                                                          <  1:252890/150‑1 (MQ=255)
                                                           ccGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCa                                                                        <  1:23899/150‑1 (MQ=255)
                                                              aGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCt                                                                     <  1:797961/150‑1 (MQ=255)
                                                                                  aaGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCgt                                                 >  1:414145/1‑150 (MQ=255)
                                                                                           ccccAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGcc                                        >  2:505030/1‑150 (MQ=255)
                                                                                             ccAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTAc                                                                                >  1:113165/1‑109 (MQ=255)
                                                                                             ccAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTAc                                                                                <  2:113165/109‑1 (MQ=255)
                                                                                                      tGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGTCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                             <  1:643845/150‑1 (MQ=255)
                                                                                                      tGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                             <  1:505030/150‑1 (MQ=255)
                                                                                                           tgAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATAc                        >  2:577435/1‑150 (MQ=255)
                                                                                                            gAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCACAGCTACTGGGTGTCGTCTTTGTACTGCCAGGTTCGTCGGTCGGCCCCCAGCTGTGCGATAcc                       >  2:662883/1‑150 (MQ=255)
                                                                                                                gCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCgg                                                                                          <  1:664591/80‑1 (MQ=255)
                                                                                                                gCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCgg                                                                                          >  2:664591/1‑80 (MQ=255)
                                                                                                                gCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGcc                   >  1:372591/1‑150 (MQ=255)
                                                                                                                   aGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTcc                <  2:414145/150‑1 (MQ=255)
                                                                                                                         gTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTgcg          >  1:708821/1‑150 (MQ=255)
                                                                                                                             gttCGCCGCTTGCCGCAGGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGcgcccg      >  1:268462/3‑150 (MQ=255)
                                                                                                                                 gccgctTGCCGCAGGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAg  >  1:272380/1‑150 (MQ=255)
                                                                                                                                     ccTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                             >  1:646551/1‑119 (MQ=255)
                                                                                                                                     ccTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGc                             <  2:646551/119‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                 
GCCGCTGCTGGTCAGTATCTGCAACGGTTCGGCCTGGCCGAAGTCGACGATCATCCGTGCCGAGGCTGGCAGATGCGTGTCGAAGTCGCTGCCCCAGGCGAATGGTGTGAGGGCCAGGGTAGTGAAGTCGCCGCCTGCGGCAAGCGGCGAACGGCTGATCGCATCACCCAGGCCATGGTAGGCCGGTGCCGGCCAGCGGTACTGGTGCAGCTTGCCCCACTGCCAGGTGCGTCGGTCGGCCCCCAGCTGTGCGATACCGGCCTCCATTGCGCCCGCCAG  >  NC_002947/1265696‑1265974

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: