Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,937 T→C 12.3% Y69Y (TAT→TAC yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,9370TC12.3% 226.0 / 18.8 89Y69Y (TAT→TACyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/48);  new base C (6/5);  total (36/53)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGAAAT  >  minE/2569873‑2570005
                                                                |                                                                    
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: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:2954639/71‑1 (MQ=255)
              tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                                  <  1:3236960/71‑1 (MQ=255)
              tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                                  <  1:688267/71‑1 (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:976676/1‑53 (MQ=255)
                tACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAc                                                                  >  1:4468382/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)
                      acATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCa                                            >  1:1599532/2‑69 (MQ=255)
                      acATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCa                                            >  1:2574627/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:614750/66‑1 (MQ=255)
                                 ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                                    <  1:1051868/66‑1 (MQ=255)
                                 ggCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                                    <  1:148293/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)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGc                                <  1:4826910/41‑1 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGc                                <  1:3639109/41‑1 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGc                 <  1:2164856/56‑1 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGa               >  1:532676/1‑58 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGTaa    <  1:2769149/69‑1 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGaa    <  1:2346235/69‑1 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGaa    <  1:715035/69‑1 (MQ=255)
                                                              tATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGaa    <  1:1000644/69‑1 (MQ=255)
                                                               aTCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaaa                   >  1:2839480/1‑53 (MQ=255)
                                                               aTCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaaa                   >  1:3161984/1‑53 (MQ=255)
                                                               aTCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGc                 <  1:3329545/55‑1 (MQ=255)
                                                               aTCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGAAAt  >  1:5274839/1‑70 (MQ=255)
                                                               aTCAACAGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaaa                   >  1:10238/1‑53 (MQ=255)
                                                                tCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaa                    <  1:2675112/51‑1 (MQ=255)
                                                                tCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGaaa   <  1:972953/68‑1 (MQ=255)
                                                                tCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGAAAt  <  1:1665210/69‑1 (MQ=255)
                                                                |                                                                    
GCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGACTGACCCGGAAAT  >  minE/2569873‑2570005

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: