Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,489,534 C→T 11.1% A177A (GCC→GCT PP_2184 → formate dehydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,489,5340CT11.1% 92.9 / 2.5 36A177A (GCC→GCTPP_2184formate dehydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base T (2/2);  total (16/20)
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

CGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGG  >  NC_002947/2489390‑2489670
                                                                                                                                                |                                                                                                                                        
cggcggTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCggg                                                                                                                                       >  6:120739/1‑148 (MQ=255)
   cgttttCGCAGGCCTGGAGGCCGGTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCTTCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGa                                                                                                                                   <  8:24896/146‑1 (MQ=255)
    ggTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGaa                                                                                                                                  <  7:285885/149‑1 (MQ=255)
         cGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTAc                                                                                                                              <  4:218360/148‑1 (MQ=255)
          gCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACgt                                                                                                                            >  6:286122/1‑149 (MQ=255)
             ggCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACgtggt                                                                                                                         <  2:275597/149‑1 (MQ=255)
                             tgGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAg                                                                                                         <  2:23987/149‑1 (MQ=255)
                                              gAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAgg                                                                                                        >  8:279496/1‑133 (MQ=255)
                                              gAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAgg                                                                                                        <  7:279496/133‑1 (MQ=255)
                                                 gTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGc                                                                                     <  5:286122/149‑1 (MQ=255)
                                                    gTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGa                                                                                  <  8:27263/149‑1 (MQ=255)
                                                           cAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCAt                                                                            >  2:128048/1‑148 (MQ=255)
                                                                 tCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCTGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATg                                                                     >  7:275833/1‑149 (MQ=255)
                                                                 tCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCTGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATg                                                                     >  7:181882/1‑149 (MQ=255)
                                                                      tCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTc                                                                                                 >  6:185833/1‑116 (MQ=255)
                                                                      tCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTc                                                                                                 <  5:185833/116‑1 (MQ=255)
                                                                        cGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAgg                                                                                                        >  3:171647/1‑107 (MQ=255)
                                                                        cGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAgg                                                                                                        <  4:171647/107‑1 (MQ=255)
                                                                                        ggtggtgCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                            >  6:162403/1‑135 (MQ=255)
                                                                                        ggtggtgCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                            <  5:162403/135‑1 (MQ=255)
                                                                                              gcggcgTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                            >  6:39058/1‑129 (MQ=255)
                                                                                              gcggcgTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                            <  5:39058/129‑1 (MQ=255)
                                                                                                   gTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGa                                   >  3:304476/1‑149 (MQ=255)
                                                                                                      cccGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCa                                >  4:131264/1‑149 (MQ=255)
                                                                                                        cggccggaATCAAGTGGCGTAGCGTGCGTGATGCTGGCGCCTGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCAtt                              <  3:42817/149‑1 (MQ=255)
                                                                                                         ggccggcATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTg                             >  4:162299/1‑149 (MQ=255)
                                                                                                          gccggcATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGc                            <  8:237381/149‑1 (MQ=255)
                                                                                                                   aaGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGc                   >  8:113256/1‑149 (MQ=255)
                                                                                                                        gCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCGGGCCTCGCCGTgg              >  7:61074/1‑149 (MQ=255)
                                                                                                                        gCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTgg              <  7:113256/149‑1 (MQ=255)
                                                                                                                        gCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTgg              <  8:61074/149‑1 (MQ=255)
                                                                                                                          gTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGt            >  8:146623/1‑149 (MQ=255)
                                                                                                                          gTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGt            <  3:131264/149‑1 (MQ=255)
                                                                                                                               gtgcgtgATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGa       <  5:120739/149‑1 (MQ=255)
                                                                                                                                    tgATGCTGGCGCTGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAgg  <  8:181882/149‑1 (MQ=255)
                                                                                                                                    tgATGCTGGCGCTGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAgg  <  8:275833/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
CGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGG  >  NC_002947/2489390‑2489670

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