Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,693,834 C→G 11.8% R87R (CGC→CGG)  PP_1489 → CheW domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,693,8340CG11.8% 94.1 / 5.7 34R87R (CGC→CGG) PP_1489CheW domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base G (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

AAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGA  >  NC_002947/1693686‑1693980
                                                                                                                                                    |                                                                                                                                                  
aaGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGcc                                                                                                                                                   >  2:167795/1‑150 (MQ=255)
aaGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGcc                                                                                                                                                   >  1:583953/1‑150 (MQ=255)
 aGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCct                                                                                                                                                  <  2:597344/150‑1 (MQ=255)
     ttGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtg                                                                                                                                              <  1:156946/150‑1 (MQ=255)
         cgcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctgg                                                                                                                                          >  2:753713/1‑150 (MQ=255)
         cgcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctgg                                                                                                                                          <  1:417441/150‑1 (MQ=255)
         cgcgTCAGCCGCTCAAGCAGCGCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGACGGCGGTGCAACGCACCAGCACCCGCctggtgctgg                                                                                                                                          <  2:411652/150‑1 (MQ=255)
                                        aCCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCt                                                                                                           >  1:256819/1‑150 (MQ=255)
                                           tgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGAt                                                                                                        >  2:110976/1‑150 (MQ=255)
                                                 ggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGa                                                                                                  <  2:583953/150‑1 (MQ=255)
                                                 ggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGa                                                                                                  <  2:256819/150‑1 (MQ=255)
                                                      cgGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAgg                                                                                             >  2:92032/1‑150 (MQ=255)
                                                           aTTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTTCCGCCTTGAGCTTTGGCCAGCCAGCGGAGCACGACACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGCTGCAGGCCAcc                                                                                        >  2:419035/1‑150 (MQ=255)
                                                                cGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCATGCCGTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCGTGGCCTGCTCCAGGATCAGGCCaagctgca                                                                                   <  1:732932/150‑8 (MQ=255)
                                                                    cATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTg                                                                               >  2:717186/1‑150 (MQ=255)
                                                                          ggCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGc                                                                         <  1:92032/150‑1 (MQ=255)
                                                                                 tGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCa                                                                                                                    <  1:578259/100‑1 (MQ=255)
                                                                                 tGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCa                                                                                                                    >  2:578259/1‑100 (MQ=255)
                                                                                        ggTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGATGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGtt                                                           >  1:145303/1‑150 (MQ=255)
                                                                                           cATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCa                                                        <  1:110976/150‑1 (MQ=255)
                                                                                                  cTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATcc                                                                                                      <  1:39424/97‑1 (MQ=255)
                                                                                                  cTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATcc                                                                                                      >  2:39424/1‑97 (MQ=255)
                                                                                                   tttccgccTTGAGCTTTGGCCAGCCGGCGGTGCCACACACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATg                                                <  2:140746/147‑1 (MQ=255)
                                                                                                   ttGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATg                                                >  2:263166/1‑150 (MQ=255)
                                                                                                          cTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGtt                                                           <  2:49562/132‑1 (MQ=255)
                                                                                                          cTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGtt                                                           >  1:49562/1‑132 (MQ=255)
                                                                                                           ttGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGac                                        >  2:361548/1‑150 (MQ=255)
                                                                                                                cTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAgg                                   <  1:361548/150‑1 (MQ=255)
                                                                                                                                 tGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTCTCTgg                  <  2:145303/150‑1 (MQ=255)
                                                                                                                                 tGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTgg                  <  1:717186/150‑1 (MQ=255)
                                                                                                                                  gCAACGCACCAGCACCAGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATTCTGGAGCAGGCCACCGAAACCCTGCGTTGACAACCGGACGAGTTCCAGCTCTATGGGGTGGACACAGGCGATGCGCCTTATCTggg                 >  1:516638/1‑150 (MQ=255)
                                                                                                                                      tgcaccagcaccaGGCGGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCg             >  1:139115/2‑150 (MQ=255)
                                                                                                                                       gCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCgg            >  1:220038/1‑150 (MQ=255)
                                                                                                                                         aCCAGCACCAGGCGGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTAt                                                 >  1:483653/1‑111 (MQ=255)
                                                                                                                                         aCCAGCACCAGGCGGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTAt                                                 <  2:483653/111‑1 (MQ=255)
                                                                                                                                                 ccGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGa  >  1:515882/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGA  >  NC_002947/1693686‑1693980

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