Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 294,191 A→G 12.2% F531L (TTC→CTC)  ybaE ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE294,1910AG12.2% 124.1 / 8.9 49F531L (TTC→CTC) ybaEpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base G (3/3);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

ACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTT  >  minE/294122‑294259
                                                                     |                                                                    
aCCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGaa                                                                      >  1:2159661/1‑70 (MQ=255)
        cgcgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGt                                                             <  1:2198781/71‑1 (MQ=255)
          cgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCa                                                           <  1:235105/71‑1 (MQ=255)
             tgtCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCg                                                        >  1:1304642/1‑71 (MQ=255)
             tgtCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCg                                                        <  1:2854586/71‑1 (MQ=255)
             tgtCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCg                                                        >  1:2218926/1‑71 (MQ=255)
                   tCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCAt                                                  <  1:1650048/71‑1 (MQ=255)
                     gCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCt                                                <  1:385674/71‑1 (MQ=255)
                       acacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCa                                                                    >  1:580387/1‑49 (MQ=255)
                        cacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATc                                                 <  1:939808/67‑1 (MQ=255)
                        cacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGt                                              >  1:1791834/1‑70 (MQ=255)
                           cGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCAt                                           >  1:1526667/1‑70 (MQ=255)
                              tCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAGCAg                                                                   <  1:399570/43‑1 (MQ=255)
                              tCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAAc                                       >  1:2905988/1‑71 (MQ=255)
                               cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAg                                                                   <  1:3343819/42‑1 (MQ=255)
                               cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAAc                                       <  1:895775/70‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                               <  1:77934/71‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                               <  1:1316530/71‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                               <  1:37296/71‑1 (MQ=255)
                                      ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                               <  1:502612/71‑1 (MQ=255)
                                            gCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGc                         <  1:3720756/71‑1 (MQ=255)
                                            gCACTAATGCGATAGTGATAGTTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGc                         <  1:2856688/71‑1 (MQ=255)
                                            gCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaaa                             >  1:2713630/1‑67 (MQ=255)
                                            gCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGc                         >  1:117733/1‑71 (MQ=255)
                                               cTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGAGTCGCATCTGTCATTAACTGGCTAAAAAcc                           >  1:959096/1‑66 (MQ=255)
                                                tAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTc                                                      >  1:3680857/1‑38 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaaa                             >  1:2811648/1‑62 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaaa                             >  1:2377613/1‑62 (MQ=255)
                                                 aaTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCttt                      >  1:2769907/1‑69 (MQ=255)
                                                  aTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTc                     <  1:329401/69‑1 (MQ=255)
                                                   tGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                               <  1:988250/58‑1 (MQ=255)
                                                    gCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCa                                                   <  1:2493110/37‑1 (MQ=255)
                                                     cGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaaaa                             >  1:3907849/1‑58 (MQ=255)
                                                     cGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTa                                 <  1:3268206/54‑1 (MQ=255)
                                                      gatagtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCa               <  1:1050592/71‑1 (MQ=255)
                                                        tagtgatagtTGAGCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCa               >  1:1570705/1‑69 (MQ=255)
                                                        tagtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCa                    <  1:2017409/64‑1 (MQ=255)
                                                          gtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAgg                  >  1:3596914/1‑64 (MQ=255)
                                                            gatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAAt         >  1:1638311/1‑71 (MQ=255)
                                                            gatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAAt         >  1:230477/1‑71 (MQ=255)
                                                             atagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAAt         <  1:1025216/70‑1 (MQ=255)
                                                             atagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATc        <  1:612628/71‑1 (MQ=255)
                                                             atagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATc        <  1:148833/71‑1 (MQ=255)
                                                              tagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCg       <  1:2576527/71‑1 (MQ=255)
                                                               agtTGAACATCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCgg      >  1:2371393/1‑71 (MQ=255)
                                                                gtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGt     <  1:663764/71‑1 (MQ=255)
                                                                 ttGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaaa                               <  1:846795/44‑1 (MQ=255)
                                                                 ttGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaa                                <  1:3859708/43‑1 (MQ=255)
                                                                 ttGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTaa                                <  1:3661643/43‑1 (MQ=255)
                                                                   gAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGc                                   >  1:2015482/1‑38 (MQ=255)
                                                                    aaCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAgg                  <  1:2475354/54‑1 (MQ=255)
                                                                    aaCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCgg      <  1:1734063/66‑1 (MQ=255)
                                                                    aaCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGtttt  <  1:2240995/70‑1 (MQ=255)
                                                                    aaCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGtttt  <  1:3586270/70‑1 (MQ=255)
                                                                     |                                                                    
ACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTT  >  minE/294122‑294259

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