Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,966 C→G 9.7% A392A (GCC→GCG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,9660CG 9.7% 167.2 / 7.7 62A392A (GCC→GCGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (4/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

AATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGT  >  minE/385897‑386034
                                                                     |                                                                    
aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGcc                                                                      >  1:406668/1‑70 (MQ=255)
 aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                                                                    >  1:2331529/1‑71 (MQ=255)
  tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                                   <  1:3649642/71‑1 (MQ=255)
  tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                                   <  1:1408146/71‑1 (MQ=255)
      gtCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACa                                                               <  1:2450192/70‑1 (MQ=255)
      aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                                   >  1:1329927/1‑67 (MQ=255)
      aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGa                                                                  >  1:1494920/1‑68 (MQ=255)
      aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACa                                                               <  1:2018197/71‑1 (MQ=255)
         cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGcc                                                            <  1:1425091/71‑1 (MQ=255)
         cTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACac                                                              >  1:2351686/1‑68 (MQ=255)
           tCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                           >  1:1677542/1‑70 (MQ=255)
            cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                           <  1:3012109/69‑1 (MQ=255)
            cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                                         >  1:3528104/1‑70 (MQ=255)
             gTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                                   <  1:3506746/60‑1 (MQ=255)
                    ttATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                                  >  1:2256305/1‑70 (MQ=255)
                    ttATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                                  >  1:3021694/1‑70 (MQ=255)
                    ttATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTg                                                 <  1:14396/71‑1 (MQ=255)
                    ttATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGaa                                                                 >  1:1355696/1‑55 (MQ=255)
                    ttATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGaa                                                                 >  1:3722328/1‑55 (MQ=255)
                     tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAAc                                                                <  1:2112680/55‑1 (MQ=255)
                      aTGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGa                                                                  <  1:1431758/52‑1 (MQ=255)
                       tGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                                  <  1:2764697/67‑1 (MQ=255)
                          cgcTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAGCAGCCCGATGGCGGTGCTGGCt                                           <  1:724819/71‑1 (MQ=255)
                          cgcTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCt                                           <  1:449226/71‑1 (MQ=255)
                           gcTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTc                                          >  1:1405036/1‑71 (MQ=255)
                            cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                                   <  1:110534/45‑1 (MQ=255)
                            cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAAGAACAGcc                                                            >  1:2128327/1‑52 (MQ=255)
                               gTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCgg                                       <  1:3141410/70‑1 (MQ=255)
                                  tCCGCCTGTATCATCGATCCGGGCGTGTATTTTGcc                                                                      >  1:1259441/1‑36 (MQ=255)
                                   ccGCCTGTATCATCGATCCGGGCGTGTATTTTGcc                                                                      >  1:375767/1‑35 (MQ=255)
                                   ccGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAgg                                   >  1:2988646/1‑70 (MQ=255)
                                    cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAggg                                  <  1:1597973/70‑1 (MQ=255)
                                    cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGa                                 >  1:254284/1‑71 (MQ=255)
                                    cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGa                                 >  1:3278000/1‑71 (MQ=255)
                                      ccTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcg                               <  1:285236/71‑1 (MQ=255)
                                         gTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcggcg                            >  1:3698815/1‑71 (MQ=255)
                                          tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGc                                            <  1:2039867/54‑1 (MQ=255)
                                          tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGc                                            <  1:1562283/54‑1 (MQ=255)
                                          tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcggcg                            <  1:2980094/70‑1 (MQ=255)
                                             catcGATCCGTGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGAtg                        >  1:2472677/1‑71 (MQ=255)
                                              atcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCgg                                       >  1:2122624/1‑55 (MQ=255)
                                                cGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTgg                     >  1:3672266/1‑71 (MQ=255)
                                                cGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTgg                     >  1:1207560/1‑71 (MQ=255)
                                                 gATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGc                                      >  1:2911890/1‑53 (MQ=255)
                                                 gATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCa                                     >  1:547902/1‑54 (MQ=255)
                                                 gATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTgg                     <  1:3684144/70‑1 (MQ=255)
                                                   tCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGt                    <  1:3016322/69‑1 (MQ=255)
                                                     cGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCt                <  1:3695775/71‑1 (MQ=255)
                                                      gggCGTGTATTTTGCCGTGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCtt               <  1:310825/71‑1 (MQ=255)
                                                      gggCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCtt               <  1:635560/71‑1 (MQ=255)
                                                         cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGa                          <  1:793820/57‑1 (MQ=255)
                                                         cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg            <  1:833616/71‑1 (MQ=255)
                                                         cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg            <  1:3750849/71‑1 (MQ=255)
                                                          gtgtATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg            <  1:3702515/70‑1 (MQ=255)
                                                             tATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGa                          >  1:1773127/1‑53 (MQ=255)
                                                             tATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcg         >  1:2694424/1‑70 (MQ=255)
                                                             tATTTTGCCATGAACAGACCGATGGCGTTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgc        >  1:3099740/1‑71 (MQ=255)
                                                              aTTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGt                    <  1:3126743/58‑1 (MQ=255)
                                                              aTTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgc        <  1:3277909/70‑1 (MQ=255)
                                                               ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg       >  1:2200650/1‑70 (MQ=255)
                                                               ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCgcgc      <  1:2028226/71‑1 (MQ=255)
                                                                tttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTTGTCGCTTATGCcgcg       >  1:3863631/1‑69 (MQ=255)
                                                                tttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg       >  1:2864977/1‑69 (MQ=255)
                                                                 ttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGAtgtg                      >  1:2577348/1‑53 (MQ=255)
                                                                   gCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggt  >  1:3931118/1‑71 (MQ=255)
                                                                     cATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTc                   <  1:1628616/52‑1 (MQ=255)
                                                                     cATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCt                <  1:3705852/55‑1 (MQ=255)
                                                                     cATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCt                <  1:3688162/55‑1 (MQ=255)
                                                                     cATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCt             <  1:1351171/58‑1 (MQ=255)
                                                                     cATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCt             <  1:870052/58‑1 (MQ=255)
                                                                     |                                                                    
AATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGT  >  minE/385897‑386034

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