Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,853,118 T→A 21.9% 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,1180TA21.9% 53.2 / 13.1 32I169N (ATT→AAT) alkBalpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (16/8);  new base A (4/3);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

TCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCT  >  NC_002947/3852971‑3853258
                                                                                                                                                   |                                                                                                                                             
tCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAtt                                                                                                                                              <  4:12107/149‑1 (MQ=255)
  aaCCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATgc                                                                                                                                            >  3:26575/1‑149 (MQ=255)
   aCCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTgcg                                                                                                                                           >  4:201706/1‑149 (MQ=255)
       ttACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGCGCAGCCTTTCGGGCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       <  1:225424/149‑1 (MQ=255)
       ttACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                       <  5:34658/149‑1 (MQ=255)
         aCCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAg                                                                                                                                      <  4:26575/148‑1 (MQ=255)
          ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGcc                                                                                                                                    <  3:201706/149‑1 (MQ=255)
            tGCCAGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCCGCGTGCCGCCAAGGCTCCGCGCATTGCGCTGAGCCAc                                                                                                                                  >  6:50530/1‑149 (MQ=255)
            tGCCAGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCAc                                                                                                                                  >  6:54389/1‑149 (MQ=255)
              ccGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCAc                                                                                                                                  <  8:197954/147‑1 (MQ=255)
                        gcCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGCTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTg                                                                                                                      >  1:85598/1‑149 (MQ=255)
                              gCCTGCACCAGGACCGCGGCGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAg                                                                                                                                      >  1:166394/1‑127 (MQ=255)
                              gCCTGCACCAGGACCGCGGCGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAg                                                                                                                                      <  2:166394/127‑1 (MQ=255)
                                                     cAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGcctt                                                                                                                                  <  6:729/108‑3 (MQ=255)
                                                     cAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGcctt                                                                                                                                  >  5:729/1‑106 (MQ=255)
                                                                                              aCCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGAt                                                                      >  1:138284/1‑126 (MQ=255)
                                                                                              aCCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGcc                                               >  7:46972/1‑149 (MQ=255)
                                                                                                  gcggTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGc                                                        >  7:72604/1‑136 (MQ=255)
                                                                                                        ttCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGc                                                            >  1:19167/1‑126 (MQ=255)
                                                                                                          ccTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcggc                                   >  4:164215/1‑149 (MQ=255)
                                                                                                           cTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcggag                                  >  5:10116/1‑147 (MQ=255)
                                                                                                            tGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcggcg                                  >  7:106955/1‑148 (MQ=255)
                                                                                                                    ggcCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCg                                                                        >  3:171186/1‑102 (MQ=255)
                                                                                                                           cgCGTGGCGACATGGCTCCGCGCTGGGGGCGGCGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTg                  <  5:50530/148‑1 (MQ=255)
                                                                                                                             gggcccGACACGTCTCCGCGCTTTGCGGTGCCCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTg                  <  5:54389/143‑1 (MQ=255)
                                                                                                                                   gACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAgg           >  1:206162/1‑149 (MQ=255)
                                                                                                                                          cTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGcc   >  4:262351/1‑149 (MQ=255)
                                                                                                                                           tCAGCGCTTTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGGGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCgg        <  7:167311/143‑1 (MQ=255)
                                                                                                                                           tCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCcgc                                                  <  1:215162/101‑1 (MQ=255)
                                                                                                                                           tCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCcgc                                                  >  2:215162/1‑101 (MQ=255)
                                                                                                                                           tCAGCGCAATGCGCTGAGCCTCGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCt  >  4:6686/1‑149 (MQ=255)
                                                                                                                                           tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGc                                                        <  4:74512/95‑1 (MQ=255)
                                                                                                                                           tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGc                                                        >  3:74512/1‑95 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGC‑GTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCT  >  NC_002947/3852971‑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: ‑