Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,862,978 G→A 10.8% G174G (GGG→GGA)  PP_5134 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,862,9780GA10.8% 106.8 / 3.8 37G174G (GGG→GGA) PP_5134hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (2/2);  total (16/21)
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.43e-01

CGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCA  >  NC_002947/5862845‑5863110
                                                                                                                                     |                                                                                                                                    
cGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCt                                                                                                                                    <  2:148014/136‑1 (MQ=255)
  gTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTc                                                                                                                                  <  2:387210/136‑1 (MQ=255)
  gTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTc                                                                                                                                  >  1:194903/1‑136 (MQ=255)
  gTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTc                                                                                                                                  <  2:25201/136‑1 (MQ=255)
  gTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTc                                                                                                                                  <  1:41891/136‑1 (MQ=255)
       ggCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCa                                                                                                                             >  2:206977/1‑136 (MQ=255)
        gCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAt                                                                                                                            <  2:94940/136‑1 (MQ=255)
        gCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAt                                                                                                                            <  2:327235/136‑1 (MQ=255)
          cAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTg                                                                                                                          <  2:409212/136‑1 (MQ=255)
           aaGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTgg                                                                                                                         >  2:39042/1‑136 (MQ=255)
                  gcgACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGtt                                                                                                                  <  1:397533/136‑1 (MQ=255)
                           tGGACGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGt                                                                                                         <  2:238382/136‑1 (MQ=255)
                             gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTgg                                                                                                       <  1:114313/136‑1 (MQ=255)
                                 tGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGTTGCTCAACTTCGTCGCCAGCAGCCAGCAGCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACTCATTGGTACGGTTCGACAGCGTGGcccc                                                                                                   <  1:1129/136‑1 (MQ=255)
                                   cTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                 <  2:68469/136‑1 (MQ=255)
                                      gCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCag                                                                                              >  1:95930/1‑136 (MQ=255)
                                      gCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCCGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCag                                                                                              >  1:392431/1‑136 (MQ=255)
                                                   aGCCGTTGCTGCCGGTGCCCAACTTCGTCGCCAGCAGCCAGCACCGCGCGCCTGAATGGCGTGAAGCCCATGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAgcg                                                                                 >  1:384977/1‑136 (MQ=255)
                                                                 gTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACACATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGaaa                                                                   <  1:120518/136‑1 (MQ=255)
                                                                                 ccagcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctgctgct                                                   >  2:214354/1‑136 (MQ=255)
                                                                                  cagcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGctgctgctg                                                  <  1:70539/136‑1 (MQ=255)
                                                                                            caCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgccc                                        <  1:214354/136‑1 (MQ=255)
                                                                                            caCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgccc                                        <  2:194903/136‑1 (MQ=255)
                                                                                               cgcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgcccgcc                                     <  2:278174/136‑1 (MQ=255)
                                                                                                  gAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGACTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCg                                                                                                  <  2:252152/72‑1 (MQ=255)
                                                                                                  gAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGACTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCg                                                                                                  <  2:250582/72‑1 (MQ=255)
                                                                                                  gAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGACTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCg                                                                                                  >  1:250582/1‑72 (MQ=255)
                                                                                                  gAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGACTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCg                                                                                                  >  1:252152/1‑72 (MQ=255)
                                                                                                    gCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGgcg                                >  1:330138/1‑136 (MQ=255)
                                                                                                                 gAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCAGCAACGGCTGa                   >  1:352057/1‑136 (MQ=255)
                                                                                                                   aGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                 >  1:22399/1‑136 (MQ=255)
                                                                                                                   aGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                 >  1:328746/1‑136 (MQ=255)
                                                                                                                   aGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                 >  1:16186/1‑136 (MQ=255)
                                                                                                                           cccGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACca         >  2:108841/1‑136 (MQ=255)
                                                                                                                               gCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcagcag     <  2:368648/136‑1 (MQ=255)
                                                                                                                               gCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcagcag     <  2:409794/136‑1 (MQ=255)
                                                                                                                                  tGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGca  >  1:182361/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCA  >  NC_002947/5862845‑5863110

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