Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,289,145	T→C	33.8%	H258H (CAT→CAC)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,289,145	0	T	C	33.8%	47.7 / 36.5	68	H258H (CAT→CAC)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base T (14/31);  new base C (12/11);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGG  >  minE/1289079‑1289213                                                                   |                                                                     gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:2615579/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:1082446/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:2015989/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:1957971/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:1946600/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:2895039/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:1768972/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:1669978/69‑1 (MQ=255) gAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:1615664/69‑1 (MQ=255)  aaCTTGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:2809980/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAATAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:605898/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:2583153/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:2294862/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:2767898/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:2845513/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:1631134/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:475583/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:73929/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAg                                                                   <  1:1060376/69‑1 (MQ=255)  aaCTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCCACCTTGCATGGAAAAAACATCAg                                                                   <  1:1730419/69‑1 (MQ=255)    cTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAGAt                                                                 <  1:2115735/69‑1 (MQ=255)         tCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAg                                                            <  1:3365199/69‑1 (MQ=255)         tCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAg                                                            <  1:107860/69‑1 (MQ=255)         tCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAg                                                            <  1:2748951/69‑1 (MQ=255)                 gAAACCGCAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:1101179/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAATCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:3432404/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:1659121/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:2963479/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:2364862/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:485854/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:2134106/69‑1 (MQ=255)                 gAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATGGCAt                                                    <  1:1025151/69‑1 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACATCa                                                                    <  1:643398/39‑1 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:2322493/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:3228550/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:275138/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:279070/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:380287/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:316270/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:3308152/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:1244341/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:1497298/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:1536276/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:1587393/1‑51 (MQ=255)                               agaagaAGCCATTTCCGACCTTGCATGGAAAAAACACCAGATGCCAGCATg                                                        >  1:198387/1‑51 (MQ=255)                                   gaagCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATg                                    >  1:2667963/1‑67 (MQ=255)                                   gaagCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATg                                    >  1:1055663/1‑67 (MQ=255)                                   gaagCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATg                                    >  1:3172583/1‑67 (MQ=255)                                       ccATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTAccgcc                                        <  1:208990/59‑1 (MQ=255)                                       ccATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTAccgcc                                        <  1:1463574/59‑1 (MQ=255)                                            tCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAg                               <  1:969877/63‑1 (MQ=255)                                                aCCTTGCATGGAAAAAACACCAGATGCCAGCATGGCATTTAATTACCGCTGATGGTCGGGGTATTAct                      <  1:75430/68‑1 (MQ=255)                                                aCCTTGCATGGAAAAAACACCAGATGCCAGCATGGCATTTAATTACCGCTGATGGTCGGGGTATTAct                      <  1:1090459/68‑1 (MQ=255)                                                aCCTTGCATGGAAAAAACACCAGATGCCAGCATGGCATTTAATTACCGCTGATGGTCGGGGTATTAct                      <  1:1355904/68‑1 (MQ=255)                                                     gCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAg                               <  1:2136942/54‑1 (MQ=255)                                                        tGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGaa              <  1:3026426/68‑1 (MQ=255)                                                        tGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGaa              <  1:2777368/68‑1 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:348814/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:2450831/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:268369/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:1441203/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:350129/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:2803741/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:1017254/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:2211737/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:165000/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:3089293/1‑69 (MQ=255)                                                                aCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGt     >  1:1754733/1‑69 (MQ=255)                                                                   tCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTggg  <  1:2168977/69‑1 (MQ=255)                                                                   tCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTggg  <  1:546801/69‑1 (MQ=255)                                                                   tCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTggg  <  1:1169343/69‑1 (MQ=255)                                                                   tCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTggg  <  1:1923152/69‑1 (MQ=255)                                                                   |                                                                     GAACTGGATCACCGCCGAAACCGAAGCGCCAGAAGAAGCCATTTCCGACCTTGCATGGAAAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGG  >  minE/1289079‑1289213

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