Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,422,357 G→A 8.9% L215L (CTG→CTA)  moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,422,3570GA 8.9% 125.9 / 2.4 45L215L (CTG→CTA) moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCGGCGACCTTGGGGCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTCCACCGGAGGCGTCTCGGTGGGC  >  NC_002947/2422223‑2422498
                                                                                                                                      |                                                                                                                                             
gcggcgACCTTGGGGCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTgg                                                                                                                                   <  2:156541/147‑1 (MQ=255)
  ggcgACCTTGGGGCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGCTCTACAACAGCAATCGCAGGCTGCTGGGCAGCTgggggc                                                                                                                               >  8:94193/1‑145 (MQ=255)
  ggcgACCTTGGGGCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGctgctg                                                                                                                                            <  4:12935/136‑1 (MQ=255)
  ggcgACCTTGGGGCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGctgctg                                                                                                                                            >  3:12935/1‑136 (MQ=255)
            gggCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTggg                                                                                                                     <  7:30842/149‑1 (MQ=255)
            gggCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTg                                                                                                                       >  5:281019/1‑147 (MQ=255)
            gggCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTg                                                                                                                       <  6:281019/147‑1 (MQ=255)
                aTGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGc                                                                                                                 <  4:154929/149‑1 (MQ=255)
                     cGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCg                                                                                                                >  6:104511/1‑145 (MQ=255)
                     cGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCg                                                                                                                <  5:104511/145‑1 (MQ=255)
                        cTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGt                                                                                                         >  3:74056/1‑149 (MQ=255)
                                         ggTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGc                                                                                        <  7:94193/149‑1 (MQ=255)
                                                 tGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGAc                                                                                <  4:74056/149‑1 (MQ=255)
                                                    cGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTg                                                                             <  5:244575/149‑1 (MQ=255)
                                                     ggTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTggg                                                                                                                     >  8:185952/1‑108 (MQ=255)
                                                     ggTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTggg                                                                                                                     <  7:185952/108‑1 (MQ=255)
                                                      gTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCgg                                                                                                <  3:124992/128‑1 (MQ=255)
                                                      gTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCgg                                                                                                >  4:124992/1‑128 (MQ=255)
                                                           gTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCa                                                                      >  1:35405/1‑149 (MQ=255)
                                                             ccACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTACTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCAACCTGCCAGAGGACCTGGACCGCACt                                                                    >  3:247481/1‑149 (MQ=255)
                                                                   gCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCgcc                                                              >  3:142310/1‑149 (MQ=255)
                                                                      acgaACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCgcctgc                                                           >  7:112854/1‑149 (MQ=255)
                                                                              gccgcgCCAGGACTGCCGCCCGGGCCTGGGCAGATCTTCACCCGCCATCTCAGGCTACGGGGCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGTCATCCTGCCAGATGACCTGGCCCGCACTCGCGCCTGCCTGGCt                                                     <  4:247481/142‑1 (MQ=255)
                                                                                 gAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAAc                                                 <  6:177754/148‑1 (MQ=255)
                                                                                   gCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTgg                                              <  4:116422/149‑1 (MQ=255)
                                                                                    ccAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTgg                                              >  2:102745/1‑148 (MQ=255)
                                                                                    ccAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTgg                                              <  1:102745/148‑1 (MQ=255)
                                                                                    ccAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGCTCAGCAGCTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTCAACTggg                                             >  2:532/1‑149 (MQ=255)
                                                                                      aGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTg                                           >  3:192720/1‑149 (MQ=255)
                                                                                      aGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTg                                           >  7:6384/1‑149 (MQ=255)
                                                                                           tGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTACGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCg                                        >  4:248857/1‑147 (MQ=255)
                                                                                               gCTCGGGCCTGGGCAGTTCTACAACAGCAATCGCTGGCTGCTGGTCAGGTGGTTGCAGCGGCTGGGTTTCGATGTGGTCGATGCCGGCATCCTGCCGGTTGACCTGGCCCGCACTCGCGCGTGCCTGGCTGAACTGGGTGGCGTAGAcc                                  <  1:532/149‑1 (MQ=255)
                                                                                               gCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGAcc                                  >  7:14305/1‑149 (MQ=255)
                                                                                                        ggggCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGGCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGt                         <  3:79705/148‑1 (MQ=255)
                                                                                                          ggCAGATCTACAACAGCAATCGCAGGCTACTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGcc                                                                          >  5:164865/1‑98 (MQ=255)
                                                                                                          ggCAGATCTACAACAGCAATCGCAGGCTACTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGcc                                                                          <  6:164865/98‑1 (MQ=255)
                                                                                                               aTCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTg                          >  5:267083/1‑141 (MQ=255)
                                                                                                               aTCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTg                          <  6:267083/141‑1 (MQ=255)
                                                                                                                 cTACAACAGCAATCGCAGGCTGCTTGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTCCACCGGa                 <  4:192720/148‑1 (MQ=255)
                                                                                                                    caacaGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTCCACCGGa                 <  8:6384/145‑1 (MQ=255)
                                                                                                                    caacaGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTACGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTCCACCgg                  <  3:248857/144‑1 (MQ=255)
                                                                                                                     aacaGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTgg                                              <  4:267615/115‑1 (MQ=255)
                                                                                                                     aacaGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTgg                                              >  3:267615/1‑115 (MQ=255)
                                                                                                                       caGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGAc                                   >  2:291962/1‑124 (MQ=255)
                                                                                                                       caGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGAc                                   <  1:291962/124‑1 (MQ=255)
                                                                                                                        aGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGCCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTCATTCTGTCCACCGGAGGCGtctc         >  5:77915/1‑149 (MQ=255)
                                                                                                                                 aGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTCCACCGGAGGCGTCTCGGTGGGc  >  5:207949/1‑147 (MQ=255)
                                                                                                                                     tgctgGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTc                        >  6:1846/1‑121 (MQ=255)
                                                                                                                                     tgctgGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTc                        <  5:1846/121‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                             
GCGGCGACCTTGGGGCATGGCCGTCTGGACGTGGTGCGCAAGGTAAAAGTGGCGGTGCTGTCCACCGGCGACGAACTGGTCGAGCCAGGCCTGCCGCTCGGGCCTGGGCAGATCTACAACAGCAATCGCAGGCTGCTGGTCAGCTGGTTGCAGCGGCTGGGTTGCGAAGTGGTCGATGCCGGCATCCTGCCGGATGACCTGGCCCGCACTCGCGCCTGCCTGGCTGAACTGGGTGGCGTAGACCTGATTCTGTCCACCGGAGGCGTCTCGGTGGGC  >  NC_002947/2422223‑2422498

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