Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,692,670 C→G 12.5% V102V (GTG→GTC iscU ← scaffold protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,692,6700CG12.5% 131.1 / 10.9 48V102V (GTG→GTCiscUscaffold protein
Reads supporting (aligned to +/- strand):  ref base C (15/27);  new base G (2/4);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

ACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCGAGAGACTTC  >  minE/1692601‑1692739
                                                                     |                                                                     
aCGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCa                                                                      >  1:1347495/1‑71 (MQ=255)
   tttGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcg                                                                   <  1:2410453/71‑1 (MQ=255)
   tttGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcg                                                                   >  1:1702770/1‑71 (MQ=255)
     tGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcggc                                                                 >  1:3355487/1‑71 (MQ=255)
                  cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGtt                                                    >  1:1497342/1‑71 (MQ=255)
                  cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGtt                                                    >  1:3424777/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)
                      aTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCa                                                        <  1:276956/63‑1 (MQ=255)
                        ggcggcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAg                                              >  1:3100382/1‑71 (MQ=255)
                           ggcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                                           <  1:2096005/71‑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:1043055/45‑1 (MQ=255)
                               ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                                 <  1:1181764/45‑1 (MQ=255)
                               ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                                 <  1:1873922/45‑1 (MQ=255)
                               ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                                 <  1:2998770/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:1956211/1‑54 (MQ=255)
                                            ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                                           >  1:3256441/1‑54 (MQ=255)
                                            ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg                             <  1:2089608/68‑1 (MQ=255)
                                            ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg                             <  1:1436081/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:1315728/65‑1 (MQ=255)
                                                  gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc                          <  1:2593588/65‑1 (MQ=255)
                                                  gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc                          <  1:1075108/65‑1 (MQ=255)
                                                    aTAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTgcgc                  >  1:1221539/1‑71 (MQ=255)
                                                        aaCAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTgc                    <  1:1600955/65‑1 (MQ=255)
                                                          cAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTtcgtc            <  1:3302325/71‑1 (MQ=255)
                                                          cAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTtcgtc            >  1:1116636/1‑71 (MQ=255)
                                                            gTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTtcgtc            >  1:2885693/1‑69 (MQ=255)
                                                               aaTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGcc                       >  1:280225/1‑55 (MQ=255)
                                                               aaTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCgaga        <  1:2815253/70‑1 (MQ=255)
                                                                aTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                                          <  1:1034173/35‑1 (MQ=255)
                                                                aTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCgagag       <  1:2215765/70‑1 (MQ=255)
                                                                     cACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCGAGAGACTTc  <  1:1095744/70‑1 (MQ=255)
                                                                     cACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCGAGAGACTTc  <  1:2423787/70‑1 (MQ=255)
                                                                     cACCGGCGGCAGTTCAAGTTCTTCAGCAATATAGGTGTTTTTg                             >  1:2470378/1‑43 (MQ=255)
                                                                     |                                                                     
ACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCGAGAGACTTC  >  minE/1692601‑1692739

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: