Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,692,652	A→G	10.3%	I108I (ATT→ATC)	iscU ←	scaffold protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,692,652	0	A	G	10.3%	169.8 / 9.7	58	I108I (ATT→ATC)	iscU	scaffold protein
Reads supporting (aligned to +/- strand):  ref base A (20/32);  new base G (2/4);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CAACTCTTATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATC  >  minE/1692583‑1692715                                                                      |                                                                cAACTCTTATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGaa                                                                >  1:3335763/1‑71 (MQ=255) cAACTCTTATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGaa                                                                >  1:2487130/1‑71 (MQ=255)  aaCTCTTATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGaa                                                                <  1:848520/70‑1 (MQ=255)      ctTATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAAt                                                               >  1:985839/1‑67 (MQ=255)         aTTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAg                                                        >  1:24903/1‑71 (MQ=255)         aTTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAg                                                        >  1:3029861/1‑71 (MQ=255)              gCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGaa                                                    <  1:1554593/70‑1 (MQ=255)               cTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACa                                                         <  1:523401/64‑1 (MQ=255)               cTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACa                                                         <  1:2412500/64‑1 (MQ=255)                ttCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                 >  1:199710/1‑71 (MQ=255)                  cACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                 <  1:955476/69‑1 (MQ=255)                  cACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                 <  1:6758/69‑1 (MQ=255)                  cACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                 <  1:2710043/69‑1 (MQ=255)                   aCGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCa                                              >  1:1347495/1‑71 (MQ=255)                      tttGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcg                                           >  1:1702770/1‑71 (MQ=255)                      tttGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcg                                           <  1:2410453/71‑1 (MQ=255)                        tGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcggc                                         >  1:3355487/1‑71 (MQ=255)                                 aGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTg                                                      <  1:2637660/49‑1 (MQ=255)                                   tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                           >  1:1823433/1‑42 (MQ=255)                                   tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                           >  1:1309593/1‑42 (MQ=255)                                   tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                           >  1:716733/1‑42 (MQ=255)                                   tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                           >  1:204809/1‑42 (MQ=255)                                     cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGtt                            >  1:3424777/1‑71 (MQ=255)                                     cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGtt                            >  1:1497342/1‑71 (MQ=255)                                     cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGt                             >  1:2907000/1‑70 (MQ=255)                                        aaTGGCGGCTTTGATCGCGTCTTCTGCCAGATTAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGttctt                         <  1:3380813/71‑1 (MQ=255)                                        aaTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGttctt                         >  1:1908823/1‑71 (MQ=255)                                         aTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAtttt                                                >  1:1544570/1‑47 (MQ=255)                                         aTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCa                                <  1:276956/63‑1 (MQ=255)                                          tGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                           <  1:2225560/35‑1 (MQ=255)                                          tGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                           <  1:1981750/35‑1 (MQ=255)                                           ggcggcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAg                      >  1:3100382/1‑71 (MQ=255)                                              ggcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                   <  1:2096005/71‑1 (MQ=255)                                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                 <  1:1166503/40‑1 (MQ=255)                                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTGAAGTTCTTCAGCaa                   <  1:2495944/70‑1 (MQ=255)                                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                   <  1:1854766/70‑1 (MQ=255)                                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                  <  1:2015791/71‑1 (MQ=255)                                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                  <  1:2523917/71‑1 (MQ=255)                                                cTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                  <  1:3324123/70‑1 (MQ=255)                                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                         <  1:2998770/45‑1 (MQ=255)                                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                         <  1:1181764/45‑1 (MQ=255)                                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                         <  1:1873922/45‑1 (MQ=255)                                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                         <  1:1043055/45‑1 (MQ=255)                                                    gATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAata                 <  1:1494418/67‑1 (MQ=255)                                                    gATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggcgg                                       >  1:1606743/1‑45 (MQ=255)                                                        gcgTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGtt         <  1:846428/71‑1 (MQ=255)                                                         cgTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGttt        <  1:3058147/71‑1 (MQ=255)                                                             ttctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGc                     <  1:366366/54‑1 (MQ=255)                                                             ttctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg     <  1:2511449/70‑1 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAg                              >  1:2668090/1‑43 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAg                      <  1:836346/51‑1 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                   >  1:3256441/1‑54 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                   >  1:1956211/1‑54 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg     <  1:1436081/68‑1 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg     <  1:2089608/68‑1 (MQ=255)                                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGa    <  1:1455435/69‑1 (MQ=255)                                                               cggCCAGGATCGAGCAGTGAATCTTGACCGGCGGCAGGGCCAGttc                           >  1:2866152/3‑46 (MQ=255)                                                                  ccAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc  <  1:2870129/68‑1 (MQ=255)                                                                     gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc  <  1:2593588/65‑1 (MQ=255)                                                                     gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc  <  1:1315728/65‑1 (MQ=255)                                                                     gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc  <  1:1075108/65‑1 (MQ=255)                                                                      |                                                                CAACTCTTATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATC  >  minE/1692583‑1692715

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