Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,899,253 T→G 10.3% I72I (ATA→ATC)  PP_4307 ← DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,899,2530TG10.3% 137.2 / 6.0 48I72I (ATA→ATC) PP_4307DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (25/18);  new base G (2/3);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAA  >  NC_002947/4899110‑4899396
                                                                                                                                               |                                                                                                                                               
aTCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGc                                                                                                                                            >  6:215578/1‑149 (MQ=255)
       cGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCGATAGCCAGcgg                                                                                                                                      <  8:128721/148‑3 (MQ=255)
        gCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCGATAGCCAGcggg                                                                                                                                     <  2:37696/148‑4 (MQ=255)
                                caccaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCg                                                                                                            >  7:289977/1‑149 (MQ=255)
                                  ccaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCggg                                                                                                          <  2:31236/149‑1 (MQ=255)
                                      cATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCa                                                                                                      <  5:109345/149‑1 (MQ=255)
                                            cGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGc                                                                                                <  4:293366/149‑1 (MQ=255)
                                                      ttGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACgcggc                                                                                      <  4:66274/149‑1 (MQ=255)
                                                      ttGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACacggc                                                                                      >  5:234808/1‑149 (MQ=255)
                                                       tGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACgcggcg                                                                                     >  4:324341/1‑149 (MQ=255)
                                                          accccgccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCa                                                                                   >  4:95810/4‑148 (MQ=255)
                                                           cacCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTCGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCaaa                                                                                 >  2:164079/1‑149 (MQ=255)
                                                              cgccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCg                                                                              >  2:679/1‑149 (MQ=255)
                                                                ccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGcc                                                                            <  5:215578/149‑1 (MQ=255)
                                                                             aaTCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACgg                                                               <  3:324341/149‑1 (MQ=255)
                                                                                    gCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCAc                                                                 >  4:44856/1‑140 (MQ=255)
                                                                                    gCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCAc                                                                 <  3:44856/140‑1 (MQ=255)
                                                                                    gCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGa                                                         >  3:14354/1‑148 (MQ=255)
                                                                                      ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                      >  3:333867/1‑149 (MQ=255)
                                                                                        gggCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACg                                                    >  8:30782/1‑149 (MQ=255)
                                                                                        gggCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACg                                                    >  8:176595/1‑149 (MQ=255)
                                                                                                gAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgcgc                                              >  6:177340/1‑147 (MQ=255)
                                                                                                   gCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGACTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGaa                                         >  1:241078/1‑149 (MQ=255)
                                                                                                      aacaacAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAacac                                      >  8:102477/1‑149 (MQ=255)
                                                                                                         aacaGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAAcacca                                    >  7:147408/1‑148 (MQ=255)
                                                                                                               agcCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAg                             <  5:177340/149‑1 (MQ=255)
                                                                                                                 ccAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCa                           >  6:144383/1‑149 (MQ=255)
                                                                                                                     ttGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCgg                       <  8:289977/149‑1 (MQ=255)
                                                                                                                      tGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCa                              <  4:20737/141‑1 (MQ=255)
                                                                                                                      tGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCa                              >  3:20737/1‑141 (MQ=255)
                                                                                                                      tGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGc                      >  8:258942/1‑149 (MQ=255)
                                                                                                                       gCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGc                      <  7:176595/148‑1 (MQ=255)
                                                                                                                         caaCAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTgg                   >  6:13349/3‑149 (MQ=255)
                                                                                                                         cGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTgg                   <  8:147408/149‑1 (MQ=255)
                                                                                                                         cGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTgg                   <  4:333867/149‑1 (MQ=255)
                                                                                                                           aCAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACACGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                 <  6:234808/149‑1 (MQ=255)
                                                                                                                             aaTGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCgg               <  7:30782/149‑1 (MQ=255)
                                                                                                                               tGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCTCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGa              >  2:77784/1‑148 (MQ=255)
                                                                                                                               tGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCAGTGCGGCAGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCCGCGACTCCGCAGCGCCAGAACACCCCATCCAGCCCCGGCTGCCCGGac             >  5:200817/1‑149 (MQ=255)
                                                                                                                                 gCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGacac           <  2:241078/149‑1 (MQ=255)
                                                                                                                                 gCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGacac           <  7:258942/149‑1 (MQ=255)
                                                                                                                                  ccACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACAcc          >  1:178404/1‑149 (MQ=255)
                                                                                                                                    aCCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCgg                       >  3:113450/1‑134 (MQ=255)
                                                                                                                                    aCCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCgg                       <  4:113450/134‑1 (MQ=255)
                                                                                                                                      cccGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCCCGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCa      >  5:334025/1‑149 (MQ=255)
                                                                                                                                       ttgCTGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                 <  8:106519/135‑1 (MQ=255)
                                                                                                                                       ttgCTGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                 >  7:106519/3‑137 (MQ=255)
                                                                                                                                          cTGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGaa  >  3:307923/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAA  >  NC_002947/4899110‑4899396

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: ‑