Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,853,118 T→A 15.5% 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,1180TA15.5% 69.5 / 3.7 31I169N (ATT→AAT) alkBalpha‑ketoglutarate‑ and Fe(II)‑dependent dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (19/7);  new base A (3/2);  total (22/9)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCT  >  NC_002947/3852970‑3853258
                                                                                                                                                     |                                                                                                                                            
gTCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAt                                                                                                                                              >  3:18982/1‑149 (MQ=255)
   aaCCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTgc                                                                                                                                           >  7:130874/1‑149 (MQ=255)
   aaCCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTgc                                                                                                                                           >  3:51196/1‑149 (MQ=255)
      cATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCt                                                                                                                                        <  4:51196/149‑1 (MQ=255)
      cATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCt                                                                                                                                        >  1:210209/1‑149 (MQ=255)
           ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                      <  4:284510/146‑1 (MQ=255)
           ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGa                                                                                                                                      >  3:284510/1‑146 (MQ=255)
           ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGc                                                                                                                                    >  5:285081/1‑148 (MQ=255)
                cGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGccttgtc                                                                                                                              <  2:210209/149‑6 (MQ=255)
                             gAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCt                                                                                                                 >  5:338197/1‑149 (MQ=255)
                               gCCTGCACCAGGACCGGGTCGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGTGATGTGCTGGTCTgg                                                                                                               <  3:313281/149‑1 (MQ=255)
                                     aCCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGc                                                                                                         >  1:218976/1‑149 (MQ=255)
                                         ggACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCg                                                                                                        >  5:113436/1‑146 (MQ=255)
                                                                   aGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTa                                                                           >  3:159168/1‑149 (MQ=255)
                                                                            tgtcggtgtcgCTGGGCTTACCGGCGGTGTTGCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTg                                                                                                                >  1:38337/1‑104 (MQ=255)
                                                                                  tgtcgCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGc                                                            >  3:37478/1‑149 (MQ=255)
                                                                                    ggcgtGGGCTTCCGGGGGTTGTT‑CTTGTGGGGGGGCGTGCCGGGTCCCTACATGGCTCCGCTCATTGCGCTGACCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGt                                                          <  4:37478/145‑1 (MQ=255)
                                                                                       cTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGc                                                                                 >  7:327714/1‑123 (MQ=255)
                                                                                               aCCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGcc                                               >  5:143335/1‑149 (MQ=255)
                                                                                                   ggggTGTT‑CGGGTGGGGCGGCCGGCAGCGGGCCGACAAGGCTCAGCGCATGGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTggg                                           <  1:30568/147‑1 (MQ=255)
                                                                                                     ggTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGc                                          >  1:305864/1‑148 (MQ=255)
                                                                                                         tt‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcg                                     >  3:311240/1‑149 (MQ=255)
                                                                                                              tGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcggcg                                  >  4:151489/1‑148 (MQ=255)
                                                                                                               gTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGCCCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGa                                >  5:244569/1‑149 (MQ=255)
                                                                                                                        ccTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCCCGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGa                       >  7:88232/1‑149 (MQ=255)
                                                                                                                                       aatggCTCTGCGATTTGGGATTAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCgg        <  4:311240/146‑1 (MQ=255)
                                                                                                                                         aGGCTCCGCGATTTGGGGTTAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGt      <  6:285081/149‑1 (MQ=255)
                                                                                                                                            cTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGcc   >  4:27608/1‑149 (MQ=255)
                                                                                                                                            cTCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGcc   >  4:336273/1‑149 (MQ=255)
                                                                                                                                             tCAGCGCATTGCGCTGAGCAACGGCGATGTGCTGGTATGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCTCAAGGCGGGTGACt  >  4:49698/1‑149 (MQ=255)
                                                                                                                                             tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTgg                                                                                                               <  5:110049/40‑1 (MQ=255)
                                                                                                                                             tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTgg                                                                                                               >  6:110049/1‑40 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GTCAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGTT‑CCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCT  >  NC_002947/3852970‑3853258

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