Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,103,089	G→C	9.3%	intergenic (+21/‑169)	phhB → / → PP_5694	pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,103,089	0	G	C	9.3%	131.0 / 3.6	43	intergenic (+21/‑169)	phhB/PP_5694	pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base C (2/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTC  >  NC_002947/5102946‑5103232                                                                                                                                                 |                                                                                                                                                actgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGc                                                                                                                                             <  7:328594/148‑1 (MQ=255)   tgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCTCTGGcccc                                                                                                                                          >  4:271912/1‑149 (MQ=255)     aaTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCTCTGGCCCCtt                                                                                                                                        <  1:221880/149‑1 (MQ=255)      aTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCTCTGGCCCCttc                                                                                                                                       <  3:271912/149‑2 (MQ=255)         gggCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                                     <  2:4209/148‑1 (MQ=255)           gCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCa                                                                                                                                  >  4:276789/1‑149 (MQ=255)            cAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTc                                                                                                                                   >  6:35494/1‑147 (MQ=255)            cAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTc                                                                                                                                   <  5:35494/147‑1 (MQ=255)                gTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATgggg                                                                                                                             <  2:76391/149‑1 (MQ=255)                 tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                                    >  1:130205/1‑141 (MQ=255)                 tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                                    <  2:130205/141‑1 (MQ=255)                 tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGc                                                                                                                            <  3:276789/149‑1 (MQ=255)                              gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCt                                                                                                                 <  1:177964/147‑1 (MQ=255)                                      aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                        >  3:69653/1‑148 (MQ=255)                                      aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                        >  7:297109/1‑148 (MQ=255)                                       cTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGc                                                                                                       >  4:321181/1‑148 (MQ=255)                                          gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCaaaa                                                                                                   <  3:321181/149‑1 (MQ=255)                                            tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgcg                                                                                                                   >  5:149498/1‑131 (MQ=255)                                            tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgcg                                                                                                                   <  6:149498/131‑1 (MQ=255)                                                         ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCgg                                                                                    >  8:23984/1‑149 (MQ=255)                                                             aaCGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGaa                                                                                 >  6:79388/1‑148 (MQ=255)                                                                          tgtgtgCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGt                                                                   <  7:23984/149‑1 (MQ=255)                                                                           gtgtgCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGt                                                                   >  8:207338/1‑148 (MQ=255)                                                                                cacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                               >  6:44723/1‑147 (MQ=255)                                                                                     acTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACtt                                                        <  7:207338/149‑1 (MQ=255)                                                                                        gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                      >  1:218181/1‑149 (MQ=255)                                                                                         aCAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGcc                                                                                                                                            <  2:245505/61‑1 (MQ=255)                                                                                         aCAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGcc                                                                                                                                            >  1:245505/1‑61 (MQ=255)                                                                                                   ggAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCg                                          >  7:281843/1‑149 (MQ=255)                                                                                                         ggCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                     >  2:34420/1‑148 (MQ=255)                                                                                                                       a‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAc                      <  4:69653/149‑1 (MQ=255)                                                                                                                            atcacCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                               <  8:80113/102‑1 (MQ=255)                                                                                                                            aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                     >  4:331605/1‑146 (MQ=255)                                                                                                                            aaCACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                     <  3:331605/146‑1 (MQ=255)                                                                                                                                ccAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTgg              <  8:281843/149‑1 (MQ=255)                                                                                                                                   gaaggaagGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcg            >  3:215773/1‑148 (MQ=255)                                                                                                                                      ggaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGa        <  5:79388/149‑1 (MQ=255)                                                                                                                                      ggaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGa        <  5:44723/149‑1 (MQ=255)                                                                                                                                       gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtttt                                   <  7:220861/121‑1 (MQ=255)                                                                                                                                       gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtttt                                   >  8:220861/1‑121 (MQ=255)                                                                                                                                       aaaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGAc       >  7:205259/2‑149 (MQ=255)                                                                                                                                          ggggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACttt    <  8:297109/149‑1 (MQ=255)                                                                                                                                            ggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtc  <  2:218181/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                ACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTC  >  NC_002947/5102946‑5103232

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑