Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,985 G→A 26.3% F486F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,9850GA26.3% 211.9 / 82.7 133F486F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (59/39);  new base A (22/13);  total (81/52)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

Only 100 of 136 total aligned reads displayed.

ACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCG  >  minE/1932916‑1933054
                                                                     |                                                                     
aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                                       <  1:2485013/70‑1 (MQ=255)
aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                                       <  1:1874147/70‑1 (MQ=255)
aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                                       >  1:1143899/1‑70 (MQ=255)
        cTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGctt                                                               <  1:905567/70‑3 (MQ=255)
        cTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                               <  1:287165/70‑1 (MQ=255)
        cTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                               <  1:1281469/70‑1 (MQ=255)
        cTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                               <  1:3173011/70‑1 (MQ=255)
         ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACa                                                                   >  1:3424891/1‑65 (MQ=255)
         ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                             >  1:2810753/1‑71 (MQ=255)
         ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                             >  1:80001/1‑71 (MQ=255)
           aTTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                            <  1:1140601/70‑1 (MQ=255)
           aTTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                            <  1:4162378/70‑1 (MQ=255)
            ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                          >  1:2969026/1‑71 (MQ=255)
               tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                                       >  1:4696342/1‑71 (MQ=255)
                 ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACg                                                      <  1:5661682/70‑1 (MQ=255)
                 ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACg                                                      <  1:5317859/70‑1 (MQ=255)
                 ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACg                                                      <  1:469845/70‑1 (MQ=255)
                    gACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGaa                                                                     >  1:1195276/1‑52 (MQ=255)
                    gACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                   >  1:4596429/1‑70 (MQ=255)
                          tCGCCGGTAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                            <  1:2761608/55‑1 (MQ=255)
                          tCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGa                                                                      <  1:2358882/45‑1 (MQ=255)
                          tCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                            <  1:4976200/55‑1 (MQ=255)
                          tCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                            <  1:5327225/55‑1 (MQ=255)
                           cGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:4480796/1‑71 (MQ=255)
                           cGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:214042/1‑71 (MQ=255)
                           cGCCGGGAAGAGAATAACGTCGCGAATAGTATGACTGTTTGTAAACAGCATAATCATTCGGt                                                    >  1:4371568/1‑62 (MQ=255)
                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:4613870/1‑70 (MQ=255)
                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:4542885/1‑70 (MQ=255)
                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:1121669/1‑70 (MQ=255)
                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:1222664/1‑70 (MQ=255)
                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                           >  1:1549438/1‑70 (MQ=255)
                            gCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACa                                                                   >  1:3372839/1‑46 (MQ=255)
                            gCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACa                                                                   >  1:1871851/1‑46 (MQ=255)
                            gCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACa                                                                   >  1:195424/1‑46 (MQ=255)
                            gCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt                                          >  1:596568/1‑71 (MQ=255)
                            gCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt                                          >  1:2515203/1‑71 (MQ=255)
                             ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                                  >  1:5675664/1‑46 (MQ=255)
                             ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                                  >  1:2576988/1‑46 (MQ=255)
                             ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                        >  1:615725/1‑56 (MQ=255)
                              cGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGa                                                 >  1:2007391/1‑62 (MQ=255)
                               gggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCaa                                           >  1:166024/1‑67 (MQ=255)
                               gggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                       <  1:1912456/71‑1 (MQ=255)
                                ggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                             <  1:1215495/64‑1 (MQ=255)
                                ggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAccc                                      <  1:2474814/71‑1 (MQ=255)
                                  aagaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                                       >  1:4509101/4‑68 (MQ=255)
                                  aagaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                                     <  1:4105807/67‑1 (MQ=255)
                                  aaCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                       >  1:605765/1‑68 (MQ=255)
                                   aCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGa                                   >  1:4756273/1‑71 (MQ=255)
                                    cAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                  >  1:818567/1‑71 (MQ=255)
                                    cAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                  >  1:2554270/1‑71 (MQ=255)
                                    cAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                  >  1:1415903/1‑71 (MQ=255)
                                     aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                             >  1:3349004/1‑43 (MQ=255)
                                     aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                             >  1:4210291/1‑43 (MQ=255)
                                     aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                             >  1:3687344/1‑43 (MQ=255)
                                     aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                   >  1:2029828/1‑53 (MQ=255)
                                     aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTccc                                      <  1:4848811/66‑1 (MQ=255)
                                      gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                >  1:1768795/1‑71 (MQ=255)
                                      gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                >  1:2071754/1‑71 (MQ=255)
                                      gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                >  1:539040/1‑71 (MQ=255)
                                      gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                >  1:4361469/1‑71 (MQ=255)
                                      gAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                                  >  1:591399/1‑69 (MQ=255)
                                       aaTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAccag                               >  1:3842704/1‑69 (MQ=255)
                                        aTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGa                                   >  1:4995649/1‑66 (MQ=255)
                                         tAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                             <  1:3501053/71‑1 (MQ=255)
                                           aCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAtt                                         <  1:5310255/57‑1 (MQ=255)
                                            cGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAtt                                         <  1:5224403/56‑1 (MQ=255)
                                             gTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCa                                            >  1:2743370/1‑52 (MQ=255)
                                              tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcgg                        <  1:5073843/71‑1 (MQ=255)
                                              tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                             <  1:1262694/50‑1 (MQ=255)
                                              tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                             <  1:405258/50‑1 (MQ=255)
                                              tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                             <  1:1028680/50‑1 (MQ=255)
                                              tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACcagcg                             <  1:2874951/66‑2 (MQ=255)
                                                gcgGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTccc                                      <  1:1094365/55‑1 (MQ=255)
                                                   gATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt                    >  1:1655341/1‑70 (MQ=255)
                                                   gATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt                    >  1:561937/1‑70 (MQ=255)
                                                    atagTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTtgg                        >  1:3781893/4‑65 (MQ=255)
                                                       gTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc                >  1:3022853/1‑70 (MQ=255)
                                                       gTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc                >  1:1587734/1‑70 (MQ=255)
                                                        tATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACcagcgg                            <  1:3806969/57‑3 (MQ=255)
                                                          tGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCAt             <  1:5646041/70‑1 (MQ=255)
                                                            gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGtt          >  1:4266918/1‑71 (MQ=255)
                                                            gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGt           >  1:5237861/1‑70 (MQ=255)
                                                             cTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc         <  1:2287609/71‑1 (MQ=255)
                                                             cTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc         <  1:2899381/71‑1 (MQ=255)
                                                             cTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACcagcgg                            >  1:3006644/1‑50 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTTCTCTCGGCGGTAAGCCATGTTCCa       <  1:227321/70‑1 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc        <  1:498485/69‑1 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc        <  1:692037/69‑1 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa       <  1:2188498/70‑1 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAg      >  1:106959/1‑71 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAg      <  1:1957555/71‑1 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAg      >  1:249944/1‑71 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGAGGTa                   >  1:163794/1‑58 (MQ=255)
                                                                  ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                  >  1:639598/1‑41 (MQ=255)
                                                                  ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg    >  1:4782943/1‑71 (MQ=255)
                                                                  ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg    >  1:3662225/1‑71 (MQ=255)
                                                                  ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg    >  1:409699/1‑71 (MQ=255)
                                                                   gTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg    <  1:4421421/70‑1 (MQ=255)
                                                                    tGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg  <  1:594551/71‑1 (MQ=255)
                                                                    tGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg  >  1:837025/1‑71 (MQ=255)
                                                                     |                                                                     
ACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCG  >  minE/1932916‑1933054

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: