Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,753,641 A→C 11.1% V95V (GTT→GTG PP_3316 ← chaperone‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,753,6410AC11.1% 80.5 / 5.5 36V95V (GTT→GTGPP_3316chaperone‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base A (15/17);  new base C (2/2);  total (17/19)
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.09e-01

TGGCGGCTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATCACACAACAG  >  NC_002947/3753502‑3753790
                                                                                                                                           |                                                                                                                                                     
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                                                           <  2:272433/150‑1 (MQ=255)
                                                                                  cGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTgg                                                           <  1:654564/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)
                                                                                                                                          cAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATCACacaaca   <  1:443773/150‑1 (MQ=255)
                                                                                                                                           aaCGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGAtcatca          <  2:112587/142‑1 (MQ=255)
                                                                                                                                           aaCGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGAtcatca          >  1:112587/1‑142 (MQ=255)
                                                                                                                                           aaCGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATCACACAACAg  >  2:447025/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
TGGCGGCTTCCTGCAACAGGGCCTCGGAGTGTTCGCTGATGCGGCTGGCCAGCCCGCTGCTGCGGCGGCGTGGCGTGCCGGACGCTGCGCGGCCAGGTTCCCCAGACGACGCTTCAGCCACTTCATCGGCATCCGTGGCAACGGATGCCGAGTCTTCTCCGATGCGCAGCAAATCGCTGCCCATGAACTCTTCGAACAAACCGCTGCGCGAGCCGAAAAGGGCTTCCAGTGGGGAAAGGGTGCGTTTTTGCTGGCGCACCAATTGGCGGTAGTGATCATCACACAACAG  >  NC_002947/3753502‑3753790

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: