Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,031,537	A→T	S421T (TCC→ACC)	manB ←	Phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,031,537	0	A	T	94.4%	46.5 / NA	18	S421T (TCC→ACC)	manB	Phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (8/9);  minor base C (0/1);  total (8/10)
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

TGTCTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAGGAGCGCAGGTTAAAGCGCCAGTCGGCA  >  REL606/2031463‑2031562                                                                           |                          tGTCTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGCATTGgaggag                          >  2:1056907/1‑76 (MQ=255)    cTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTgcgc                       >  2:510068/1‑76 (MQ=255)    cTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTgcgc                       >  2:455603/1‑76 (MQ=255)    cTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTgcgc                       >  2:231553/1‑76 (MQ=255)    cTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTgcgc                       >  1:439277/1‑76 (MQ=255)    cTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTgcgc                       <  2:1105292/76‑1 (MQ=255)      tttCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAg                     >  2:1530243/1‑76 (MQ=255)      tttCTCCCATCAGCGGTACACCCCCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAg                     <  1:950435/76‑1 (MQ=14)       ttcttcCATCATCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg                    >  1:43378/1‑76 (MQ=11)       ttcttcCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg                    >  1:655258/1‑76 (MQ=255)         cttcCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt                  <  1:1570369/76‑1 (MQ=255)         cttcCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt                  >  2:1584203/1‑76 (MQ=255)         cttcCATCAGCGGTACATCACCGCGCGATTCCACACTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt                  >  1:359395/1‑76 (MQ=11)           tcCATCAGCGGTACATCCCCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaa                <  1:1123992/76‑1 (MQ=11)           tcCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaa                <  2:845780/76‑1 (MQ=255)           tcCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaa                <  1:636921/76‑1 (MQ=255)              aTCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCATNTTAAAgc             <  2:1433434/76‑1 (MQ=11)              aTCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc             <  1:348727/76‑1 (MQ=255)               tCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGCGCGCAGGTTAAAgcg            <  1:156940/76‑1 (MQ=255)                   cGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAg        <  1:1567148/76‑1 (MQ=255)                         aTCACCGCGCGATTCCACATTCCACCGGCCCGCCGGTACGGTGTTGTAGGAGCGCATGTTAATGCGTCAGTCGGCa  <  1:1711268/76‑1 (MQ=17)                                                                           |                          TGTCTTTTCTTCCATCAGCGGTACATCACCGCGCGATTCCACATTCAACCGCACCACCGGTTCGGTGTTGGAGGAGCGCAGGTTAAAGCGCCAGTCGGCA  >  REL606/2031463‑2031562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑