Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,753,638 G→T 11.4% A96A (GCC→GCA PP_3316 ← chaperone‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,753,6380GT11.4% 98.8 / 2.8 35A96A (GCC→GCAPP_3316chaperone‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base T (2/2);  total (17/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 = 9.71e-01

TCCGTGGCACGCTTGGCGGCTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATCAC  >  NC_002947/3753489‑3753783
                                                                                                                                                     |                                                                                                                                                 
tCCGTGGCACGCTTGGCGGCTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTg                                                                                                                                                   >  1:272433/1‑150 (MQ=255)
             tGGCGGCTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCg                                                                                                                                      <  2:25822/150‑1 (MQ=255)
                 ggcTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAg                                                                                                                                    >  1:213820/1‑148 (MQ=255)
                 ggcTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAg                                                                                                                                    <  2:213820/148‑1 (MQ=255)
                  gcTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTGCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGtct                                                                                                                                 >  1:103977/1‑150 (MQ=255)
                   cTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCGCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCTCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGtctt                                                                                                                                <  2:102230/150‑1 (MQ=255)
                        tGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCg                                                                                                                           >  2:286433/1‑150 (MQ=255)
                           aaCAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATg                                                                                                                        >  2:788183/1‑150 (MQ=255)
                                   cTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCaa                                                                                                                >  1:687602/1‑150 (MQ=255)
                                   cTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGGTCCCCAGACGACGCTTAAGCCACTTCATCGGCATCCGTGGGAACGGATGCCGAGTCTTCTCCGATGCGCAGCaa                                                                                                                >  1:687601/1‑150 (MQ=255)
                                            ttCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGcc                                                                                                       >  2:654564/1‑150 (MQ=255)
                                            ttCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGcc                                                                                                       >  1:504778/1‑150 (MQ=255)
                                               gCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCAt                                                                                                    >  2:256169/1‑150 (MQ=255)
                                                  gATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAAc                                                                                                                                              <  2:852471/105‑1 (MQ=255)
                                                  gATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAAc                                                                                                                                              >  1:852471/1‑105 (MQ=255)
                                                               aGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTTGCCACGGATGcc                                                                                                                                       >  2:290283/1‑99 (MQ=255)
                                                               aGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTTGCCACGGATGcc                                                                                                                                       <  1:290283/99‑1 (MQ=255)
                                                                 cccGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTTGCCACGGATGc                                                                                                                                        >  1:249131/1‑96 (MQ=255)
                                                                 cccGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTTGCCACGGATGc                                                                                                                                        <  2:249131/96‑1 (MQ=255)
                                                                                   tggcgtgCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAg                                                                >  2:443773/1‑150 (MQ=255)
                                                                                    ggcgtgCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAACCAGCTGCTCATGAACTGTTCGAACAAACCGCTGCGCGAGCCGAAAAgg                                                               <  1:788183/150‑1 (MQ=255)
                                                                                         gCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTc                                                          <  1:402839/150‑1 (MQ=255)
                                                                                               cGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTgg                                                    <  1:654564/150‑1 (MQ=255)
                                                                                               cGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTgg                                                    <  2:272433/150‑1 (MQ=255)
                                                                                                  tGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGa                                                 >  1:131397/1‑150 (MQ=255)
                                                                                                    cgcgGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGCCCGGTGCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGaaa                                               <  1:581369/150‑1 (MQ=255)
                                                                                                          cAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGc                                         >  2:802195/1‑150 (MQ=255)
                                                                                                                    aGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTg                               <  2:103977/150‑1 (MQ=255)
                                                                                                                        gacgCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGgcgc                           >  2:459210/1‑150 (MQ=255)
                                                                                                                                  ccATTTCATCGGCTTACGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAGCCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGc                 <  2:687601/150‑1 (MQ=255)
                                                                                                                                  ccACTTCATCGGCATGCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGc                 <  2:687602/150‑1 (MQ=255)
                                                                                                                                        cATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGt           <  2:504778/150‑1 (MQ=255)
                                                                                                                                         aTCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAg                                              <  2:235325/114‑1 (MQ=255)
                                                                                                                                         aTCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAg                                              >  1:235325/1‑114 (MQ=255)
                                                                                                                                             gCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGatca      <  2:343445/150‑1 (MQ=255)
                                                                                                                                                 ccGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATcac  <  1:459210/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TCCGTGGCACGCTTGGCGGCTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATCAC  >  NC_002947/3753489‑3753783

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: