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

CAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTG  >  NC_002947/3852972‑3853259
                                                                                                                                                   |                                                                                                                                             
cAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGCCCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTg                                                                                                                                             <  7:323278/149‑1 (MQ=255)
 aaCCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTg                                                                                                                                             <  3:254176/148‑1 (MQ=255)
   ccATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCTATgcgc                                                                                                                                          >  5:72761/1‑149 (MQ=255)
         ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTGCTTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCAATGCGCTGAGc                                                                                                                                     <  1:102485/149‑1 (MQ=255)
         ccTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCTATGCGCTGAGc                                                                                                                                     <  6:72761/148‑1 (MQ=255)
             ccGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACgg                                                                                                                                <  5:207208/149‑1 (MQ=255)
                ggCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGa                                                                                                                             >  1:84108/1‑149 (MQ=255)
                            aGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGAGCAGAGCCACGGCGATGTGCTGGTCTg                                                                                                                 >  2:184104/1‑149 (MQ=255)
                              ccTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTggg                                                                                                               >  5:100673/1‑149 (MQ=255)
                                  cACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTgg                                                                                                           >  7:19790/1‑149 (MQ=255)
                                             cgacgaGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCt                                                                                                >  2:142593/1‑149 (MQ=255)
                                                   gCAGGATTTCGGTCAGCCGATCGTGTCTGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGACTGCAGCTTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCt                                                                                          >  5:373159/1‑149 (MQ=255)
                                                                 aGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTa                                                                            >  7:244467/1‑149 (MQ=255)
                                                                   ccGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑GTCCTGTTAGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACCGCGATGTGCTGGTCTGGGCTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTAcc                                                                          >  7:367633/1‑149 (MQ=255)
                                                                           gtcggtgtcgCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAAc                                                                  >  1:168029/1‑149 (MQ=255)
                                                                           gtcggtgtcgCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAAc                                                                  >  7:75850/1‑149 (MQ=255)
                                                                             cggtgtcgCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGAGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACAGCTTGCGCTTCCACGGCGTCCTACCGATCAAAccc                                                                >  1:93117/1‑149 (MQ=255)
                                                                              ggtgtcgCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCg                                                               >  5:177573/1‑149 (MQ=255)
                                                                               gtgtcgCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCgg                                                              >  4:162158/1‑149 (MQ=255)
                                                                                     cTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCa                                                        >  8:307317/1‑149 (MQ=255)
                                                                                                         tCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAAcgg                                     >  8:258872/1‑149 (MQ=255)
                                                                                                           cTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACCGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGCCGAGCTACCGATCAACCCCGGCGTGCATCCGCGCCTGGGCGCAcggcg                                   >  5:299722/1‑149 (MQ=255)
                                                                                                             gTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGa                                 >  3:388936/1‑149 (MQ=255)
                                                                                                              tgcggcgcCCTGCCGCGGCCCCACACGCCTCAGCGCTTGGCGCTGCGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGa                                 <  2:84108/146‑1 (MQ=255)
                                                                                                                    ggcCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCt                          <  7:307317/149‑1 (MQ=255)
                                                                                                                     gcCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTg                         >  7:74883/1‑149 (MQ=255)
                                                                                                                        tGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCACACCCGGCGTGCATCCGCGCCTGGGCGCACGGCGGATCAACCTGAcc                      >  5:283153/1‑149 (MQ=255)
                                                                                                                        tGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCg                                                               >  3:127257/1‑108 (MQ=255)
                                                                                                                         gCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCt                     >  4:388831/1‑149 (MQ=255)
                                                                                                                              gTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACAGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTgcg                 >  8:59303/1‑148 (MQ=255)
                                                                                                                                   gACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGc           >  6:248059/1‑149 (MQ=255)
                                                                                                                                     cAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCgg         >  3:83149/1‑149 (MQ=255)
                                                                                                                                      aaGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCggg        >  7:152174/1‑149 (MQ=255)
                                                                                                                                      aaGGCGGAGCGCTTCGCGCGGTGCCACGGCGATGTGCTGGTCTGGGGTGGCGGAGACCGGTTGCGCTGCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCgg         <  8:367633/148‑1 (MQ=255)
                                                                                                                                         gcccAGCGCTTTGCGCTGACCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGc     <  6:177573/146‑1 (MQ=255)
                                                                                                                                         gCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGc     <  7:258872/149‑1 (MQ=255)
                                                                                                                                          cTCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGcc    >  8:321998/1‑149 (MQ=255)
                                                                                                                                           tctgcgcTTGGCGTTGCGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCt   <  8:185361/146‑1 (MQ=255)
                                                                                                                                           tCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGc                                                             >  3:369333/1‑91 (MQ=255)
                                                                                                                                           tCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGc                                                             <  4:369333/91‑1 (MQ=255)
                                                                                                                                           tCAGCGCAATGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCt   >  5:163553/1‑149 (MQ=255)
                                                                                                                                             aGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTg  >  6:281870/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CAACCATTACCTGCCGGGCACGCGCCTGAGCCTGCACCAGGACCGCGACGAGCAGGATTTCGGTCAGCCGATCGTGTCGGTGTCGCTGGGCTTACCGGCGGTG‑TTCCTGTTCGGCGGCCTGCAGCGTGCCGACAAGGCTCAGCGCATTGCGCTGAGCCACGGCGATGTGCTGGTCTGGGGTGGCGTAGACCGCTTGCGCTTCCACGGCGTGCTACCGATCAAACCCGGCGTGCATCCGCGCCTGGGCGAACGGCGGATCAACCTGACCTTGCGCAAGGCGGGTGCCTG  >  NC_002947/3852972‑3853259

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