Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,934 C→T 16.4% L68L (CTC→CTT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,9340CT16.4% 184.2 / 24.3 79L68L (CTC→CTTyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/41);  new base T (8/5);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCC  >  minE/2569864‑2569999
                                                                      |                                                                 
aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                                   <  1:17661/71‑1 (MQ=255)
         gttGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                                 >  1:4930081/3‑64 (MQ=255)
         gttGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                                 >  1:5410307/3‑64 (MQ=255)
         gCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCg                                                          <  1:3910979/71‑1 (MQ=255)
            gggTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAg                                                       >  1:2167260/1‑71 (MQ=255)
            gggTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAg                                                       >  1:5049589/1‑71 (MQ=255)
             ggTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAg                                                       <  1:408269/70‑1 (MQ=255)
                gAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCggg                                                   >  1:2862595/1‑69 (MQ=255)
                gAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTg                                                   >  1:2055722/1‑71 (MQ=255)
                 aGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTgg                                                  >  1:1817937/1‑71 (MQ=255)
                 aGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTg                                                   <  1:179309/70‑1 (MQ=255)
                 aGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTg                                                   <  1:4988941/70‑1 (MQ=255)
                     tGTTTAGTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgg                                              >  1:738884/1‑71 (MQ=255)
                     tGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTg                                               <  1:264557/70‑1 (MQ=255)
                      gTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGt                                                    <  1:3236798/64‑1 (MQ=255)
                      gTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                             >  1:4756302/1‑71 (MQ=255)
                       tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                            <  1:688267/71‑1 (MQ=255)
                       tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                            <  1:3236960/71‑1 (MQ=255)
                       tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                            <  1:2954639/71‑1 (MQ=255)
                        ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                                      >  1:692925/1‑60 (MQ=255)
                        ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                                      >  1:4911116/1‑60 (MQ=255)
                        ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                                      >  1:3863784/1‑60 (MQ=255)
                         tACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAc                                                            >  1:4468382/1‑53 (MQ=255)
                         tACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAc                                                            >  1:976676/1‑53 (MQ=255)
                         tACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                             >  1:2295864/1‑68 (MQ=255)
                         tACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTa                                          >  1:4382838/1‑71 (MQ=255)
                           ctctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCgg                                                         <  1:4119199/54‑1 (MQ=255)
                             ttacATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCa                                      <  1:985011/68‑1 (MQ=255)
                              tacATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCgggg                                                  <  1:600294/56‑1 (MQ=255)
                              tGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTa                                                                 <  1:4915046/43‑1 (MQ=255)
                               acATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCa                                      >  1:2574627/2‑69 (MQ=255)
                               acATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCa                                      >  1:1599532/2‑69 (MQ=255)
                                cATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGt                                    <  1:2119049/70‑1 (MQ=255)
                                  tGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAg                                                       <  1:318272/49‑1 (MQ=16)
                                  tGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAg                                                       <  1:34380/49‑1 (MQ=16)
                                  tGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAcc                                                           <  1:4021559/45‑1 (MQ=255)
                                  tGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTcac                                 <  1:3253850/71‑1 (MQ=255)
                                  tGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTcac                                 <  1:1913031/71‑1 (MQ=255)
                                   gTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATgg                                              >  1:2496781/1‑57 (MQ=255)
                                   gTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTca                                  >  1:4231000/1‑69 (MQ=255)
                                    ttATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGa                                                        >  1:999846/1‑46 (MQ=16)
                                    ttATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTa                                          >  1:2835761/1‑60 (MQ=255)
                                    ttATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACAt                               <  1:1248647/71‑1 (MQ=255)
                                     tATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                             >  1:4415900/1‑56 (MQ=255)
                                     tATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                              >  1:5597047/1‑71 (MQ=255)
                                      aTCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCgggga                                                 >  1:2988230/1‑50 (MQ=255)
                                      aTCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGc                                       <  1:468464/61‑1 (MQ=255)
                                          ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                              <  1:1051868/66‑1 (MQ=255)
                                          ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                              <  1:148293/66‑1 (MQ=255)
                                          ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                              <  1:614750/66‑1 (MQ=255)
                                          ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                              <  1:4453950/66‑1 (MQ=255)
                                          ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                              <  1:638176/66‑1 (MQ=255)
                                          ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCg                         >  1:3336548/1‑71 (MQ=255)
                                           gCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGtt                            <  1:4769139/67‑1 (MQ=255)
                                           gCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGt                        <  1:1753047/71‑1 (MQ=255)
                                             cGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCg                      <  1:4391573/71‑1 (MQ=255)
                                             cGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCg                      <  1:3495332/71‑1 (MQ=255)
                                              gTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGc                     <  1:505330/71‑1 (MQ=255)
                                               tttCTGTGGATCTGGGCGACGCTCTATCAACCGGa                                                        <  1:1133980/35‑1 (MQ=255)
                                               tttCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTg                                               <  1:718837/44‑1 (MQ=255)
                                               tttCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTAt                                         <  1:2822988/50‑1 (MQ=255)
                                                ttCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                            <  1:3754158/46‑1 (MQ=255)
                                                 tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGt                                                    >  1:5377266/1‑37 (MQ=255)
                                                 tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTg                                                   <  1:1495127/38‑1 (MQ=255)
                                                 tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGc                                                 >  1:5674379/1‑40 (MQ=255)
                                                 tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTcac                                 >  1:3816473/1‑56 (MQ=255)
                                                 tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAg                   >  1:578893/1‑70 (MQ=255)
                                                    gtgGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGgag                >  1:5551162/1‑70 (MQ=255)
                                                     tgGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGga                 <  1:3482460/68‑1 (MQ=255)
                                                     tgGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaa              <  1:3762746/71‑1 (MQ=255)
                                                      ggATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGt                                           <  1:5142237/41‑1 (MQ=255)
                                                      ggATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaa              <  1:3769422/70‑1 (MQ=255)
                                                       gATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAg            >  1:145935/1‑71 (MQ=255)
                                                        aTCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTAt                                         <  1:4576953/41‑1 (MQ=255)
                                                        aTCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTAt                                         <  1:55658/41‑1 (MQ=255)
                                                          cTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGt                        >  1:4983688/1‑56 (MQ=255)
                                                             ggCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGc           <  1:1425475/66‑1 (MQ=255)
                                                             ggCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGc           <  1:789601/66‑1 (MQ=255)
                                                             ggCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCg          <  1:2495885/67‑1 (MQ=255)
                                                                 aCGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGAccc  <  1:143264/71‑1 (MQ=255)
                                                                      |                                                                 
ATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCC  >  minE/2569864‑2569999

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: