Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,159 G→A 25.0% 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,1590GA25.0% 48.6 / 13.2 38intergenic (+65/‑82)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (12/13);  new base A (4/4);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

ATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGA  >  NC_002947/5485011‑5485305
                                                                                                                                                    |                                                                                                                                                  
aTGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtgg                                                                                                                                                    <  1:171254/149‑1 (MQ=255)
  gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAAt                                                                                                                                                  >  3:29180/1‑149 (MQ=255)
        tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                            >  7:163428/1‑149 (MQ=255)
             aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                        <  8:163428/148‑1 (MQ=255)
               gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAtt                                                                                                                                     >  5:148830/1‑149 (MQ=255)
                  aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGGCCGCGCATACGAAGGTTGCATCACCCCCGCATTctc                                                                                                                                  >  8:167259/1‑147 (MQ=255)
                  aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                  >  1:72998/1‑149 (MQ=37)
                   tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                                 >  7:259962/1‑149 (MQ=37)
                     gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCTGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                <  8:259962/148‑1 (MQ=25)
                     gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCggg                                                                                                                               <  2:72998/149‑1 (MQ=25)
                     gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                               <  6:148830/149‑1 (MQ=255)
                     gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                <  6:66676/148‑1 (MQ=255)
                      ctgctgCCGACGAGGCCGGCATCTCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGc                                                                                                                              <  4:29180/149‑1 (MQ=25)
                              gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATAAGATGGTGGATTCACCGCCGCATTAGCGGGCCAGCCCGc                                                                                                                      >  7:72773/1‑149 (MQ=255)
                                      cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAAATCCTGTGGGAGCGGGCCTGGCCGCGAATAACATGGTGGGTTAAACGCCCCAATAGCGGGCAACCCCGCTACcacaa                                                                                                              >  8:174040/1‑148 (MQ=255)
                                                           tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCATTCCCGCCAAAAAGAAGGTGGATTAACCCCCGCAATAGCCCGCAACCACGCTACCCCACTGTTAGGCCCCAGCCACGGGa                                                                                         >  4:225540/1‑149 (MQ=255)
                                                               aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGGATAAGATGGGGGGTTGACCGCCGCATTCGCCGGCAAGCCCGCGCCCACAGGGTTCGGCGCAAGCCACGGGATGAc                                                                                     >  6:296637/1‑149 (MQ=255)
                                                                 gcccaaCTTCCGCCACAAATTACGCGGATCCCGCCGCATGCGGACACCCTGGGGCTGCGGGCGGGCCCGCTAATACCATGGTCGATCCCCCCCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                                   <  8:189955/144‑1 (MQ=255)
                                                                  cgcCACTTCCGCCACTAATTACTCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                <  2:213332/103‑1 (MQ=255)
                                                                  cgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                >  1:213332/1‑103 (MQ=255)
                                                                          ccGCCACTAATTACGCGTATCCCGAGGCAAGTCTCGATCCTGGGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCCCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                           <  5:135027/148‑1 (MQ=38)
                                                                           cGCCACTAATTCCGCGGATCCCGAGGCACGGCGAGAGCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                  <  1:209992/124‑1 (MQ=37)
                                                                            gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCACGCGAATAAGATGGTGGGTTCACCGCCGCCTTAGCGGGCAAGCCCCCGCCCACACTGGTCGGCGCAAGCCCCCGGATCACGTGAATGGAgggt                                                                        >  4:239799/1‑147 (MQ=255)
                                                                                cTAATTCCGCGGATCCCGAGGCAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                                    <  1:207298/149‑1 (MQ=255)
                                                                                 tAATTCCGCGTATCCCGTGGCATGGCGAGACCCTGTGGTAGGGGGCTGGCCCGCTAATACTATGGTGTATTCCCCGCCCCATGCGGGGCAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                   <  4:271087/149‑1 (MQ=255)
                                                                                       cgcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                             <  3:225540/149‑1 (MQ=255)
                                                                                             tCCCGTGGCAAGTCGAGATCGTGTGGGAGCGGGCTTGCCCGCTAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                        <  8:72773/148‑1 (MQ=255)
                                                                                               cctaGGAAAGTCGAGATCGTGTGGGAGCGGGCTTCCCGGCTAATACGATGGTGTATTCACCGCCGAATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                      <  7:167259/145‑1 (MQ=255)
                                                                                               ccGAGGCAAGTCGAGATCCTGTGGGAGCCGGCTTGCCCGCTAATACGATGGGGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                      <  6:302412/148‑1 (MQ=255)
                                                                                                   ggCAAGTCGATATCCTGGGGCAGCGGGCTTCCCGGACAATACTATGTTGAATTCACCGCCCAATTCGCGGTAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                                 <  3:51429/149‑1 (MQ=255)
                                                                                                                       gcgCGGGGTTGCCGGCTAATACGTTGGTGGATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                <  3:129838/144‑1 (MQ=255)
                                                                                                                            ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                       <  7:19575/118‑1 (MQ=255)
                                                                                                                            ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                       >  8:19575/1‑118 (MQ=255)
                                                                                                                            ggCTTGCCCGCCAATACGATGGTGCATTCACCGCCCCATTCCCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGccc                        <  3:22785/149‑1 (MQ=255)
                                                                                                                               ttGCCCGCGAATACGATGGTGGATTCACAGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                     >  7:176724/1‑149 (MQ=255)
                                                                                                                                  cccGCTAATACGATGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                  <  5:327752/149‑1 (MQ=255)
                                                                                                                                    cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                 >  3:304334/1‑148 (MQ=37)
                                                                                                                                    cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                 >  8:237297/1‑148 (MQ=255)
                                                                                                                                           ccGATGGGGCATGCCCCGCCGCATCCGCGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg         <  3:45024/148‑1 (MQ=255)
                                                                                                                                                  tcgATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa  <  5:117092/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGA  >  NC_002947/5485011‑5485305

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