Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,853,118 T→A 17.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,1180TA17.4% 49.0 / 4.1 23I169N (ATT→AAT) alkBalpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (15/4);  new base A (2/2);  total (17/6)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTG  >  NC_002947/3852971‑3853259
                                                                                                                                                   |                                                                                                                                             
tCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAtt                                                                                                                                              <  3:247601/149‑1 (MQ=255)
tCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAtt                                                                                                                                              >  5:143378/1‑149 (MQ=255)
          ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGc                                                                                                                                     <  2:351717/148‑1 (MQ=255)
           cTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGACTGCAGCGTGCCGACAAGGCTCAGCGCATAGCGCTGAGCCa                                                                                                                                   >  4:61610/1‑149 (MQ=255)
            tGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       <  7:155406/144‑1 (MQ=255)
            tGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       >  8:155406/1‑144 (MQ=255)
                    aCGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCt                                                                                                                                         >  1:39968/1‑134 (MQ=255)
                    aCGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCt                                                                                                                                         <  2:39968/134‑1 (MQ=255)
                                  gCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTg                                                                                                            >  6:179155/1‑149 (MQ=255)
                                    tccAGGACCGCTACGAGCAGGATTTCGGTCAGCCTATCGTGTCGGGGTCGCTGGGCTTACGGGCGGTGTTCCTGTTCGGCGGCCTGCTGGGTGCCGACAAGGCTCAGCGCATTGTGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGc                                                                                                          <  4:152880/148‑1 (MQ=255)
                                      cAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGTGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGc                                                                                                          >  3:152880/1‑147 (MQ=255)
                                           ccGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGAcc                                                                                                   >  2:97959/1‑149 (MQ=255)
                                                  gaGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTgcg                                                                                            >  5:4415/1‑149 (MQ=255)
                                                                    ccGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTg                                                                                                                      >  3:350470/1‑105 (MQ=255)
                                                                     cGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTg                                                                                                                      >  3:70220/1‑104 (MQ=255)
                                                                                   tcgCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGa                                                                        >  5:156086/1‑136 (MQ=255)
                                                                                                  ggggTGTGCGTGTTGGGCGGCCTGCTGCGGGCCCAAATGCCTCAGCGCTTGGGGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTggg                                            <  5:179155/147‑1 (MQ=255)
                                                                                                                  gcggcCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAAc                            >  3:46293/1‑149 (MQ=255)
                                                                                                                               tGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTgcgc                >  7:333252/1‑148 (MQ=255)
                                                                                                                                   gACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCACCCTGACCTTGCGCAAGGc           >  3:93129/1‑149 (MQ=255)
                                                                                                                                   gACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGc           >  1:184139/1‑149 (MQ=255)
                                                                                                                                   gACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCATGGGCGAAAGGCGGATCAACCTGACCTTGCGCAAGGc           >  3:249161/1‑149 (MQ=255)
                                                                                                                                          cTCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGcc    >  1:306378/1‑149 (MQ=255)
                                                                                                                                             aGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCt   <  2:306378/147‑1 (MQ=255)
                                                                                                                                             aGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTg  >  1:142674/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTG  >  NC_002947/3852971‑3853259

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: