Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,103,090	A→T	9.3%	intergenic (+22/‑168)	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,090	0	A	T	9.3%	111.5 / 3.6	43	intergenic (+22/‑168)	phhB/PP_5694	pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/21);  new base T (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 = 7.92e-01

ACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  NC_002947/5102946‑5103237                                                                                                                                                  |                                                                                                                                                    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                                                                                                                                         <  2:130205/141‑1 (MQ=255)                 tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttttt                                                                                                                                         >  1:130205/1‑141 (MQ=255)                 tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGc                                                                                                                                 <  3:276789/149‑1 (MQ=255)                              gtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCt                                                                                                                      <  1:177964/147‑1 (MQ=255)                                      aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                             >  7:297109/1‑148 (MQ=255)                                      aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                             >  3:69653/1‑148 (MQ=255)                                       cTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGc                                                                                                            >  4:321181/1‑148 (MQ=255)                                          gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCaaaa                                                                                                        <  3:321181/149‑1 (MQ=255)                                            tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgcg                                                                                                                        <  6:149498/131‑1 (MQ=255)                                            tCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgcg                                                                                                                        >  5:149498/1‑131 (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)                                                                                                                                                  agGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGAAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcactca  >  6:64016/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAA‑GTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  NC_002947/5102946‑5103237

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