Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,168 C→T 15.8% intergenic (+74/‑73) 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,1680CT15.8% 85.0 / 7.7 38intergenic (+74/‑73)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (12/20);  new base T (3/3);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

AAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTC  >  NC_002947/5485024‑5485315
                                                                                                                                                |                                                                                                                                                   
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                                                                 >  8:19575/1‑118 (MQ=255)
                                                                                                               ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                                 <  7:19575/118‑1 (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)
                                                                                                                                         ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCa        >  2:275171/1‑149 (MQ=255)
                                                                                                                                           cACCGCCGCTTGCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc       <  3:239799/148‑1 (MQ=255)
                                                                                                                                             ccgccgCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc                    <  7:95847/133‑1 (MQ=255)
                                                                                                                                             ccgccgCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc                    >  8:95847/1‑133 (MQ=255)
                                                                                                                                              cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGt   >  8:174570/1‑149 (MQ=255)
                                                                                                                                               gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc  <  7:174040/149‑1 (MQ=255)
                                                                                                                                               gccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc  >  5:150872/1‑149 (MQ=255)
                                                                                                                                                ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc  >  8:61509/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTC  >  NC_002947/5485024‑5485315

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: