Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,837 C→G 11.7% T349T (ACC→ACG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8370CG11.7% 157.9 / 11.1 60T349T (ACC→ACGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base G (2/5);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCT  >  minE/385770‑385907
                                                                   |                                                                      
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)
           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:2903705/1‑71 (MQ=255)
                                   tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                  >  1:969233/1‑71 (MQ=255)
                                   tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                  >  1:2899356/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)
                                                              gggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg        >  1:1923530/1‑70 (MQ=255)
                                                               ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg        <  1:2062311/69‑1 (MQ=255)
                                                               ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGa       <  1:3520558/70‑1 (MQ=255)
                                                                   cACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCt  <  1:3416980/71‑1 (MQ=255)
                                                                   |                                                                      
CGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCT  >  minE/385770‑385907

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: