Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,471 A→G 5.2% M561V (ATG→GTG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,4710AG 5.2% 222.1 / 2.1 77M561V (ATG→GTG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (25/48);  new base G (2/2);  total (27/50)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGG  >  minE/386405‑386540
                                                                  |                                                                     
gATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGgc                                                                   <  1:2984461/71‑1 (MQ=255)
    cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                                                                <  1:5492450/70‑1 (MQ=255)
    cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                                                                <  1:5444251/70‑1 (MQ=255)
     gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGAt                                                                      <  1:2398335/63‑1 (MQ=255)
     gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                                                                <  1:2758244/69‑1 (MQ=255)
     gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                                                                <  1:146810/69‑1 (MQ=255)
      cTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGAt                                                                      <  1:490666/62‑1 (MQ=255)
         ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATg                                                            <  1:4941133/69‑1 (MQ=255)
         ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGc                                                           <  1:1693235/70‑1 (MQ=255)
                ttAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGcc                                                   <  1:4829369/71‑1 (MQ=255)
                  aaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgt                                                 >  1:3141143/1‑71 (MQ=255)
                  aaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgt                                                 >  1:1956186/1‑71 (MQ=255)
                   acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCtgt                                                      >  1:5453234/1‑65 (MQ=255)
                   acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtg                                                >  1:22249/1‑71 (MQ=255)
                     acTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtggc                                              <  1:2690552/71‑2 (MQ=255)
                      ctctGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtggtg                                             <  1:4502642/71‑1 (MQ=255)
                         tgtgGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGgtgt                                            <  1:1909524/69‑1 (MQ=255)
                         tgtgGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGgtgt                                            <  1:1902847/69‑1 (MQ=255)
                           tgGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                         <  1:2220919/70‑1 (MQ=255)
                           tgGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTgttgtt                                        <  1:3904189/71‑1 (MQ=255)
                             gCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTGTGTGCCGTGGTGTTGTTc                                       >  1:1524396/1‑70 (MQ=255)
                             gCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTc                                       <  1:1081807/70‑1 (MQ=255)
                              ccGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCaa                                     <  1:1004850/71‑1 (MQ=255)
                              ccGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCaa                                     <  1:5119035/71‑1 (MQ=255)
                               cGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtg                                                <  1:1543935/59‑1 (MQ=255)
                               cGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtg                                                <  1:4678787/59‑1 (MQ=255)
                               cGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCaa                                     <  1:897326/70‑1 (MQ=255)
                                 cTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGAt                                  <  1:4079546/71‑1 (MQ=255)
                                 cTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGAt                                  <  1:4090526/71‑1 (MQ=255)
                                 cTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGAt                                  <  1:1507950/71‑1 (MQ=255)
                                    tttGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGa                                <  1:3445311/70‑1 (MQ=255)
                                    tttGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGa                                <  1:4236191/70‑1 (MQ=255)
                                      cggTATCTCCAACCAGATGCTGGCAGCCGTAGCGCTGATACTGGGCACCGTTGTGCTGATTAAGATGaagc                             >  1:42089/2‑69 (MQ=255)
                                       ggTATCTCCAACCAGATGCTGGCAGCCGTAGCGCTGGTACTGGGCACCGTTGTGCTGATTAAGATGaagc                             <  1:3029975/70‑3 (MQ=255)
                                       ggTATCTCCAACCAGATGCTGGCAGCCGTAGCGCTGGTACTGGGCACCGTTGTGCTGATTAAGATGAAGCg                            >  1:4953888/1‑71 (MQ=255)
                                         tATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGAt                                  <  1:1215862/63‑1 (MQ=255)
                                         tATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGAt                                  <  1:3022940/63‑1 (MQ=255)
                                         tATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGt                           >  1:1260441/1‑70 (MQ=255)
                                         tATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGt                           >  1:4084974/1‑70 (MQ=255)
                                          aTTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTc                          >  1:169844/1‑70 (MQ=255)
                                           tctccAACCAGATGCTGGCAGCCGTAGCGCTGGTACTGGGCACCGTTGTGCTGATTAAGATGaa                               <  1:83311/61‑1 (MQ=255)
                                             gCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACg                      >  1:2495636/1‑71 (MQ=255)
                                             gCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACg                      >  1:5034588/1‑71 (MQ=255)
                                             gCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACg                      >  1:2519130/1‑71 (MQ=255)
                                               cAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTTGTGTTGTTCAAGATGa                                <  1:3648264/59‑1 (MQ=255)
                                               cAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGa                                <  1:2647953/59‑1 (MQ=255)
                                                 aCCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTAc                  <  1:4907755/71‑1 (MQ=255)
                                                 aCCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTAc                  <  1:3394191/71‑1 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                              >  1:1921479/1‑58 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                              >  1:38327/1‑58 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                              >  1:3826783/1‑58 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                              >  1:3349362/1‑58 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                              >  1:883964/1‑58 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGaaa                              >  1:399368/1‑58 (MQ=255)
                                                  ccAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACg                 >  1:4225943/1‑71 (MQ=255)
                                                    aGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                         <  1:3550681/45‑1 (MQ=255)
                                                    aGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                         <  1:3463705/45‑1 (MQ=255)
                                                    aGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGtgttgt                                         <  1:4303031/45‑1 (MQ=255)
                                                     gATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGa                                <  1:665549/53‑1 (MQ=255)
                                                     gATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGa                                <  1:5393335/53‑1 (MQ=255)
                                                     gATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGa                                <  1:2260098/53‑1 (MQ=255)
                                                      aTGCTGGCAGGGATGGCGCTGATGCTGTGTGCCGTGGTGTTGTTCAAGATg                                 <  1:2513001/51‑1 (MQ=255)
                                                      aTGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATg                                 <  1:709535/51‑1 (MQ=255)
                                                      aTGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATg                                 <  1:3714851/51‑1 (MQ=255)
                                                        gCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTa                   >  1:2188059/1‑63 (MQ=255)
                                                          tGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATTAAACGTCAACGTTACGCCTGGGTg         >  1:5646910/1‑71 (MQ=255)
                                                          tGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTc                          >  1:4326924/1‑54 (MQ=255)
                                                          tGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTg         >  1:4514974/1‑71 (MQ=255)
                                                          tGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTg         >  1:439857/1‑71 (MQ=255)
                                                            gCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGAt                                  <  1:496917/44‑1 (MQ=255)
                                                            gCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTg         <  1:4342024/69‑1 (MQ=255)
                                                               gggATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCt    <  1:1337188/71‑1 (MQ=255)
                                                               gggATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCt    <  1:532881/71‑1 (MQ=255)
                                                               gggATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCt    <  1:3191149/71‑1 (MQ=255)
                                                                 gATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGtt                    >  1:4364511/1‑53 (MQ=255)
                                                                 gATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTa                   >  1:1600408/1‑54 (MQ=255)
                                                                 gATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  <  1:5376025/71‑1 (MQ=255)
                                                                  aTGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  >  1:1027559/1‑70 (MQ=255)
                                                                  aTGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  >  1:2022396/1‑70 (MQ=255)
                                                                  aTGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  >  1:3415129/1‑70 (MQ=255)
                                                                  aTGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  >  1:5302921/1‑70 (MQ=255)
                                                                  aTGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  >  1:4361443/1‑70 (MQ=255)
                                                                  aTGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTgg  >  1:4953811/1‑70 (MQ=255)
                                                                  |                                                                     
GATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGG  >  minE/386405‑386540

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: