Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,268 A→G 26.7% D28D (GAT→GAC 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,2680AG26.7% 179.3 / 67.8 101D28D (GAT→GACcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (14/13);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 109 total aligned reads displayed.

ACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTT  >  minE/1211198‑1211338
                                                                      |                                                                      
aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca                                                                        >  1:2657771/1‑71 (MQ=255)
aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca                                                                        >  1:3192446/1‑71 (MQ=255)
 cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca                                                                        >  1:3527977/1‑70 (MQ=255)
    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGTAGAAGTGTACGAACACATCt                                                                     <  1:2715024/70‑1 (MQ=255)
    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                                                    >  1:3407707/1‑71 (MQ=255)
    aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGAGTACGAACACATCt                                                                     >  1:103504/1‑70 (MQ=255)
     aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAGGTGTACGAACACATCttt                                                                   <  1:947240/71‑1 (MQ=255)
      cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACTCATCTTtg                                                                  >  1:1761059/1‑71 (MQ=255)
      cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtg                                                                  >  1:3923645/1‑71 (MQ=255)
       gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                                                 >  1:637407/1‑71 (MQ=255)
          cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                                               <  1:1901325/70‑1 (MQ=255)
            gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                                            >  1:582459/1‑71 (MQ=255)
             aCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                                            <  1:543206/70‑1 (MQ=255)
             aCCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGcc                                                             <  1:1559666/69‑1 (MQ=25)
              ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt                                                           >  1:3633317/1‑70 (MQ=255)
               cTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt                                                           >  1:1249309/1‑69 (MQ=255)
                 tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGcc                                                             >  1:440163/1‑65 (MQ=19)
                 tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTcttc                                                       >  1:3353712/1‑68 (MQ=16)
                   aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                     >  1:2785260/1‑71 (MQ=255)
                   aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                     <  1:77269/71‑1 (MQ=255)
                   aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                     <  1:1982487/71‑1 (MQ=255)
                     ccAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGa                                                   <  1:559992/71‑1 (MQ=255)
                       agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg                                                  >  1:1540361/1‑70 (MQ=255)
                       agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                                                 <  1:2122453/71‑1 (MQ=255)
                         agTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                                               >  1:629527/1‑71 (MQ=255)
                         agTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTaa                                               >  1:2564493/1‑71 (MQ=255)
                           ttttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                             <  1:1300803/71‑1 (MQ=255)
                             ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                             >  1:689997/1‑69 (MQ=255)
                             ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                             >  1:3587239/1‑69 (MQ=255)
                               gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg                                          <  1:2114138/70‑1 (MQ=255)
                                aaaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAcc                                        <  1:647470/68‑1 (MQ=12)
                                 aaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTc                                                                      <  1:103943/38‑1 (MQ=16)
                                 aaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTc                                                                      <  1:1407624/38‑1 (MQ=16)
                                 aaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTtgctg                                                               >  1:2345410/3‑47 (MQ=255)
                                 aaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAccg                                       <  1:737923/69‑2 (MQ=255)
                                 aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCGGCCGGAGTAATGAAGcc                                        >  1:2443707/1‑70 (MQ=255)
                                 aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGcc                                        >  1:217784/1‑70 (MQ=255)
                                  gCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                                                    <  1:2521295/41‑1 (MQ=255)
                                  gCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                                                 <  1:3169696/60‑1 (MQ=255)
                                  accATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTct                                                         >  1:3407354/2‑51 (MQ=255)
                                   ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                     <  1:2626606/71‑1 (MQ=255)
                                    cATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                                                >  1:486188/1‑43 (MQ=255)
                                    cATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                                                >  1:1345245/1‑43 (MQ=255)
                                    cATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                     <  1:247712/70‑1 (MQ=255)
                                      ttACCCTGGATAGCGGAGAAGTGTACGATCACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                                      <  1:2296236/67‑1 (MQ=255)
                                      ttACCCTGGATAGCGGAGAAGTGTACGAACACATCt                                                                     <  1:3368449/36‑1 (MQ=255)
                                      ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                                            <  1:591721/45‑1 (MQ=255)
                                        aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                                                <  1:1464705/39‑1 (MQ=255)
                                        aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                                                <  1:1899480/39‑1 (MQ=255)
                                        aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATc                                                          >  1:1419261/1‑45 (MQ=255)
                                        aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCt                                 <  1:3298623/70‑1 (MQ=255)
                                         cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                <  1:1499862/70‑1 (MQ=255)
                                          ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                              <  1:1762428/71‑1 (MQ=255)
                                           cTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtttg                              >  1:3132258/1‑69 (MQ=255)
                                           cTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCat                                                        >  1:2991813/1‑43 (MQ=12)
                                           cTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg                                                        <  1:3017016/44‑1 (MQ=255)
                                            tGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                <  1:1583846/67‑1 (MQ=255)
                                            tGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGc                                         >  1:525164/1‑58 (MQ=12)
                                            tGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt                               <  1:3433160/68‑1 (MQ=12)
                                            tGGATAGCGGAGAAGTGTACGAACACATCTTTTCTGCCATCAGCCGGAGTAATGAAGCCaaa                                     <  1:70980/62‑1 (MQ=255)
                                            tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                             <  1:3924194/54‑1 (MQ=255)
                                             ggATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                     <  1:3065592/61‑1 (MQ=255)
                                              gATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtttgg                             <  1:3510652/68‑3 (MQ=255)
                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGAc                           <  1:831281/70‑1 (MQ=255)
                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGAc                           <  1:2987128/70‑1 (MQ=255)
                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGACGGAGTAATGAAGCCAAAACCTTTAGAc                           <  1:3581191/70‑1 (MQ=255)
                                              gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGACGGAGTAATGAAGCCAAAACCTTTAGAc                           <  1:3585649/70‑1 (MQ=255)
                                               aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc                         >  1:95075/1‑71 (MQ=255)
                                               aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc                         >  1:855873/1‑71 (MQ=255)
                                               aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc                         >  1:1876307/1‑71 (MQ=255)
                                                  gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                   >  1:3017519/1‑58 (MQ=255)
                                                  gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt                      >  1:74359/1‑71 (MQ=255)
                                                   cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCttt                               >  1:3227927/1‑61 (MQ=255)
                                                   cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCttt                               >  1:2541735/1‑61 (MQ=255)
                                                   cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                     <  1:2159560/71‑1 (MQ=255)
                                                    ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg                                                    <  1:309257/39‑1 (MQ=255)
                                                    agagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATg                                             <  1:1597601/45‑1 (MQ=21)
                                                     gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCgaat                                    >  1:3068446/1‑53 (MQ=18)
                                                      agaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAc                  >  1:2135013/1‑71 (MQ=255)
                                                      agaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                  >  1:3304559/1‑71 (MQ=255)
                                                         aGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCAc               <  1:1681386/71‑1 (MQ=255)
                                                          gtgtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAccga                                      >  1:685861/1‑45 (MQ=255)
                                                          gtgtAAGAACACATCTTTGCTGCCATCAGCCGGAGt                                                 <  1:3463070/36‑1 (MQ=255)
                                                           tgtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                <  1:1521991/52‑1 (MQ=255)
                                                           tgtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACtt             <  1:3909223/71‑1 (MQ=12)
                                                            gtACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCTAAACCTTTAGACTCGTTGAACCACtt             >  1:868425/1‑70 (MQ=255)
                                                            gtACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACtt             >  1:845793/1‑70 (MQ=255)
                                                             tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTaaa                   >  1:578390/1‑63 (MQ=9)
                                                             tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTa            <  1:1194781/70‑1 (MQ=255)
                                                             tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATACtttg                              >  1:883755/1‑51 (MQ=9)
                                                               cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcg         >  1:1793018/1‑70 (MQ=255)
                                                               cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCa                                       >  1:676270/1‑41 (MQ=255)
                                                               cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                   >  1:1429503/1‑45 (MQ=255)
                                                                 aaCACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGa                            >  1:1510217/1‑50 (MQ=255)
                                                                  acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                   >  1:306497/1‑42 (MQ=255)
                                                                  acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc                         <  1:272189/52‑1 (MQ=255)
                                                                   cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc     >  1:3287333/1‑71 (MQ=255)
                                                                    acaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt        <  1:2200473/67‑1 (MQ=255)
                                                                    aaaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc     >  1:263882/3‑70 (MQ=255)
                                                                      aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt  <  1:1820088/71‑1 (MQ=255)
                                                                      |                                                                      
ACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTT  >  minE/1211198‑1211338

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: