Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	863,899	T→A	11.5%	intergenic (+16/+20)	hemH → / ← uraA	ferrochelatase/uracil permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	863,899	0	T	A	11.5%	132.0 / 8.8	52	intergenic (+16/+20)	hemH/uraA	ferrochelatase/uracil permease
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base A (3/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGC  >  NC_002947/863751‑864043                                                                                                                                                     |                                                                                                                                                  gCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCtt                                                                                                                                                   >  2:93310/1‑149 (MQ=255)     gAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCg                                                                                                                                               <  3:4625/149‑1 (MQ=255)        gagaTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGG‑cc                                                                                                                                           <  3:23062/149‑1 (MQ=255)             aGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTc                                                                                                                                      >  7:138351/1‑149 (MQ=255)               gCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGG‑CCCCTTc                                                                                                                                      <  3:64936/147‑1 (MQ=255)               gCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGG‑CCCCTTc                                                                                                                                      >  4:64936/1‑147 (MQ=255)                 aTGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGAc                                                                                                                                  >  5:231743/1‑149 (MQ=255)                             gAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATc                                                                                                                      >  8:204710/1‑149 (MQ=255)                             gAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATc                                                                                                                      <  1:93310/149‑1 (MQ=255)                              aGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATcc                                                                                                                     <  2:63683/149‑1 (MQ=255)                                   tttGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGctt                                                                                                                >  5:232734/1‑149 (MQ=255)                                         gAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGG‑CCCCtt                                                                                                                                       >  6:198995/1‑120 (MQ=255)                                         gAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGG‑CCCCtt                                                                                                                                       <  5:198995/120‑1 (MQ=255)                                               ggcggcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTc                                                                                                      >  7:250103/1‑147 (MQ=255)                                                     cAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGt                                                                                               >  2:171211/1‑148 (MQ=255)                                                       ggAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGAGG‑CCCCTTCCGACCGGTCATACATACAGCTTCTTGTGCTTCCAGCGGTCaa                                                                                            >  2:122875/1‑148 (MQ=255)                                                          gctggtgctggtGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAc                                                                                         <  6:231743/149‑1 (MQ=255)                                                            tggtgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCCTACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAccc                                                                                       >  5:24870/1‑149 (MQ=255)                                                              gtgctggtgcCGTGCCTTAATGATCACCCGGAGTGGGTTTGGGTGTTGGCTTATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCgg                                                                                     <  6:24870/149‑1 (MQ=255)                                                                    gtgccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAAt                                                                               >  1:243281/1‑149 (MQ=255)                                                                     tgccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCGTTGCGG‑CCCCTTGCGACCGGTCATACATCCAGCTTCTTGTGCTTGCCGCTGGCATTACCCGGCAGAAt                                                                               >  1:254915/1‑148 (MQ=255)                                                                       ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg                                                                            <  1:171211/149‑1 (MQ=255)                                                                            ccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCa                                                                       <  8:138351/149‑1 (MQ=255)                                                                             cTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTTATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAg                                                                      <  3:258900/149‑1 (MQ=255)                                                                              tGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAGGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAgg                                                                     >  1:42411/1‑149 (MQ=255)                                                                                  tGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCaa                                                                 >  1:27344/1‑149 (MQ=255)                                                                                       aCCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTc                                                                                              <  6:97573/115‑1 (MQ=255)                                                                                       aCCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTc                                                                                              >  5:97573/1‑115 (MQ=255)                                                                                       aCCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGc                                                             >  4:79116/1‑148 (MQ=255)                                                                                                  gTGAGGGTGTTGGCTTATATGTGTGAAAAGGCCTGAATCGACGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGca                                                 <  1:122875/149‑1 (MQ=255)                                                                                                     aGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGATGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGcaca                                               <  2:42411/148‑1 (MQ=255)                                                                                                         tgtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGc                                           >  6:9257/1‑148 (MQ=255)                                                                                                                tGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATa                                                               <  2:75548/121‑1 (MQ=255)                                                                                                                tGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATa                                                               >  1:75548/1‑121 (MQ=255)                                                                                                                  tcttgtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt                                  <  2:48623/145‑1 (MQ=255)                                                                                                                   tatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCa                                <  3:79116/149‑1 (MQ=255)                                                                                                                     tgtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCa                                                                             >  5:264625/1‑102 (MQ=255)                                                                                                                     tgtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCa                                                                             <  6:264625/102‑1 (MQ=255)                                                                                                                     tgtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAgg                              <  5:9257/149‑1 (MQ=255)                                                                                                                      gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGACCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTAAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt                                  >  1:48623/1‑145 (MQ=255)                                                                                                                      gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt                                  <  6:223076/144‑1 (MQ=255)                                                                                                                      gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt                                  >  5:223076/1‑144 (MQ=255)                                                                                                                      gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGc                             >  5:259072/1‑149 (MQ=255)                                                                                                                         tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGcccc                          <  7:204710/149‑1 (MQ=255)                                                                                                                          gAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACAACCAGCTTGTTGTGCTTGCAGCTGTCCTTACCCGGCAGAATCAGGTTCAGGGCAATAGCCAGAATGGCGCACAGGGCCATGCCCTTCAGGCCCCa                         >  3:12501/1‑149 (MQ=255)                                                                                                                           aaaaGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                         <  6:259072/148‑1 (MQ=255)                                                                                                                                         gAAGGGGCCGCAAAGCGG‑CCCCTTCCGACCGGTCATAAATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGt           >  8:82436/1‑148 (MQ=255)                                                                                                                                           aGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAc        <  8:250103/149‑1 (MQ=255)                                                                                                                                               gCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCt    <  1:58958/149‑1 (MQ=255)                                                                                                                                                ccGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTg   >  6:169353/1‑149 (MQ=255)                                                                                                                                                ccGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTg   <  2:243281/149‑1 (MQ=255)                                                                                                                                                ccGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTg   <  2:254915/149‑1 (MQ=255)                                                                                                                                                 cGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGc  >  4:127640/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGG‑CCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGC  >  NC_002947/863751‑864043

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