Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,994 T→A 29.4% V483V (GTA→GTT) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,9940TA29.4% 225.0 / 119.9 146V483V (GTA→GTT) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (62/41);  new base A (30/13);  total (92/54)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 147 total aligned reads displayed.

TTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAA  >  minE/1932925‑1933063
                                                                     |                                                                     
ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                                      >  1:2810753/1‑71 (MQ=255)
ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                                      >  1:80001/1‑71 (MQ=255)
  aTTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                                     <  1:4162378/70‑1 (MQ=255)
  aTTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                                     <  1:1140601/70‑1 (MQ=255)
   ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                                   >  1:2969026/1‑71 (MQ=255)
        ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACg                                                               <  1:469845/70‑1 (MQ=255)
        ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACg                                                               <  1:5317859/70‑1 (MQ=255)
           gACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                            >  1:4596429/1‑70 (MQ=255)
             cGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGa                                                          >  1:1911385/1‑70 (MQ=255)
                 tCGCCGGTAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAc                                                                     <  1:2761608/55‑1 (MQ=255)
                  cGCCGGGAAGAGAATAACGTCGCGAATAGTATGACTGTTGTAAACAGCATAATCATTCGGt                                                             >  1:3371417/1‑61 (MQ=255)
                   gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAtt                                                                 >  1:1258370/1‑56 (MQ=255)
                   gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                                    >  1:1121669/1‑70 (MQ=255)
                   gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                                    >  1:1549438/1‑70 (MQ=255)
                   gCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt                                                   >  1:596568/1‑71 (MQ=255)
                    ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                                 >  1:615725/1‑56 (MQ=255)
                    ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                                 >  1:1964495/1‑56 (MQ=255)
                     cGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGa                                                          >  1:2007391/1‑62 (MQ=255)
                      gggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCaa                                                    >  1:166024/1‑67 (MQ=255)
                      gggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCaa                                                    >  1:538936/1‑67 (MQ=255)
                      gggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                                <  1:1912456/71‑1 (MQ=255)
                       ggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                                      <  1:1215495/64‑1 (MQ=255)
                         aagaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                                                >  1:4509101/4‑68 (MQ=255)
                         aagaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                                              <  1:4105807/67‑1 (MQ=255)
                          agaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGa                                            >  1:4320936/3‑71 (MQ=255)
                           cAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                           >  1:2554270/1‑71 (MQ=255)
                            aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                                      >  1:3687344/1‑43 (MQ=255)
                            aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                                      >  1:3349004/1‑43 (MQ=255)
                            aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                                      >  1:4210291/1‑43 (MQ=255)
                             gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                         >  1:539040/1‑71 (MQ=255)
                             gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                         >  1:4361469/1‑71 (MQ=255)
                             gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                         >  1:2071754/1‑71 (MQ=255)
                             gAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                         >  1:602986/1‑71 (MQ=255)
                             gAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                                           >  1:591399/1‑69 (MQ=255)
                              aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct                                         >  1:1113569/1‑70 (MQ=255)
                              aaTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAccag                                        >  1:3842704/1‑69 (MQ=255)
                               aTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGa                                            >  1:4995649/1‑66 (MQ=255)
                               aTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                                           >  1:4415567/1‑67 (MQ=255)
                                tAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                                      <  1:5687052/71‑1 (MQ=255)
                                  aCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAtt                                                  <  1:5310255/57‑1 (MQ=255)
                                  aCGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAtac                                                 >  1:3015437/1‑56 (MQ=255)
                                   cGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAtt                                                  <  1:5224403/56‑1 (MQ=255)
                                    gTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCa                                                     >  1:4889263/1‑52 (MQ=255)
                                     tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcgg                                 <  1:5073843/71‑1 (MQ=255)
                                     tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                                      <  1:4781095/50‑1 (MQ=255)
                                     tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                                      <  1:1028680/50‑1 (MQ=255)
                                     tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                                      <  1:1262694/50‑1 (MQ=255)
                                     tCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAcc                                                      <  1:405258/50‑1 (MQ=255)
                                      cgcgAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACcagcg                                      <  1:4651124/65‑2 (MQ=255)
                                       gcgGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTccc                                               <  1:1094365/55‑1 (MQ=255)
                                         ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAGTTCCCAGACCTGCTGTcggcgg                              <  1:2206873/70‑1 (MQ=255)
                                         ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                                        <  1:2438293/60‑1 (MQ=255)
                                         ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt                             <  1:834487/71‑1 (MQ=255)
                                          gATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                                                >  1:2037821/1‑35 (MQ=255)
                                          gATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                                      >  1:1153979/1‑61 (MQ=255)
                                          gATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt                             >  1:1655341/1‑70 (MQ=255)
                                           atagTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTtgg                                 >  1:3781893/4‑65 (MQ=255)
                                              gTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc                         >  1:1587734/1‑70 (MQ=255)
                                               tATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACcagcgg                                     <  1:2880802/57‑3 (MQ=255)
                                                 tGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCAt                      <  1:5646041/70‑1 (MQ=255)
                                                   gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGtt                   >  1:259517/1‑71 (MQ=255)
                                                    cTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc                  <  1:2287609/71‑1 (MQ=255)
                                                    cTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACcagcgg                                     >  1:3006644/1‑50 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTTCTCTCGGCGGTAAGCCATGTTCCa                <  1:227321/70‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTa                            >  1:4117786/1‑58 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTa                            >  1:4531051/1‑58 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc                 <  1:692037/69‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc                 <  1:5648073/69‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc                 <  1:2511018/69‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc                 <  1:1677124/69‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTcc                 <  1:498485/69‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa                <  1:2188498/70‑1 (MQ=255)
                                                       ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAg               >  1:249944/1‑71 (MQ=255)
                                                         ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                           >  1:639598/1‑41 (MQ=255)
                                                         ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg             >  1:3662225/1‑71 (MQ=255)
                                                         ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg             >  1:4782943/1‑71 (MQ=255)
                                                          gTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg             <  1:4421421/70‑1 (MQ=255)
                                                           tGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg           >  1:837025/1‑71 (MQ=255)
                                                            aaaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAgcgcag          >  1:2897954/2‑69 (MQ=255)
                                                            aaaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAgcgcag          >  1:3994867/2‑69 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGAc                                           >  1:4682425/1‑37 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGAc                                           >  1:4576652/1‑37 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGAc                                           >  1:4574457/1‑37 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGAc                                           >  1:4565460/1‑37 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt         <  1:4710148/71‑1 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt         <  1:244594/71‑1 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt         <  1:3268170/71‑1 (MQ=255)
                                                             aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt         <  1:3781630/71‑1 (MQ=255)
                                                             aaCAGCATAATCATTCGGTCGATACCAATACCCAGa                                            >  1:4090682/1‑36 (MQ=16)
                                                             aaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTtggcggca                            >  1:2556513/1‑50 (MQ=255)
                                                              aCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTaa                           <  1:4545840/52‑1 (MQ=255)
                                                               cAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcgg                                 >  1:1046314/1‑45 (MQ=255)
                                                               cAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCAt                      >  1:5240619/1‑56 (MQ=255)
                                                               cAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCAt                      >  1:4860255/1‑56 (MQ=255)
                                                                aGCATTACCATACGGTCGATACCAATTCCCAGACct                                         >  1:2953900/1‑36 (MQ=255)
                                                                 gCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg           >  1:3103808/1‑65 (MQ=255)
                                                                 gCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg           >  1:1486122/1‑65 (MQ=255)
                                                                 gCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg           >  1:1454453/1‑65 (MQ=255)
                                                                  cATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACg     >  1:379263/1‑70 (MQ=255)
                                                                     tACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTaa  >  1:1038074/1‑70 (MQ=255)
                                                                     |                                                                     
TTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAA  >  minE/1932925‑1933063

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: