Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,677,313	T→C	42.8%	D206G (GAC→GGC)	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,313	0	T	C	42.8%	18.8 / 58.2	77	D206G (GAC→GGC)	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base T (9/35);  new base C (13/20);  total (22/55)
Fisher's exact test for biased strand distribution p-value = 8.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTT  >  minE/2677254‑2677376                                                            |                                                                tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:1035465/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:35410/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:364094/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:3056273/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:2969661/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:272739/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:846985/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:1164931/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:2353152/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCg                                                         <  1:975116/68‑1 (MQ=255) tAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCACTAGCg                                                         <  1:2755582/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTATCGTATc                                                     <  1:435918/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:2240195/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:798744/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:321868/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:143213/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:372919/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:1566201/68‑1 (MQ=255)     aCGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:3257266/68‑1 (MQ=255)       gTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCac                                                  <  1:3292178/69‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2094274/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2340988/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:1107969/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2394190/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2405307/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:800955/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:22394/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2574127/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:102402/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:1278301/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2736733/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:128064/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2862455/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:709376/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:536461/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:3140427/67‑1 (MQ=255)                   cggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:1424381/67‑1 (MQ=255)                   cAGTTTGCAGTATGTAATTAAACAGGCCTTTGATTTGTGTGTCATTAGCg                                                         <  1:3163215/50‑1 (MQ=255)                    ggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:2597244/67‑1 (MQ=255)                    ggTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:146350/67‑1 (MQ=255)                    aGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCA‑TAGCGTATCCACTACTTAATAAGCAg                                    <  1:638106/69‑1 (MQ=255)                     gTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAg                                    <  1:977956/69‑1 (MQ=255)                     gTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAg                                    <  1:1992846/69‑1 (MQ=255)                     gTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAg                                    <  1:969978/69‑1 (MQ=255)                     gTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAg                                    <  1:3361474/69‑1 (MQ=255)                     gTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAg                                    <  1:324404/69‑1 (MQ=255)                      tttGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAg                                    <  1:38090/68‑1 (MQ=255)                              aTGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGc                           <  1:3367738/69‑1 (MQ=255)                              aTGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGc                           <  1:123007/69‑1 (MQ=255)                               tGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGc                           >  1:2856707/1‑68 (MQ=255)                                 tAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     <  1:1517202/40‑1 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:1808136/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:647630/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:1792161/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:1735880/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:1717779/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:3363385/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:1372905/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:483444/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:1866685/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:3024914/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        >  1:2269440/1‑67 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctc                         >  1:2225529/1‑66 (MQ=255)                                   aTTAAACAGCCCTTTGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGCAGGACATTTGctct                        >  1:2746436/1‑67 (MQ=255)                                     tAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     >  1:3107281/1‑36 (MQ=255)                                     tAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     >  1:753988/1‑36 (MQ=255)                                     tAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     >  1:2124622/1‑36 (MQ=255)                                     tAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATc                                                     >  1:1311605/1‑36 (MQ=255)                                                ttGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAgg         >  1:1051168/1‑69 (MQ=255)                                                ttGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAgg         >  1:2362666/1‑69 (MQ=255)                                                ttGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAgg         >  1:1175562/1‑69 (MQ=255)                                                ttGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCACATCAgg         >  1:1593488/1‑69 (MQ=255)                                                 tGATTTGTCTGTCATTAGCGTATCCACTTCTTAATAAGCAGGCCATTTGctct                        <  1:337832/53‑1 (MQ=255)                                                 tGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        <  1:2776199/53‑1 (MQ=255)                                                 tGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        <  1:2819898/53‑1 (MQ=255)                                                 tGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGctct                        <  1:1928656/53‑1 (MQ=255)                                                 tGATTTGTCTGCCATTAGCGTATCCACTACTTAATAAGc                                      <  1:907019/39‑1 (MQ=255)                                                      tgtctgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTcgcg    >  1:3379366/1‑68 (MQ=255)                                                      tgtctgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGt   >  1:2096951/1‑69 (MQ=255)                                                      tgtctgccTTTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2323167/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTGCCAGCCCCAt             <  1:2368670/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:898624/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2278526/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:699303/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:1350958/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:497430/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:3119984/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2828176/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2048977/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2057799/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2493956/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:2264300/59‑1 (MQ=255)                                                      tgtctgccATTAGCGTATCCACTACTTAATAAGCAGGCAATTTGCTCTACCAGCCCCAt             <  1:1306480/59‑1 (MQ=255)                                                       gtctgccATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCAt             <  1:1705740/58‑1 (MQ=255)                                                        tctgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGtt  <  1:690958/68‑1 (MQ=255)                                                            |                                                                TAAAACGTACAGCGTCGCCAGTTTGCAGTATGTAATTAAACAGCCCTTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTT  >  minE/2677254‑2677376

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: ‑