Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,280 A→G 26.5% D24D (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,2800AG26.5% 160.4 / 65.3 98D24D (GAT→GACcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (34/38);  new base G (17/9);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 103 total aligned reads displayed.

GACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCA  >  minE/1211210‑1211350
                                                                      |                                                                      
gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                                                        >  1:582459/1‑71 (MQ=255)
 aCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                                                        <  1:543206/70‑1 (MQ=255)
  ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt                                                                       >  1:3633317/1‑70 (MQ=255)
   cTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt                                                                       >  1:1249309/1‑69 (MQ=255)
     tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTcttc                                                                   >  1:3353712/1‑68 (MQ=16)
       aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                                 <  1:1982487/71‑1 (MQ=255)
       aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                                 <  1:77269/71‑1 (MQ=255)
       aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                                 >  1:2785260/1‑71 (MQ=255)
        gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg                                                                <  1:1518876/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:3587239/1‑69 (MQ=255)
                 ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                                         >  1:689997/1‑69 (MQ=255)
                 ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                                                        <  1:3211159/70‑1 (MQ=255)
                   gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg                                                      <  1:2114138/70‑1 (MQ=255)
                    aaaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAcc                                                    <  1:647470/68‑1 (MQ=12)
                     aaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAccg                                                   <  1:737923/69‑2 (MQ=255)
                     aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCGGCCGGAGTAATGAAGcc                                                    >  1:2443707/1‑70 (MQ=255)
                     aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGcc                                                    >  1:217784/1‑70 (MQ=255)
                      gCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                                                             <  1:3169696/60‑1 (MQ=255)
                      accATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTct                                                                     >  1:3407354/2‑51 (MQ=255)
                       ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                                 <  1:2626606/71‑1 (MQ=255)
                        cATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                                 <  1:247712/70‑1 (MQ=255)
                          ttACCCTGGATAGCGGAGAAGTGTACGATCACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                                                  <  1:2296236/67‑1 (MQ=255)
                          ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                                                        <  1:591721/45‑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: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:2541735/1‑61 (MQ=255)
                                       cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCttt                                           >  1:3227927/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:1887077/1‑71 (MQ=255)
                                          agaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                              >  1:3516745/1‑71 (MQ=255)
                                          agaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                              >  1:3304559/1‑71 (MQ=255)
                                           gaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAcc                             <  1:699068/71‑1 (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)
                                                  aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAc                      >  1:530704/1‑71 (MQ=255)
                                                   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)
                                                       cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc                 <  1:2904577/71‑1 (MQ=255)
                                                        acaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                    <  1:2200473/67‑1 (MQ=255)
                                                        aaaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc                 >  1:263882/3‑70 (MQ=255)
                                                         cgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCt                                             >  1:1569698/3‑41 (MQ=12)
                                                          gtCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCtt              >  1:1462865/2‑71 (MQ=255)
                                                          aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt              <  1:1820088/71‑1 (MQ=255)
                                                            cTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCtt              >  1:2838813/1‑69 (MQ=255)
                                                            cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGaa                               >  1:1231530/1‑52 (MQ=255)
                                                            cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCttt             >  1:3855000/1‑70 (MQ=255)
                                                               tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTa                        >  1:3592555/1‑56 (MQ=255)
                                                               tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTa                        >  1:111644/1‑56 (MQ=255)
                                                               tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATc         >  1:2031568/1‑71 (MQ=255)
                                                               tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATc         >  1:2793454/1‑71 (MQ=255)
                                                               tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATc         >  1:1210538/1‑71 (MQ=255)
                                                                 ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCtt       <  1:501487/71‑1 (MQ=255)
                                                                 ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCtt       <  1:83884/71‑1 (MQ=255)
                                                                 ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCtt       <  1:1277899/71‑1 (MQ=255)
                                                                  tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                 <  1:1698098/44‑1 (MQ=255)
                                                                     cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGcc   <  1:3587172/71‑1 (MQ=255)
                                                                      aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa  <  1:3058124/71‑1 (MQ=255)
                                                                      |                                                                      
GACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCA  >  minE/1211210‑1211350

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: