Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,202 A→T 8.7% V50V (GTT→GTA cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,211,2020AT 8.7% 192.2 / 9.1 69V50V (GTT→GTAcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (3/3);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTT  >  minE/1211135‑1211272
                                                                   |                                                                      
gATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACg                                                                     >  1:3090982/1‑71 (MQ=255)
gATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACg                                                                     <  1:431366/71‑1 (MQ=255)
 aTTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGt                                                                    <  1:3835973/71‑1 (MQ=255)
     gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                 <  1:1105075/70‑1 (MQ=255)
     gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                 <  1:748010/70‑1 (MQ=255)
     gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                 <  1:879768/70‑1 (MQ=255)
         aGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCt                                                            <  1:3007819/71‑1 (MQ=255)
          gATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCtt                                                           <  1:935294/71‑1 (MQ=255)
            tAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCtt                                                           >  1:2449383/1‑69 (MQ=255)
             aGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGa                                                               >  1:3587162/1‑64 (MQ=255)
               cTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGa                                                               <  1:853926/62‑1 (MQ=255)
                    aCGTTAACAGCTGCCGGACCTTTGTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGt                                                 >  1:3797890/1‑71 (MQ=255)
                    aCGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGt                                                 >  1:2123489/1‑71 (MQ=255)
                        tAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCa                                                         <  1:1509052/59‑1 (MQ=255)
                         aaCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                             <  1:873703/70‑1 (MQ=255)
                         aaCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                             <  1:1399248/70‑1 (MQ=255)
                          aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                             >  1:1160875/1‑69 (MQ=255)
                          aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                           <  1:359650/71‑1 (MQ=255)
                           cAGCTGCCGGACCTTTCTGGCCTTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGtt                                                <  1:3585270/65‑1 (MQ=255)
                           cAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGa                                            >  1:1860093/1‑69 (MQ=255)
                           cAGCTGCCGGACCTTTCTGGCCGTCATGAAATTAGAACTCAACGt                                                                    >  1:2425342/1‑45 (MQ=255)
                              cTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCaa                                                                       >  1:2057414/1‑39 (MQ=255)
                              cTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                       >  1:1055219/1‑71 (MQ=255)
                               tGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                      >  1:886128/1‑71 (MQ=255)
                               tGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                      >  1:2308363/1‑71 (MQ=255)
                                 ccGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                      <  1:3759316/69‑1 (MQ=255)
                                 ccGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTa                                    >  1:2699652/1‑71 (MQ=255)
                                 ccGGACCTTTCTGGCCGTCCTGAATCTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAtt                                     <  1:531293/70‑1 (MQ=255)
                                 ccGGACCTTTCTGGACGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTa                                    <  1:2236277/71‑1 (MQ=255)
                                   ggACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                       >  1:2800256/1‑66 (MQ=255)
                                    gACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCaga                                                   <  1:1283729/53‑1 (MQ=255)
                                      ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                           >  1:2902704/1‑59 (MQ=255)
                                      ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                               >  1:899542/1‑71 (MQ=255)
                                      ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                               >  1:778458/1‑71 (MQ=255)
                                        tttCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                 <  1:1051872/35‑1 (MQ=255)
                                         ttCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                       >  1:1354117/1‑60 (MQ=255)
                                          tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTc                                                          >  1:3051221/1‑40 (MQ=255)
                                          tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                             <  1:690552/69‑1 (MQ=255)
                                           cTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATa                           <  1:3792426/70‑1 (MQ=255)
                                           cTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATa                           <  1:725742/70‑1 (MQ=255)
                                              gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGc                         <  1:3019907/69‑1 (MQ=255)
                                              gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                        <  1:2478583/70‑1 (MQ=255)
                                                 gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                     >  1:1478149/1‑70 (MQ=255)
                                                    cTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGcc                                        <  1:3159318/48‑1 (MQ=255)
                                                     tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGttt                                               >  1:2751592/1‑40 (MQ=255)
                                                     tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGttt                                               >  1:186154/1‑40 (MQ=255)
                                                     gtgaTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAgtg                <  1:477945/68‑1 (MQ=255)
                                                        aTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTAc             >  1:1059989/1‑71 (MQ=255)
                                                         tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                               >  1:217490/1‑52 (MQ=255)
                                                          ttCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTa              <  1:706655/68‑1 (MQ=255)
                                                           tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                             <  1:1915725/36‑1 (MQ=255)
                                                           tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                   >  1:156016/1‑62 (MQ=255)
                                                           tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                   >  1:2571572/1‑62 (MQ=255)
                                                           tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGa           <  1:176980/70‑1 (MQ=255)
                                                            cGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGaa          >  1:1637821/1‑70 (MQ=255)
                                                             gAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAac         <  1:1652031/70‑1 (MQ=255)
                                                             gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                         <  1:1687002/38‑1 (MQ=255)
                                                             gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                         <  1:1179802/38‑1 (MQ=255)
                                                             gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                         <  1:1420971/38‑1 (MQ=255)
                                                              aaCTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAacac       >  1:372179/1‑71 (MQ=255)
                                                              aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                               >  1:1054772/1‑47 (MQ=255)
                                                              aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac       <  1:24290/71‑1 (MQ=255)
                                                              aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca        >  1:2495753/1‑70 (MQ=255)
                                                               aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac       <  1:3755859/70‑1 (MQ=255)
                                                               aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca      >  1:3192446/1‑71 (MQ=255)
                                                               aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca      >  1:2657771/1‑71 (MQ=255)
                                                                cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                  <  1:2087680/42‑1 (MQ=255)
                                                                cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                        <  1:933161/52‑1 (MQ=255)
                                                                cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca      >  1:3527977/1‑70 (MQ=255)
                                                                 tCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                   <  1:605378/56‑1 (MQ=255)
                                                                   aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                  >  1:1478301/1‑39 (MQ=255)
                                                                   aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGTAGAAGTGTACGAACACATCt   <  1:2715024/70‑1 (MQ=255)
                                                                   aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt  >  1:3407707/1‑71 (MQ=255)
                                                                   aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGAGTACGAACACATCt   >  1:103504/1‑70 (MQ=255)
                                                                   |                                                                      
GATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTT  >  minE/1211135‑1211272

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: