Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,857,758 A→C 11.2% E84A (GAG→GCG)  PP_5134 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,857,7580AC11.2% 78.4 / 2.2 33E84A (GAG→GCG) PP_5134hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/14);  new base C (2/2);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGG  >  NC_002947/5857620‑5857904
                                                                                                                                          |                                                                                                                                                  
gCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGTCCAGGCGCTCGCCTGGGCg                                                                                                                                         <  1:194777/150‑1 (MQ=255)
gCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACg                                                                                                                                         >  2:164733/1‑150 (MQ=255)
                     tgGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtt                                                                                                                    <  1:390704/150‑1 (MQ=255)
                            tGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCa                                                                                                             <  2:392852/150‑1 (MQ=255)
                              cACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAgg                                                                                                           <  1:18738/150‑1 (MQ=255)
                                      tCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGaa                                                                                                   <  2:476908/150‑1 (MQ=255)
                                            gACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAgcgc                                                                                               <  2:67710/148‑1 (MQ=255)
                                            gACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAgcgc                                                                                               >  1:67710/1‑148 (MQ=255)
                                                       tttGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGCGCGCCTGGACGGCCCGGCGCGGCTGGTGCGGTGTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGa                                                                                  >  2:185753/1‑150 (MQ=255)
                                                       tttGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGa                                                                                  >  2:403467/1‑150 (MQ=255)
                                                       tttGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGa                                                                                  >  2:298572/1‑150 (MQ=255)
                                                            tGGAGCTGTACGACACGCTCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCagg                                                                             >  2:162760/1‑150 (MQ=255)
                                                                    tACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAg                                                                     >  1:213147/1‑150 (MQ=255)
                                                                            gCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTgcg                                                             >  1:7853/1‑150 (MQ=255)
                                                                              ccGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGcc                                                           >  1:255449/1‑150 (MQ=255)
                                                                                  ccTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAAct                                                       <  2:159580/150‑1 (MQ=255)
                                                                                       aaGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTgg                                                  <  1:162760/150‑1 (MQ=255)
                                                                                          aTGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCa                                               >  2:326576/1‑150 (MQ=255)
                                                                                           tGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAg                                              <  1:185753/150‑1 (MQ=255)
                                                                                                tCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCa                                                                                                             >  1:291348/1‑82 (MQ=255)
                                                                                                tCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCa                                                                                                             <  2:291348/82‑1 (MQ=255)
                                                                                                           gACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtttt                                                                                                                  >  1:439630/1‑66 (MQ=255)
                                                                                                           gACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGtttt                                                                                                                  <  2:439630/66‑1 (MQ=255)
                                                                                                                 cTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGAGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATc                        <  1:164733/150‑1 (MQ=255)
                                                                                                                      aCGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATgc                   <  2:255449/150‑1 (MQ=255)
                                                                                                                        gCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAAc                                                        >  2:466263/1‑111 (MQ=255)
                                                                                                                        gCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAAc                                                        <  1:466263/111‑1 (MQ=255)
                                                                                                                             gAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTc                                                                                                                 <  2:312291/49‑1 (MQ=255)
                                                                                                                             gAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTc                                                                                                                 >  1:312291/1‑49 (MQ=255)
                                                                                                                               gcgcCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGcc          >  1:364736/1‑150 (MQ=255)
                                                                                                                                 gtccAGGCGCTCGCCTGGGCGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTAt                              <  2:331299/126‑1 (MQ=255)
                                                                                                                                 gtccAGGCGCTCGCCTGGGCGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTAt                              >  1:331299/3‑128 (MQ=255)
                                                                                                                                       gCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCAAACAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTgg  >  1:36111/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GCCCCAGCACTACTCGCCATGTGGAAGGTGCACGGCTGTCGGTGGACGGTGAACCTTTGCTGGAGCTGTACGACACGCCCGGCCTGGAAGATGCCATCGCCCTGCTCGACTACCTTGAACGCCTGGAGCGCCCAGGCGAGCGCCTGGACGGCCCTGCCCGGCTGGAGCGTTTTCTGCAGGGCAGCGAAGCGCGCCAACGCTTCGAGCAGGAGGCCAAGGTATTGCGCCAACTGCTGGCCAGTCATGCCGGCCTGTATGTGATCGATGCGCGCGAGCCGGTGCTGG  >  NC_002947/5857620‑5857904

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