Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 795,629 G→A 19.0% A590T (GCC→ACC)  PP_0680 → ATP‑dependent protease‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947795,6290GA19.0% 43.3 / 4.9 21A590T (GCC→ACC) PP_0680ATP‑dependent protease‑like protein
Reads supporting (aligned to +/- strand):  ref base G (7/10);  new base A (2/2);  total (9/12)
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.34e-01

ACGGGTGCTGGACGACATGCTCGCCGGCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTCACCGGTTACCTGGG  >  NC_002947/795499‑795735
                                                                                                                                  |                                                                                                          
aCGGGTGCTGGACGACATGCTCGCCCGCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCAcc                                                                                                       >  1:405342/1‑136 (MQ=255)
          gacgacATGCTCGCCGGCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCcg                                                                                             >  2:452421/1‑136 (MQ=255)
              acATGCTCGCCGGCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCcggcgg                                                                                         <  2:269600/136‑1 (MQ=255)
                    tCGCCGGCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCgg                                                                                   <  1:149021/136‑1 (MQ=255)
                         ggCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTg                                                                              <  1:235868/136‑1 (MQ=255)
                                         acacGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCAATTATTCAGCGGTATTGTCGATATCGAggttta                                                              <  1:36627/136‑6 (MQ=255)
                                                      ggCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGc                                                 <  1:286675/136‑1 (MQ=255)
                                                       gCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGa                                                              <  2:461788/122‑1 (MQ=255)
                                                       gCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGa                                                              >  1:461788/1‑122 (MQ=255)
                                                            cGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCg                                           <  1:452421/136‑1 (MQ=255)
                                                              gCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGAt                                         <  2:194646/136‑1 (MQ=255)
                                                                     cAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCCACCTTGGCCAGCCGATCCACTcc                                  >  1:293248/1‑136 (MQ=255)
                                                                     cAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTcc                                  >  1:295470/1‑136 (MQ=255)
                                                                                ccGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGtgat                       >  2:447789/1‑136 (MQ=255)
                                                                                ccGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGACCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGtgat                       >  1:327813/1‑136 (MQ=255)
                                                                                ccGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGACCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGtgat                       >  1:325180/1‑136 (MQ=255)
                                                                                       tGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTc                >  1:405010/1‑136 (MQ=255)
                                                                                                  aCTCGGCGTTCGGCATGCCCGCGCGTATTTCGACCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTCACCGGTTACCt     <  2:325180/136‑1 (MQ=255)
                                                                                                  aCTCGGCGTTCGGCATGCCCGCGCGTATTTCGACCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTCACCGGTTACCt     <  2:327813/136‑1 (MQ=255)
                                                                                                    tCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTCACCGGTTACCCggg  <  2:389334/136‑1 (MQ=255)
                                                                                                     cGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTCACCGGTTACCTggg  <  2:405342/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                          
ACGGGTGCTGGACGACATGCTCGCCGGCATCATCCTGATCGACACGGAGGGCGCGGCCATCGGCAAGTGCAACGGGCTGACCGTGCTTGAAGTGGGTGACTCGGCGTTCGGCATGCCCGCGCGTATTTCGGCCACCGTCTACCCCGGCGGCAGCGGTATTGTCGATATCGAGCGTGAGGTCAACCTTGGCCAGCCGATCCACTCCAAAGGGGTGATGATCCTCACCGGTTACCTGGG  >  NC_002947/795499‑795735

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