Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,167 G→A 14.7% intergenic (+73/‑74) 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,1670GA14.7% 79.2 / 6.1 35intergenic (+73/‑74)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base A (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 = 7.76e-01

TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGC  >  NC_002947/5485022‑5485312
                                                                                                                                                 |                                                                                                                                                   
tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGCATAAGATGGTGGATTCAccgccc                                                                                                                                                  >  2:50684/1‑148 (MQ=255)
   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                                                                                                                                            >  3:127792/1‑103 (MQ=18)
                                                    aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATGCGGCGGTGAATCCACCATCGTATTcg                                                                                                                                            <  4:127792/103‑1 (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)
                                                                                                                                                 |                                                                                                                                                   
TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCC‑ACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCA‑CGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGC  >  NC_002947/5485022‑5485312

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