Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,831 T→C 20.3% S347S (TCT→TCC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8310TC20.3% 118.0 / 25.8 59S347S (TCT→TCCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base C (3/9);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATG  >  minE/385766‑385900
                                                                 |                                                                     
ttCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCggta                                                                   <  1:516376/70‑3 (MQ=255)
 tCCCGTTCATGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAc                                                                  >  1:2241936/1‑70 (MQ=255)
   ccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                <  1:2624561/70‑3 (MQ=255)
   ccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                <  1:198973/70‑3 (MQ=255)
    cGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAc                                                               >  1:3610705/1‑70 (MQ=255)
    cGTTCCTGTTGATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTacgacg                                                              >  1:3111843/1‑71 (MQ=255)
         cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAGGCCGaa                                                         <  1:513094/71‑1 (MQ=255)
           gTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                                           >  1:871837/1‑67 (MQ=255)
           gTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                <  1:345865/62‑3 (MQ=255)
               aTCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCt                                                   <  1:2748866/71‑1 (MQ=255)
               aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                    <  1:3476404/70‑1 (MQ=255)
               aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCt                                                   >  1:793095/1‑71 (MQ=255)
                tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                  <  1:1647770/71‑1 (MQ=255)
                tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                  <  1:652499/71‑1 (MQ=255)
                   ccATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                    <  1:3153357/66‑1 (MQ=255)
                    cATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGGTGATCTCTTCTGGGACCACGc                                                             >  1:973459/1‑56 (MQ=255)
                     aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                                             <  1:983120/71‑1 (MQ=255)
                     aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                                             <  1:609374/71‑1 (MQ=255)
                      tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                                            >  1:3689463/1‑71 (MQ=255)
                        gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatga                                          <  1:2135801/71‑1 (MQ=255)
                        gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                          <  1:1849348/71‑1 (MQ=255)
                         ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                                         <  1:2278416/71‑1 (MQ=255)
                         ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                                         <  1:2030088/71‑1 (MQ=255)
                          cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                        >  1:2297925/1‑71 (MQ=255)
                          cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                        <  1:1613947/71‑1 (MQ=255)
                             tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaacc                                      >  1:3047714/1‑67 (MQ=255)
                             tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaac                                       <  1:3812158/69‑3 (MQ=255)
                             tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGc                                                             <  1:277114/47‑1 (MQ=255)
                             tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                                     <  1:2970180/71‑1 (MQ=255)
                               gcgcGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcga                                              <  1:910831/60‑1 (MQ=255)
                               gcgcGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAg                                   <  1:3158300/71‑1 (MQ=255)
                                    gtgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGt                               >  1:2231017/1‑70 (MQ=255)
                                       tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           >  1:2899356/1‑71 (MQ=255)
                                       tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           >  1:969233/1‑71 (MQ=255)
                                       tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           >  1:2903705/1‑71 (MQ=255)
                                        cTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           >  1:3010349/1‑70 (MQ=255)
                                         tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg                              <  1:1400676/66‑1 (MQ=255)
                                         tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                         >  1:611994/1‑71 (MQ=255)
                                           gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                         <  1:3108896/69‑1 (MQ=255)
                                           gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc                       >  1:2382780/1‑71 (MQ=255)
                                            cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg                              >  1:3858316/1‑63 (MQ=255)
                                            cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCg                      >  1:805865/1‑71 (MQ=255)
                                             ttCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcg                                               >  1:1392380/1‑45 (MQ=255)
                                             ttCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCgg                     >  1:1162794/1‑71 (MQ=255)
                                               ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCTAACGAAGGGCAGGCGTGCTTTAt                        >  1:1428709/1‑66 (MQ=255)
                                               ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                                       <  1:230622/51‑1 (MQ=255)
                                               ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGc                    >  1:819680/1‑70 (MQ=255)
                                                   gcgcTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTACgg               >  1:3022534/1‑71 (MQ=255)
                                                     gcTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg             <  1:3509131/71‑1 (MQ=255)
                                                      cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg            >  1:3016323/1‑71 (MQ=255)
                                                      cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg            <  1:1662811/71‑1 (MQ=255)
                                                      cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg             <  1:1583949/70‑1 (MQ=255)
                                                       tGATCTCTTCTGGGACCACTCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                          <  1:2659821/56‑1 (MQ=255)
                                                       tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                          <  1:1582497/56‑1 (MQ=255)
                                                       tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg            >  1:374929/1‑70 (MQ=255)
                                                        gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                          <  1:88656/39‑1 (MQ=255)
                                                        gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                          <  1:1591502/39‑1 (MQ=255)
                                                         aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTggg           <  1:3076262/69‑1 (MQ=255)
                                                          tctcTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg             >  1:696242/1‑66 (MQ=255)
                                                              ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa    >  1:3139172/1‑71 (MQ=255)
                                                               tctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt       >  1:515839/1‑67 (MQ=255)
                                                                ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc                             >  1:3580757/1‑44 (MQ=255)
                                                                ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATg  >  1:2695806/1‑71 (MQ=255)
                                                                 |                                                                     
TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATG  >  minE/385766‑385900

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: