Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,692,655 A→C 10.3% S107S (TCT→TCG iscU ← scaffold protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,692,6550AC10.3% 146.2 / 13.6 58S107S (TCT→TCGiscUscaffold protein
Reads supporting (aligned to +/- strand):  ref base A (20/32);  new base C (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

ATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGC  >  minE/1692591‑1692723
                                                                |                                                                    
aTTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAg                                                                >  1:3029861/1‑71 (MQ=255)
aTTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAg                                                                >  1:24903/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:716733/1‑42 (MQ=255)
                          tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                                   >  1:1823433/1‑42 (MQ=255)
                          tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                                   >  1:1309593/1‑42 (MQ=255)
                          tCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGaa                                                                   >  1:204809/1‑42 (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)
                                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:1181764/45‑1 (MQ=255)
                                         ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                 <  1:2998770/45‑1 (MQ=255)
                                         ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                 <  1:1043055/45‑1 (MQ=255)
                                         ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                 <  1:1873922/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:1315728/65‑1 (MQ=255)
                                                            gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc          <  1:1075108/65‑1 (MQ=255)
                                                            gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc          <  1:2593588/65‑1 (MQ=255)
                                                              aTAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTgcgc  >  1:1221539/1‑71 (MQ=255)
                                                                |                                                                    
ATTTTGCTTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGC  >  minE/1692591‑1692723

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: