Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,527 A→G 10.8% intergenic (+19/‑73) PP_1080 → / → grxD molybdopterin‑binding oxidoreductase/monothiol glutaredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,241,5270AG10.8% 96.0 / 3.4 37intergenic (+19/‑73)PP_1080/grxDmolybdopterin‑binding oxidoreductase/monothiol glutaredoxin
Reads supporting (aligned to +/- strand):  ref base A (18/15);  new base G (2/2);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

ACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGCC  >  NC_002947/1241384‑1241670
                                                                                                                                                |                                                                                                                                               
aCCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGgcg                                                                                                                                            >  1:68040/1‑149 (MQ=255)
           gCCGTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGTGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGCcttgc                                                                                                                                 <  8:212207/149‑3 (MQ=255)
           gCCCTGGCGCAACCCGGCGTGAGTGCCAATGAGCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTc                                                                                                                                  <  4:30641/148‑1 (MQ=255)
                   gcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATa                                                                                                                             >  8:215394/1‑145 (MQ=255)
                   gcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATa                                                                                                                             <  7:215394/145‑1 (MQ=255)
                                       aTGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACt                                                                                                     <  5:39377/149‑1 (MQ=255)
                                               accGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGa                                                                                             >  8:74543/1‑149 (MQ=255)
                                                         gcTTTCGGGACCAGGTTTCCGGAACCGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                   <  1:28854/149‑1 (MQ=255)
                                                           tttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                   >  4:254283/1‑147 (MQ=255)
                                                           tttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAAGCCGAGCGCCGTGCTCGGCctt                                                                                                                                   <  6:151934/99‑1 (MQ=38)
                                                           tttCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAAGCCGAGCGCCGTGCTCGGCctt                                                                                                                                   >  5:151934/1‑99 (MQ=38)
                                                                    cAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCTGATAGTTCAGCCGAGGTGCTACa                                                                        <  2:68040/149‑1 (MQ=255)
                                                                     aGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACAt                                                                       >  1:47811/1‑149 (MQ=255)
                                                                       gTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATg                                                                      >  5:6082/1‑148 (MQ=255)
                                                                        tttCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGa                                                                    >  1:268607/1‑149 (MQ=255)
                                                                             ggAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGTAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGatat                                                                 >  1:215541/1‑148 (MQ=255)
                                                                                       cgcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGAt                                                     <  5:253684/149‑1 (MQ=255)
                                                                                          tCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCa                                                   >  8:146212/1‑148 (MQ=255)
                                                                                            aaTGGCGTACCCGTGCAGGTCGCGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGa                                                                    <  2:181707/129‑1 (MQ=255)
                                                                                            aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCa                                                   >  7:63992/1‑146 (MQ=255)
                                                                                              tGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAagag                                              <  3:254283/149‑1 (MQ=255)
                                                                                                gCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATCCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAagag                                              >  5:76826/1‑147 (MQ=255)
                                                                                                     cccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGAtt                                        <  6:6082/148‑1 (MQ=255)
                                                                                                       cGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGcc                                     <  2:268607/149‑1 (MQ=255)
                                                                                                         tGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCcaa                                   <  2:199011/149‑1 (MQ=255)
                                                                                                                      cGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAaga                                               <  1:170362/124‑1 (MQ=255)
                                                                                                                      cGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAaga                                               >  2:170362/1‑124 (MQ=255)
                                                                                                                       gCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGc                     >  2:43158/1‑149 (MQ=255)
                                                                                                                            aGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTa                 >  5:198332/1‑148 (MQ=255)
                                                                                                                            aGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTAc                <  2:215541/149‑1 (MQ=255)
                                                                                                                                 tcggaaagCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGa            >  5:237212/8‑148 (MQ=255)
                                                                                                                                  ggCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGaaa          <  7:74543/149‑1 (MQ=255)
                                                                                                                                          cGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAag                                                >  4:100359/1‑103 (MQ=255)
                                                                                                                                          cGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAag                                                <  3:100359/103‑1 (MQ=255)
                                                                                                                                          cGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGTCTCGcc  >  4:32241/1‑149 (MQ=255)
                                                                                                                                          cGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGcc  >  3:9879/1‑149 (MQ=255)
                                                                                                                                          cGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGcc  >  3:228838/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGCC  >  NC_002947/1241384‑1241670

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: