Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,600,978	C→T	intergenic (+30/‑223)	bcp → / → hyfA	thiol peroxidase/hydrogenase 4 component A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,600,978	0	C	T	94.4%	46.8 / NA	18	intergenic (+30/‑223)	bcp/hyfA	thiol peroxidase/hydrogenase 4 component A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/10);  minor base G (0/1);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGCTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATG  >  NC_000913/2600732‑2601203                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  gctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGc                                                                                                                                                                                                                               <  2:109515/251‑1 (MQ=255)         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCa                                                                                                                                                                                                                       <  2:143111/251‑1 (MQ=255)                                gACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGa                                                                                                                                                                                                <  1:233582/251‑1 (MQ=255)                                     ccaGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCAACGATGGCATTCAGCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACTCGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGAGCAAGATcaca                                                                                                                                                                                           <  2:87437/251‑1 (MQ=255)                                                   aaTTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTcccc                                                                                                                                                                             <  1:36965/251‑1 (MQ=255)                                                             cTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACCCgg                                                                                                                                                                   >  2:92413/1‑251 (MQ=255)                                                                                         aCCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGACTCCATTTTTTGCTGAACGtt                                                                                                                                       >  1:234210/1‑251 (MQ=255)                                                                                                                               cGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACTCGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTaa                                                                                                 <  2:152241/251‑1 (MQ=255)                                                                                                                                   gATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGCCGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTACATTCCGTGCTGGCTGCGTTTGAGGCAAGCATACCTCACTTATCGTGATCAAGATCACATTCTCGCTTTCCACTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATAc                                                                                              >  1:266750/1‑250 (MQ=255)                                                                                                                                   gATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATa                                                                                               <  1:171368/249‑1 (MQ=255)                                                                                                                                            aTCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACAACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCATTTAATACATAACTTAtt                                                                                    >  2:95989/1‑251 (MQ=255)                                                                                                                                                         ttGAGGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTTCAGTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGGTATCAATATCTCATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACTTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACAGAACTTCATGAATATATTg                                                                          <  2:204404/248‑1 (MQ=255)                                                                                                                                                                          cgCAATCCCCACGACGTTGGGCGGAACTGTCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCTGTTGCGGCCAGCATACCTCACTTCCCGTGATCAAGATCACATTCTCGCCTTCCCCTGCTACACGGGTGTCGAATCCATTTTTTGCTGAACGGTAATGACCATCTTTTTTGTAGCGTTCATAATCCAGTTAATACATAACTTATTGAATATATGGAGGTAATCTGAATGCCATcc                                                      <  2:237689/250‑1 (MQ=255)                                                                                                                                                                                      gacgTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAatat                                          >  1:13868/1‑251 (MQ=255)                                                                                                                                                                                                           aaaGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCa                     <  1:191111/251‑1 (MQ=255)                                                                                                                                                                                                           aaaGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCAAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATc                      >  1:149110/1‑250 (MQ=255)                                                                                                                                                                                                                acacGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTc                >  1:1122/1‑251 (MQ=255)                                                                                                                                                                                                                               ttGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATTCCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTTCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATg  <  1:88417/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  GCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGCTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATG  >  NC_000913/2600732‑2601203

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