Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	863,900	T→A	9.3%	intergenic (+17/+19)	hemH → / ← uraA	ferrochelatase/uracil permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	863,900	0	T	A	9.3%	115.5 / 3.6	43	intergenic (+17/+19)	hemH/uraA	ferrochelatase/uracil permease
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base A (2/2);  total (21/22)
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

AGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGC  >  NC_002947/863757‑864043                                                                                                                                                |                                                                                                                                                 agagATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGc                                                                                                                                             <  7:259410/149‑1 (MQ=255)   agaTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGccc                                                                                                                                           >  2:23590/1‑149 (MQ=255)             gCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTGCGACCg                                                                                                                                 >  1:148557/1‑149 (MQ=255)                         gCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCa                                                                                                                     <  5:105222/149‑1 (MQ=255)                           aGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTTGCTGATATGTGTTAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAg                                                                                                                    <  7:43260/148‑1 (MQ=255)                             tttGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGctt                                                                                                                 >  7:83974/1‑149 (MQ=255)                                 tCGTGGCGGGCGGCCAGGTGCGGGGGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTTTTAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGcttctt                                                                                                              <  4:271535/148‑1 (MQ=255)                                 tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtg                                                                                                             >  2:280010/1‑149 (MQ=255)                                     ggcgggcgGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCt                                                                                                         >  6:168297/1‑149 (MQ=255)                                              ccAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCCCGTTTGGGCCCCGTACCACCCGGGATACAACCACCATCATGTGCTTCAAACTGt                                                                                                >  1:302462/1‑149 (MQ=255)                                                         tggtggtgcGGTGCCTGAATGATCACCGGGTGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGc                                                                                     <  6:189847/146‑1 (MQ=255)                                                                 ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg                                                                             <  8:83974/149‑1 (MQ=255)                                                                  cgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAgg                                                                            <  4:215022/149‑1 (MQ=255)                                                                              aTCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCaa                                                                  >  5:295086/1‑147 (MQ=255)                                                                               tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGTGGCCCCTTCCGACCGGTCATACATCCAGcttctt                                                                                                              <  8:240578/102‑1 (MQ=255)                                                                                  cccGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCca                                                            >  2:206280/1‑149 (MQ=255)                                                                                   ccGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCca                                                            >  4:70945/1‑148 (MQ=255)                                                                                    cGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCcac                                                           >  8:198028/1‑148 (MQ=255)                                                                                           ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATa                                                                <  5:173731/136‑1 (MQ=255)                                                                                           ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTACGACCGGTCATACATCCAGCTTATTGTGCTTACAGCTGTAATTACCCCGCAGAAACAGGTTCAGGGCaaaa                                                                >  6:173731/1‑134 (MQ=255)                                                                                              gAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTc                                                                                                                                       >  8:17707/1‑62 (MQ=37)                                                                                              gAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTc                                                                                                                                       <  7:17707/62‑1 (MQ=37)                                                                                                      tGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGAt                                        <  1:280010/148‑1 (MQ=255)                                                                                                        gCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGcc                                     >  2:211257/1‑149 (MQ=255)                                                                                                         cTGATATGTGTGAAACGGCCGGACTCCACGGGGCCGATTTGCGGCCCCGTCCGACCGGTCAGACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGccc                                    <  6:216523/149‑1 (MQ=255)                                                                                                             tatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCa                                <  7:198028/149‑1 (MQ=255)                                                                                                              atGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCa                                <  8:65227/148‑1 (MQ=255)                                                                                                                gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGTCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGgc                                                    <  4:247900/126‑1 (MQ=255)                                                                                                                gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGgc                                                    >  3:247900/1‑126 (MQ=255)                                                                                                                gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCa                                >  7:40025/1‑146 (MQ=255)                                                                                                                gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCa                                <  8:40025/146‑1 (MQ=255)                                                                                                                 tgtgAATAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGGTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGc                             <  3:196774/148‑1 (MQ=255)                                                                                                                  gtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCCGCCCCTTACGACCGGTAATAAATACAGCTTATTTTGCTTACAGCTGTAATTAACCGGCAGAAACAGGTTCAGGTCAAAAGCCACAA‑TGACGCACAGGGCGATACCCTTCAgtc                             >  4:196774/1‑145 (MQ=255)                                                                                                                     aaaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCGTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAg                        >  7:269130/1‑149 (MQ=255)                                                                                                                      aaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAg                        <  5:168297/148‑1 (MQ=255)                                                                                                                          gCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATTGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCAt                    <  1:206280/149‑1 (MQ=255)                                                                                                                            cTGAATCGAAGGGGCCGCTTTGCGGCCCCTTGCGACCGGGAAAACATCCAGCTTAATGTGCTTACAGCTAGCATTACCCGCGAGAAACAGGTTAAGCGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTACAGGCCCCACTCATCCg                 >  6:127094/1‑149 (MQ=255)                                                                                                                                aTCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGc              <  1:211257/148‑1 (MQ=255)                                                                                                                                  ggAAGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCa                                <  8:213378/127‑1 (MQ=255)                                                                                                                                  ggAAGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCa                                >  7:213378/2‑128 (MQ=255)                                                                                                                                         gCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCt    >  1:53364/1‑149 (MQ=255)                                                                                                                                          ccGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTg   <  8:35438/149‑1 (MQ=255)                                                                                                                                           cGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGTGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGc  >  8:276893/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                 AGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAA‑TGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGC  >  NC_002947/863757‑864043

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