Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,683,278 T→G 13.2% Y145S (TAT→TCT)  yhgF ← predicted transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,683,2780TG13.2% 66.6 / 12.7 38Y145S (TAT→TCT) yhgFpredicted transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base T (17/16);  new base G (2/3);  total (19/19)
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.80e-01

ATCAGGATATAGCGCGCGCCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGCCAACGGCTC  >  minE/2683212‑2683342
                                                                  |                                                                
aTCAGGATATAGCGCGCGCCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACata                                                                 >  1:1838651/1‑68 (MQ=255)
 tCAGGATATAGCGCGCGCCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACat                                                                  <  1:682619/66‑1 (MQ=255)
      atatAGCGCGCGCCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAgc                                                        <  1:1992037/71‑1 (MQ=255)
         taGCGCGCGCCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGCGCAgcggc                                                     <  1:1915487/71‑1 (MQ=255)
            cgcgcgcCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGAc                                                  <  1:77253/71‑1 (MQ=255)
                       cAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCgtgt                                       >  1:299290/1‑71 (MQ=255)
                       cAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCgtg                                        <  1:118044/70‑1 (MQ=255)
                       cAGCGCGGCTTTGGTATCTG‑‑ACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCgtgtgt                                     >  1:1589796/1‑71 (MQ=255)
                          cgcgGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCgtgtgt                                     <  1:556031/70‑1 (MQ=255)
                             ggCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACg                                 <  1:1239745/71‑1 (MQ=255)
                              gCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATtgtga                                                           <  1:246588/44‑2 (MQ=255)
                              gCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACgg                                <  1:574346/71‑1 (MQ=255)
                                  tGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATc                             >  1:1289139/1‑70 (MQ=255)
                                     tATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACg                                 <  1:1461271/63‑1 (MQ=255)
                                     tATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACg                                 <  1:1373054/63‑1 (MQ=255)
                                       tCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACgg                                <  1:2178940/62‑1 (MQ=255)
                                       tCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCca                       >  1:1698869/1‑71 (MQ=255)
                                        cTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTcc                        >  1:466551/1‑69 (MQ=255)
                                         tGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGAt                              >  1:623388/1‑62 (MQ=255)
                                         tGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCca                       >  1:1668230/1‑69 (MQ=255)
                                         tGCCACGCCTTTATCGGCATAAACAGATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCca                       >  1:1282266/1‑69 (MQ=255)
                                          gCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCAcag                    >  1:531462/1‑71 (MQ=255)
                                            cACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCAcagca                  <  1:1184183/71‑1 (MQ=255)
                                              cGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAgg                <  1:514106/71‑1 (MQ=255)
                                                ccTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGt               >  1:1840297/1‑70 (MQ=255)
                                                ccTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGt               >  1:1940773/1‑70 (MQ=255)
                                                  tttATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCgtgtgt                                     <  1:993030/46‑1 (MQ=255)
                                                  tttATCGGCATAAACAGATTGTGCAGCGGCGACTTCTGGCGTGTGTGACgg                                <  1:1613880/51‑1 (MQ=255)
                                                  tttATCGGCATAAACAGATTGTGCAGCGGCGACTTCTGGCGTGTGTGACgg                                <  1:1857307/51‑1 (MQ=255)
                                                  tttATCGGCATAAACAGATTGTGCAGCGGCGACTTCTGGCGTGTGTGACgg                                <  1:1753776/51‑1 (MQ=255)
                                                    tATCGGCATAAACATATTGTGCAGCTGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGc           >  1:1491683/1‑70 (MQ=255)
                                                    tATCGGCATAAACAGATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGc           >  1:713127/1‑70 (MQ=255)
                                                     aTCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGcc          >  1:593307/1‑70 (MQ=255)
                                                        ggCATAAACATATTGTGCAGCGGCCACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGCTCAGCCAACg      >  1:659635/1‑71 (MQ=255)
                                                            tAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGGCAGCCAACGGCTc  <  1:265892/71‑1 (MQ=255)
                                                             aaaCATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGCCAACGGCTc  <  1:1198074/70‑1 (MQ=255)
                                                              aaCATATTGTGCAGCGGCGACTTCTGGCGTGTGTGa                                   >  1:771346/1‑36 (MQ=255)
                                                              aaCATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCca                       >  1:1806122/1‑48 (MQ=255)
                                                              aaCATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCca                       >  1:1004702/1‑48 (MQ=255)
                                                              aaCATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGCCAACGGCTc  >  1:1588521/1‑69 (MQ=255)
                                                                  tatTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGCCAACGGCTc  <  1:452108/65‑1 (MQ=255)
                                                                  tatTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGCCAACGGCTc  <  1:1747720/65‑1 (MQ=255)
                                                                  |                                                                
ATCAGGATATAGCGCGCGCCATCCAGCGCGGCTTTGGTATCTGCCACGCCTTTATCGGCATAAACATATTGTGCAGCGGCGACTTCTGGCGTGTGTGACGGATCGCTCCACAGCAGGTCAGCCAACGGCTC  >  minE/2683212‑2683342

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: