Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,485,159	G→A	28.1%	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_002947	5,485,159	0	G	A	28.1%	51.6 / 19.1	40	intergenic (+65/‑82)	purH/purD	bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base A (5/5);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-03

TGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAC  >  NC_002947/5485012‑5485306                                                                                                                                                    |                                                                                                                                                        tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa                                                                                                                                                        <  5:147724/149‑1 (MQ=255)     tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCa                                                                                                                                                    >  7:177685/1‑149 (MQ=255)         gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccat                                                                                                                                                >  1:68692/1‑147 (MQ=38)           tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgccg                                                                                                                                              >  6:303516/1‑149 (MQ=255)           tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGAATACGATGGTGGATTAAccgccg                                                                                                                                              >  1:292825/1‑149 (MQ=255)            gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcg                                                                                                                                              >  1:91241/1‑144 (MQ=38)                tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATTGTGGATTCACCGCCG‑CATTc                                                                                                                                        <  7:185316/149‑1 (MQ=255)                tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTc                                                                                                                                        <  5:303516/149‑1 (MQ=255)                    cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                                    <  2:91241/149‑1 (MQ=37)                     gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                                    <  5:136809/148‑1 (MQ=37)                       tgctgcCGACGAGGCCGGCATCGCGATGGCTTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGGCCGCGAATAAGGTGGTGGGTTGACCCGCG‑CATTAGCGGGc                                                                                                                                  >  7:256922/1‑148 (MQ=255)                            ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc                                                                                                                                  >  4:22445/1‑143 (MQ=255)                            ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc                                                                                                                                  <  3:22445/143‑1 (MQ=255)                                cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTc                                                                                                                        <  3:215101/149‑1 (MQ=255)                                                 tGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcgc                                                                                                        <  8:177685/148‑1 (MQ=255)                                                         aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAATAAGATGGGGGGTTAACCGCCG‑CCTTAGCGGGCAACCCCCCGCCCCCAGTGGTAGGCGGCA‑GCCA‑‑Cgg                                                                                            >  2:298663/1‑149 (MQ=255)                                                                     cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                                    <  7:20116/100‑1 (MQ=17)                                                                     cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg                                                                                                                                    >  8:20116/1‑100 (MQ=17)                                                                             ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTcgc                                                                                                                                      >  5:11485/1‑90 (MQ=255)                                                                             ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTcgc                                                                                                                                      <  6:11485/90‑1 (MQ=255)                                                                                        gcgTATCCCGTGGCAAGTCGAGATCCTGTGGGAGGGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                              <  7:219067/148‑1 (MQ=255)                                                                                        gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCATAACATGGTGGTTTAAACCCCGCCATAGCCGGGCAA‑CCCCCGCCCCCACAGTTGGGCGCCA‑GCCC‑‑CCCGGTCACATGAATACAGGGTTTTACaataaa                                                              >  8:219067/1‑143 (MQ=255)                                                                                                  aGGCAAGTCGAGATCCTGTGGGCGCGGGCTTGCCCGCGAATACGAGGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat                                                   <  8:256922/149‑1 (MQ=255)                                                                                                     cAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCgg                                                <  3:241020/149‑1 (MQ=255)                                                                                                             gaTCCGGTGGGCGCGGGCTTGCCCGCGAATACTATGGTGTATTCACCGCCC‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg                                        <  2:10098/149‑1 (MQ=39)                                                                                                               tCCTGTGGGCGCGGTGTTGCCCGATAATAGTATGGTGTATTCACCGCCT‑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)                                                                                                                                                   ggATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCCCAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc  >  6:266307/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAC  >  NC_002947/5485012‑5485306

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