Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 489,290 G→C 14.7% R111G (CGC→GGC)  surA ← peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947489,2900GC14.7% 83.5 / 5.6 33R111G (CGC→GGC) surApeptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base C (3/2);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGAATGCGCTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTCGGCGGCACGCC  >  NC_002947/489151‑489431
                                                                                                                                           |                                                                                                                                             
acCTGAATGCGCTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCGGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATg                                                                                                                                     <  1:546433/150‑1 (MQ=255)
  cTGAATGCGCTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGTCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGt                                                                                                                                   >  1:177351/1‑150 (MQ=255)
  cTGAATGCGCTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGt                                                                                                                                   >  1:158741/1‑150 (MQ=255)
     aaTGCGCTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGcc                                                                                                                                >  1:182706/1‑150 (MQ=255)
        gcgcTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAAt                                                                                                                             <  2:158741/150‑1 (MQ=255)
                                  cgcgGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGa                                                                                                   <  2:182706/150‑1 (MQ=255)
                                    cgGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATg                                                                                                 <  1:243415/150‑1 (MQ=255)
                                     ggCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCCCAGCGCAATggcat                                                                                                                                 >  2:453016/1‑114 (MQ=255)
                                     ggCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCCCAGCGCAATggcat                                                                                                                                 <  1:453016/117‑4 (MQ=255)
                                          atgatCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATg                                                                                                                            <  2:448366/117‑1 (MQ=255)
                                          atgatCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATg                                                                                                                            >  1:448366/1‑117 (MQ=255)
                                           tgatCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATAc                                                                                          <  2:321227/150‑1 (MQ=255)
                                                   tcGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAgcgc                                                                                  <  1:835349/150‑1 (MQ=255)
                                                   tcGCGCTTGACCTGCCCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAgcgc                                                                                  <  1:515959/150‑1 (MQ=255)
                                                                               cgAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTcga                                                      <  2:496903/150‑1 (MQ=255)
                                                                                    gACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATc                                                 >  2:803257/1‑150 (MQ=255)
                                                                                                             cgcgGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCa                        <  1:803257/150‑1 (MQ=255)
                                                                                                                        tCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAg                                      >  2:241849/1‑125 (MQ=255)
                                                                                                                        tCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAg                                      <  1:241849/125‑1 (MQ=255)
                                                                                                                        tCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTc             >  1:837478/1‑150 (MQ=255)
                                                                                                                          cAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCa                                                                                        >  2:580970/1‑73 (MQ=255)
                                                                                                                          cAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCa                                                                                        <  1:580970/73‑1 (MQ=255)
                                                                                                                          cAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTcgg           >  2:473210/1‑150 (MQ=255)
                                                                                                                          cAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGACTGGTTCAGCTCTTCGTCGGTGATGCGGATAACAGAGCGCTCGCCGATCTTCAACTGCAGGGTATCGACGATCAGGCGTTACAGGACCTGCTGTTCCAGGGTGCTGGTagg           >  2:204354/1‑150 (MQ=255)
                                                                                                                           aGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCaga                                                                                      <  2:269178/74‑1 (MQ=255)
                                                                                                                           aGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCaga                                                                                      >  1:269178/1‑74 (MQ=255)
                                                                                                                            gCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTcggcg         <  2:446698/150‑1 (MQ=255)
                                                                                                                             cgctaTGGCATTGCGCTGGGCTATGGTGCAAATCGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGACCGCTCGCCGAACTGCAGCTGCAGGTTCTCGACCATCAGGCGTTCCCGGACCTGCTGTTCCAGGGCCCTGGTcggcgg        <  2:29014/146‑1 (MQ=255)
                                                                                                                                   accATTGCCCAGCGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTcggcg         <  2:89874/142‑1 (MQ=255)
                                                                                                                                   accATTGCCCAGCGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTcggcg         >  1:89874/2‑143 (MQ=255)
                                                                                                                                   accATTGCCCAGCGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTCGGCGGCAcgcc  >  2:570389/2‑150 (MQ=255)
                                                                                                                                       ttGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAgc                                                                                    <  2:112149/64‑1 (MQ=255)
                                                                                                                                       ttGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAgc                                                                                    >  1:112149/1‑64 (MQ=255)
                                                                                                                                           |                                                                                                                                             
ACCTGAATGCGCTCGGCAACGCGGCGTTGGCGCACGCGGCTGATGATCATCTCGCGCTTGACCTGCTCGCGAGCGTCGTCGAACGACAGGCCGTCACGGGCCAGCGCCGCGCGGAACTGGTCCAGCGACATGCCATTGCGCTGGGCAATGGTGCCAATGGCCTGGTTCAGCTCTTCGTCGGTGATGCGGATACCAGAGCGCTCGCCGATCTGCAGCTGCAGGTTCTCGACGATCAGGCGTTCCAGGACCTGCTGTTCCAGGGCGCTGGTCGGCGGCACGCC  >  NC_002947/489151‑489431

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: ‑