Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,928 G→C 14.4% A66A (GCG→GCC yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,9280GC14.4% 227.3 / 26.1 90A66A (GCG→GCCyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base C (8/5);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

CATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCG  >  minE/2569863‑2569991
                                                                 |                                                               
cATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACg                                                              >  1:4651877/1‑69 (MQ=255)
cATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGc                                                             >  1:876759/1‑70 (MQ=255)
 aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                           <  1:17661/71‑1 (MQ=255)
          gttGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                         >  1:5410307/3‑64 (MQ=255)
          gttGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                         >  1:4930081/3‑64 (MQ=255)
          gCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:1167625/1‑57 (MQ=255)
          gCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCg                                                  <  1:3910979/71‑1 (MQ=255)
          gCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGAGCGACg                                                              <  1:734327/59‑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:4988941/70‑1 (MQ=255)
                  aGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTg                                           <  1:179309/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:2954639/71‑1 (MQ=255)
                        tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                    <  1:3236960/71‑1 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:4199155/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:3554903/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:3539379/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:1533267/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:1526812/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:1318107/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                >  1:1296281/1‑42 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                              >  1:3863784/1‑60 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                              >  1:692925/1‑60 (MQ=255)
                         ttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                              >  1:4911116/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:34380/49‑1 (MQ=16)
                                   tGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAg                                               <  1:318272/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:4453950/66‑1 (MQ=255)
                                           ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                      <  1:148293/66‑1 (MQ=255)
                                           ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                      <  1:638176/66‑1 (MQ=255)
                                           ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                      <  1:614750/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:55658/41‑1 (MQ=255)
                                                         aTCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTAt                                 <  1:4576953/41‑1 (MQ=255)
                                                           cTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGt                >  1:4983688/1‑56 (MQ=255)
                                                              ggCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGc   <  1:789601/66‑1 (MQ=255)
                                                              ggCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGc   <  1:1425475/66‑1 (MQ=255)
                                                              ggCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCg  <  1:2495885/67‑1 (MQ=255)
                                                                 |                                                               
CATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCG  >  minE/2569863‑2569991

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: