Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,171 C→T 15.2% intergenic (+77/‑70) 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,1710CT15.2% 73.1 / 6.8 35intergenic (+77/‑70)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base T (2/3);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

AGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAG  >  NC_002947/5485025‑5485318
                                                                                                                                                  |                                                                                                                                                     
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)
                            agcgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTCCGCGGATCCCGAGCCAAGTCCAGATCCTGGGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGtt                                                                                                                        <  7:160363/147‑1 (MQ=255)
                                     gTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGATATCCTGTGGGAGCGGGCTTGCCCGCGAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCaa                                                                                                                <  1:94164/148‑1 (MQ=255)
                                          cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGCCAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGgc                                                                                                                    <  5:73315/139‑1 (MQ=255)
                                               gCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATc                                                                                                     <  6:304417/149‑1 (MQ=255)
                                                 aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATGCGGCGGTGAATCCACCATCGTATTcg                                                                                                                                                  <  4:127792/103‑1 (MQ=18)
                                                 aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATGCGGCGGTGAATCCACCATCGTATTcg                                                                                                                                                  >  3:127792/1‑103 (MQ=18)
                                                            ccGCCAATAATTCCGCGTATCCCGAGGCATGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTGGGTGTATTCCCCGCCGCTTTCGGGGGAAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                         <  4:38986/148‑1 (MQ=37)
                                                                         cgcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCGGCTAATACGTTGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                           <  5:224257/149‑1 (MQ=255)
                                                                                                      tgGGAGCGGGCTTGCCCGCGTATAGGTTGGTGGATTCACGGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGGTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  1:50684/149‑1 (MQ=39)
                                                                                                            cGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGc                                       <  3:248864/149‑1 (MQ=255)
                                                                                                             gggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGcc                                      >  4:143441/1‑149 (MQ=255)
                                                                                                              ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGcac                                     >  6:250845/1‑147 (MQ=255)
                                                                                                               gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCcct                                    >  6:313128/1‑149 (MQ=255)
                                                                                                                      cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCcctggcctg                              >  2:159915/1‑148 (MQ=255)
                                                                                                                       gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATTAAATTTTTGATAATCTGCAGCGGTGGCCGTGAGCA‑CGCcctggcctgg                             >  8:252233/1‑148 (MQ=25)
                                                                                                                          aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAg                         >  7:182315/1‑149 (MQ=255)
                                                                                                                            tACGATGGTGGATTCACCGCCGCATTCGCGGGAAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGTATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTc                       <  7:199899/149‑1 (MQ=255)
                                                                                                                               gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGcc                    >  6:289073/1‑149 (MQ=255)
                                                                                                                                 tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCa                  <  1:93307/149‑1 (MQ=255)
                                                                                                                                 cggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCa                  <  7:252233/148‑1 (MQ=38)
                                                                                                                                   gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAg                 >  5:289266/1‑148 (MQ=255)
                                                                                                                                      gATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACCGCCCTGGCCTGGAAAGTCGCCCAGGa               <  1:159915/148‑1 (MQ=255)
                                                                                                                                           aCCGCCGCATTCGCGGGCAAGCCCGCTCCCAACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcg         >  1:27610/1‑149 (MQ=255)
                                                                                                                                           aCCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgc        >  4:139460/1‑149 (MQ=255)
                                                                                                                                                  cATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCgag  >  8:50966/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAG  >  NC_002947/5485025‑5485318

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