Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,150 A→G 26.2% intergenic (+56/‑91) 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,1500AG26.2% 62.6 / 23.3 42intergenic (+56/‑91)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base A (11/20);  new base G (5/6);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

GTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC  >  NC_002947/5485006‑5485297
                                                                                                                                                 |                                                                                                                                                       
gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCTCTAATTACGCTGATGCCGTTGCTTGCGCCGAACACTGTGGGAGCGGGCTTGACCGCGAATGCg                                                                                                                                                       <  2:46371/148‑1 (MQ=21)
 ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                    <  4:210743/149‑1 (MQ=255)
      tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa                                                                                                                                               <  5:147724/149‑1 (MQ=255)
          tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCa                                                                                                                                           >  7:177685/1‑149 (MQ=255)
              gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccat                                                                                                                                       >  1:68692/1‑147 (MQ=38)
                tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgccg                                                                                                                                     >  6:303516/1‑149 (MQ=255)
                tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCGTGCCCGCGAATACGATGGTGGATTAAccgccg                                                                                                                                     >  1:292825/1‑149 (MQ=255)
                 gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcg                                                                                                                                     >  1:91241/1‑144 (MQ=38)
                     tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATTGTGGATTCACCGCCG‑CATTc                                                                                                                               <  7:185316/149‑1 (MQ=255)
                     tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTc                                                                                                                               <  5:303516/149‑1 (MQ=255)
                         cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                           <  2:91241/149‑1 (MQ=37)
                          gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                           <  5:136809/148‑1 (MQ=37)
                            tgctgcCGACGAGGCCGGCATCGCGATGGCTTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCGTGGCCGCGAATAAGGTGGTGGGTTGACCCGCG‑CATTAGCGGGc                                                                                                                         >  7:256922/1‑148 (MQ=255)
                                 ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc                                                                                                                         >  4:22445/1‑143 (MQ=255)
                                 ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc                                                                                                                         <  3:22445/143‑1 (MQ=255)
                                     cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTc                                                                                                               <  3:215101/149‑1 (MQ=255)
                                                      tGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcgc                                                                                               <  8:177685/148‑1 (MQ=255)
                                                              aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCCCGAATAAGATGGGGGGTTAACCGCCG‑CCTTAGCGGGCAACCCCCCGCCCCCAGTGGTAGGCGGCA‑GCCA‑‑Cgg                                                                                   >  2:298663/1‑149 (MQ=255)
                                                                          cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                           <  7:20116/100‑1 (MQ=17)
                                                                          cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                           >  8:20116/1‑100 (MQ=17)
                                                                                  ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTcgc                                                                                                                             >  5:11485/1‑90 (MQ=255)
                                                                                  ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTcgc                                                                                                                             <  6:11485/90‑1 (MQ=255)
                                                                                             gcgTATCCCGTGGCAAGTCGAGATC‑CTGTGGGAGGGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                     <  7:219067/148‑1 (MQ=255)
                                                                                             gcgGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCGTTCCCCCGCATAACATGGTGGTTTAAACCCCGCCATAGCCGGGCAA‑CCCCCGCCCCCACAGTTGGGCGCCA‑GCCC‑‑CCCGGTCACATGAATACAGGGTTTTACaataaa                                                     >  8:219067/1‑143 (MQ=255)
                                                                                                       aGGCAAGTCGAGATC‑CTGTGGGCGCGGGCTTGCCCGCGAATACGAGGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat                                          <  8:256922/149‑1 (MQ=255)
                                                                                                          cAAGGCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCgg                                       <  3:241020/149‑1 (MQ=255)
                                                                                                                  gaTC‑CGGTGGGCGCGGGCTTGCCCGCGAATACTATGGTGTATTCACCGCCC‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg                               <  2:10098/149‑1 (MQ=39)
                                                                                                                    tC‑CTGTGGGCGCGGTGTTGCCCGATAATAGTATGGTGTATTCACCGCCT‑ATTTCGGGGGCATGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                              <  1:323597/148‑1 (MQ=255)
                                                                                                                         ttgGGCGCGGGTTTCCCCGCCAATAGGTTGTTGCTTTCACCGCCG‑CATTCCGGGGCATGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                          <  3:221770/147‑1 (MQ=255)
                                                                                                                          tgGGAGCGGGCTTGCCCGCTAATACGATGGTGTATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAACCCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                          <  1:112623/148‑1 (MQ=38)
                                                                                                                            ggAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCACCCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                        <  1:298663/149‑1 (MQ=255)
                                                                                                                              aGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGAGGGCAAGCCCGCTCCCACAGTGTTAGGCGCAA‑GCCA‑‑CGGGATCACTTGAAACGAGGTTTTTACATGAAAGTTTTGATAATAGGCAGCGGTGGGCGTGAGCAc                    >  4:50684/1‑149 (MQ=38)
                                                                                                                               gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGtttt                                                              <  2:7792/106‑1 (MQ=255)
                                                                                                                               gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGtttt                                                              >  1:7792/1‑106 (MQ=255)
                                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                        >  2:17295/1‑134 (MQ=25)
                                                                                                                                         ccGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                        <  1:17295/134‑1 (MQ=25)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAC‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                                     >  4:97707/1‑121 (MQ=255)
                                                                                                                                         ccGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                                     <  3:97707/121‑1 (MQ=255)
                                                                                                                                           gcgAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                              <  3:289308/126‑1 (MQ=25)
                                                                                                                                           gcgAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                              >  4:289308/1‑126 (MQ=25)
                                                                                                                                            ctaaTACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa      <  3:50684/147‑1 (MQ=255)
                                                                                                                                                tACGATGGTGGATCCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc  <  5:105081/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
GTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC  >  NC_002947/5485006‑5485297

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