Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,032,758	G→A	L14L (CTA→TTA)	manB ←	Phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,032,758	0	G	A	92.3%	42.4 / ‑0.5	26	L14L (CTA→TTA)	manB	Phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (12/12);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GTTTTCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAGTTTCCCGCGAATATCATAAGCTTTAA  >  REL606/2032685‑2032784                                                                          |                           gTTTTCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAAtt                          <  2:724128/76‑1 (MQ=255)  ttttCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAAttt                         <  1:1901586/76‑1 (MQ=255)  ttttCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAAttt                         >  1:222342/1‑76 (MQ=255)  ttttCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAAttt                         >  2:191996/1‑76 (MQ=255)   tttCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAAtttt                        <  2:1032011/76‑2 (MQ=255)    ttCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAAttttc                       >  2:635991/1‑74 (MQ=255)     tCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTGAGTTCTTCGCCTAATTTtcc                      <  2:1598280/76‑1 (MQ=14)      cGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATGTTCATTCAGTTCTTCGCCTAATTTtcca                     >  1:1113302/1‑75 (MQ=11)      aggTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAGTTTCCcg                     >  1:1219216/2‑76 (MQ=255)        gTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCcacg                   <  1:1991951/76‑1 (MQ=11)        gTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCcacg                   >  2:1512598/1‑76 (MQ=11)        gTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCcacg                   >  2:1072324/1‑76 (MQ=11)          ttGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGaa                 >  1:1953900/1‑76 (MQ=11)          ttGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGaa                 >  1:831719/1‑76 (MQ=11)           taagaAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAat                <  2:209221/74‑1 (MQ=11)           tGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAat                >  2:885915/1‑76 (MQ=11)           tGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAat                >  1:704963/1‑76 (MQ=11)           tGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAat                >  1:1240070/1‑76 (MQ=11)           tGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAat                <  1:1195033/76‑1 (MQ=11)             agaAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAatat              <  2:1985549/76‑1 (MQ=11)              gaAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAATATc             >  1:1151080/1‑76 (MQ=11)               aaaTTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAATATCa            <  2:397413/76‑1 (MQ=11)               aaaTTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAATATCa            <  2:397430/76‑1 (MQ=11)                 aTTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAATATCATa          <  1:247266/76‑1 (MQ=11)                  ttCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTTCCACGAATATCATag         <  1:508252/76‑2 (MQ=11)                         tAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAGTTTCCCGCGAATATCATAAGCTTTaa  >  1:306409/1‑76 (MQ=255)                                                                          |                           GTTTTCGGTTTGAGAAATTCGCCATAGGCACGCCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAGTTTCCCGCGAATATCATAAGCTTTAA  >  REL606/2032685‑2032784

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