Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,853,118 T→A 13.4% I169N (ATT→AAT)  alkB → alpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,853,1180TA13.4% 90.5 / 5.1 36I169N (ATT→AAT) alkBalpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (26/5);  new base A (3/2);  total (29/7)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCT  >  NC_002947/3852980‑3853258
                                                                                                                                          |                                                                                                                                             
aCCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       <  5:285297/147‑1 (MQ=255)
aCCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAg                                                                                                                                      <  7:51242/148‑1 (MQ=255)
aCCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGc                                                                                                                                     >  2:285571/1‑149 (MQ=255)
aCCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAg                                                                                                                                      <  3:70659/148‑1 (MQ=255)
  cTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       >  6:332613/1‑145 (MQ=255)
  cTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       <  5:332613/145‑1 (MQ=255)
        ggCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGa                                                                                                                             >  2:38915/1‑149 (MQ=255)
               gcCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTg                                                                                                                      >  3:57931/1‑149 (MQ=255)
                      ccTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTggg                                                                                                              >  6:281540/1‑149 (MQ=255)
                        tGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGt                                                                                                            >  6:253356/1‑149 (MQ=255)
                             cAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCg                                                                                                        >  3:27438/1‑148 (MQ=255)
                                  ccGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGAcc                                                                                                  >  5:194940/1‑149 (MQ=255)
                                  ccGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAg                                                                                                                                      <  6:224117/114‑1 (MQ=255)
                                  ccGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAg                                                                                                                                      >  5:224117/1‑114 (MQ=255)
                                      gacgaGCAGGATTTCGGTCAGCTGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCtt                                                                                              >  1:330935/1‑149 (MQ=255)
                                                          accGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTAc                                                                          >  6:319586/2‑149 (MQ=255)
                                                           ccGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTAcc                                                                         >  8:133137/1‑149 (MQ=255)
                                                             gATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCGGGGGTGGCGTAGACCGCTAGCGCTTCCACGGAGTGCTACCGa                                                                       >  8:316382/1‑149 (MQ=255)
                                                                    tcggtgtcgCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGc                                                                             >  1:247279/1‑136 (MQ=255)
                                                                      ggtgtcgCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCg                                                              >  6:201474/1‑149 (MQ=255)
                                                                                ggCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATcc                                                    >  7:220565/1‑149 (MQ=255)
                                                                                  cTTACCGGGGGTTTTCCTGTGGGGGGGCCTGCGGGGGGCCGACATGGATCTGCGATTTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGTGTGCTACCGATCAAACCCGGCGTGCATCcgc                                                  <  1:38915/149‑1 (MQ=255)
                                                                                   ttACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCcgcg                                                 >  3:103224/1‑149 (MQ=255)
                                                                                       cggcggTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTg                                             >  8:48947/1‑149 (MQ=255)
                                                                                       cggcggTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTg                                             >  6:332073/1‑149 (MQ=255)
                                                                                             tgtTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGaa                                       >  1:149265/1‑149 (MQ=255)
                                                                                               ttCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcg                                     >  3:143313/1‑149 (MQ=255)
                                                                                                                   gCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTgc                 >  1:78087/1‑149 (MQ=255)
                                                                                                                   gCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTgc                 >  3:4703/1‑149 (MQ=255)
                                                                                                                        ccGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTGCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATAAACATGACCTTccgc               >  1:259997/1‑147 (MQ=255)
                                                                                                                         cGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCt                    >  5:149916/1‑140 (MQ=255)
                                                                                                                         cGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAgg           >  6:139231/1‑149 (MQ=255)
                                                                                                                               ggCTCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTg     >  3:93497/1‑149 (MQ=255)
                                                                                                                                  tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGccg  >  7:255584/1‑148 (MQ=255)
                                                                                                                                   cAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGc    <  2:318494/146‑1 (MQ=255)
                                                                                                                                   cAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGc    >  1:318494/1‑146 (MQ=255)
                                                                                                                                          |                                                                                                                                             
ACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCT  >  NC_002947/3852980‑3853258

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: