Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,159 G→A 19.9% intergenic (+65/‑82) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,485,1590GA19.9% 69.3 / 11.6 36intergenic (+65/‑82)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base A (4/3);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

TGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGA  >  NC_002947/5485012‑5485305
                                                                                                                                                   |                                                                                                                                                   
tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa                                                                                                                                                    <  7:95351/149‑1 (MQ=255)
      aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAcc                                                                                                                                              <  3:68622/149‑1 (MQ=255)
          tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGCATAAGATGGTGGATTCAccgccc                                                                                                                                          >  2:50684/1‑148 (MQ=255)
           gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGGCCGCGAATAACATGGTGGATTAAcccacc                                                                                                                                          >  6:224257/1‑144 (MQ=255)
             aGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAccatcgtat                                                                                                                                       >  6:112470/1‑142 (MQ=38)
             aGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAt                                                                                                                                       <  6:229827/149‑1 (MQ=255)
                 aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                   <  5:112470/149‑1 (MQ=37)
                 aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                                        <  4:199576/144‑1 (MQ=255)
                 aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                                        >  3:199576/1‑144 (MQ=255)
                     ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                                >  5:304417/1‑148 (MQ=255)
                          gcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCGAAAAAGAAGGTGGATTAACCGCCGCAATAGCGGGCAAGc                                                                                                                           >  4:297647/1‑148 (MQ=255)
                          gcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                           <  3:297647/148‑1 (MQ=255)
                           ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATAAGATGGGGGGTTGACCGCCGCCATAGCGGGCCAGccc                                                                                                                         >  3:38986/1‑149 (MQ=255)
                           ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                           <  4:295392/147‑1 (MQ=255)
                           ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                           <  3:61924/147‑1 (MQ=255)
                                gaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGGCCGCGCATAAGATTGTGGGTTTCCCCCCGGCTTAGCGGGGCAGGCCGCTcg                                                                                                                    >  3:101972/1‑148 (MQ=255)
                                         agcgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTCCGCGGATCCCGAGCCAAGTCCAGATCCTGGGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGtt                                                                                                           <  7:160363/147‑1 (MQ=255)
                                                  gTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGATATCCTGTGGGAGCGGGCTTGCCCGCGAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                   <  1:94164/148‑1 (MQ=255)
                                                       cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGCCAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgc                                                                                                       <  5:73315/139‑1 (MQ=255)
                                                            gCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATc                                                                                        <  6:304417/149‑1 (MQ=255)
                                                             cATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAccat                                                                                                                                            <  4:174026/96‑3 (MQ=25)
                                                             cATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAccat                                                                                                                                            >  3:174026/1‑94 (MQ=25)
                                                              aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATGCGGCGGTGAATCCACCATCGTATTcg                                                                                                                                    <  4:127792/103‑1 (MQ=18)
                                                              aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATGCGGCGGTGAATCCACCATCGTATTcg                                                                                                                                    >  3:127792/1‑103 (MQ=18)
                                                                         ccGCCAATAATTCCGCGTATCCCGAGGCATGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTGGGTGTATTCCCCGCCGCTTTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                            <  4:38986/148‑1 (MQ=37)
                                                                                      cgcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCGGCTAATACGTTGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                              <  5:224257/149‑1 (MQ=255)
                                                                                                                   tgGGAGCGGGCTTGCCCGCGTATAGGTTGGTGGATTCACGGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGGTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                 <  1:50684/149‑1 (MQ=39)
                                                                                                                         cGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGc                          <  3:248864/149‑1 (MQ=255)
                                                                                                                          gggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGcc                         >  4:143441/1‑149 (MQ=255)
                                                                                                                           ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGcac                        >  6:250845/1‑147 (MQ=255)
                                                                                                                            gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCcct                       >  6:313128/1‑149 (MQ=255)
                                                                                                                                   cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCcctggcctg                 >  2:159915/1‑148 (MQ=255)
                                                                                                                                    gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATTAAATTTTTGATAATCTGCAGCGGTGGCCGTGAGCA‑CGCcctggcctgg                >  8:252233/1‑148 (MQ=25)
                                                                                                                                       aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAg            >  7:182315/1‑149 (MQ=255)
                                                                                                                                         tACGATGGTGGATTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGTATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTc          <  7:199899/149‑1 (MQ=255)
                                                                                                                                            gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGcc       >  6:289073/1‑149 (MQ=255)
                                                                                                                                              tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCa     <  1:93307/149‑1 (MQ=255)
                                                                                                                                              cggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCa     <  7:252233/148‑1 (MQ=38)
                                                                                                                                                gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAg    >  5:289266/1‑148 (MQ=255)
                                                                                                                                                   gATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACCGCCCTGGCCTGGAAAGTCGCCCAGGa  <  1:159915/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGA  >  NC_002947/5485012‑5485305

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