Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,161 G→A 10.4% R94R (CGC→CGT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,1610GA10.4% 269.3 / 17.2 96R94R (CGC→CGTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (45/41);  new base A (7/3);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGC  >  minE/1934092‑1934228
                                                                     |                                                                   
gTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGc                                                                    <  1:1932544/71‑1 (MQ=255)
gTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGc                                                                    <  1:1793692/71‑1 (MQ=255)
     aTTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGAcc                                                                >  1:143681/1‑70 (MQ=255)
     aTTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCg                                                               >  1:869084/1‑71 (MQ=255)
     aTTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCg                                                               >  1:770755/1‑71 (MQ=255)
       tttttGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGc                                                              >  1:1670600/1‑70 (MQ=255)
       tttttGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGc                                                              >  1:3141339/1‑70 (MQ=255)
            gAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACg                                                         <  1:1929597/70‑1 (MQ=255)
            gAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGt                                                        >  1:753054/1‑71 (MQ=255)
             aaCTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGt                                                        <  1:305061/70‑1 (MQ=255)
                tgatCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCACCGACATCCTg                                                    >  1:2110839/4‑71 (MQ=255)
                tGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTg                                                    >  1:427440/1‑71 (MQ=255)
                  cTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTg                                                    <  1:78737/69‑1 (MQ=255)
                  cTCGTTATAAACGCCTTCCGGGAGATCCTCACGGGCAACTTACAGCTGAATGCGACCGCCAACGTCCTg                                                    <  1:1628814/69‑1 (MQ=255)
                   tCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCa                                                  <  1:3229031/70‑1 (MQ=255)
                   tCGTTATAAACACCTTCTGGCAGGCTATCTCTTGCAACGTACAGTTGAATACGGCCACCGACATCCTGCaa                                                 >  1:3533858/1‑70 (MQ=255)
                    cGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAgg                                                >  1:605132/1‑71 (MQ=255)
                      ttATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGt                                              <  1:3138423/71‑1 (MQ=255)
                       tataAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGt                                              >  1:84639/1‑70 (MQ=255)
                        ataAACGCCTTCCGGGAGATTGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTaa                                            >  1:1146836/1‑71 (MQ=255)
                           aaCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         <  1:1110733/71‑1 (MQ=255)
                           aaCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         <  1:1251231/71‑1 (MQ=255)
                           aaCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         <  1:414129/71‑1 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGAcc                                                                <  1:3948456/47‑1 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGc                                                              >  1:2067766/1‑49 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCaa                                                           >  1:2360674/1‑52 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCaa                                                           >  1:20534/1‑52 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCaa                                                           >  1:3085135/1‑52 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCt                                                     >  1:1558403/1‑58 (MQ=255)
                            aCGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         >  1:2757308/1‑70 (MQ=255)
                             cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         >  1:2351674/1‑69 (MQ=255)
                             cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         >  1:1195468/1‑69 (MQ=255)
                             cGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         >  1:3482992/1‑69 (MQ=255)
                                cTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGAc                                    >  1:3822468/1‑71 (MQ=255)
                                  tCCGGGAGGTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCa                                                  >  1:1264911/1‑55 (MQ=255)
                                  tCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCa                                                  >  1:1887158/1‑55 (MQ=255)
                                  tCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCa                                                  >  1:4883/1‑55 (MQ=255)
                                   ccGGGAGATCGTCACGGGCAACGTACAGCTGTATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                  <  1:1339389/70‑1 (MQ=255)
                                   ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACg                                   <  1:1693437/69‑1 (MQ=255)
                                   ccGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                 >  1:1586997/1‑71 (MQ=255)
                                     gggAGATCGTCACGGGCAACGTACAGCTGAATGCg                                                                   <  1:1532557/35‑1 (MQ=255)
                                     gggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         <  1:3897061/61‑1 (MQ=255)
                                     gggAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCtt                                >  1:2888995/1‑70 (MQ=255)
                                         gaTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAcc                            >  1:1706607/1‑70 (MQ=255)
                                         gaTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAcc                            >  1:3931835/1‑70 (MQ=255)
                                         gaTCGTCACGGGCAACGTACAGCTGAAAGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAccc                           <  1:436605/71‑1 (MQ=255)
                                          aTCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAccc                           <  1:2413681/70‑1 (MQ=255)
                                          aTCGTCACGGGCAACGTACAGCTGAATACGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCa                          >  1:3835260/1‑71 (MQ=255)
                                           tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                       <  1:1893285/57‑1 (MQ=255)
                                           tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                       <  1:3817657/57‑1 (MQ=255)
                                           tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                       <  1:1682024/57‑1 (MQ=255)
                                           tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                         <  1:3872507/71‑1 (MQ=255)
                                           tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                         <  1:46477/71‑1 (MQ=255)
                                           tCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                         <  1:654561/71‑1 (MQ=255)
                                            cGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                         >  1:1667779/1‑70 (MQ=255)
                                            cGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                         >  1:2058786/1‑70 (MQ=255)
                                            cGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCat                         >  1:1019782/1‑70 (MQ=255)
                                             gTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTAccc                           >  1:1479895/1‑67 (MQ=255)
                                              tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAggg                                               <  1:2661970/46‑1 (MQ=255)
                                              tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGa                                         <  1:2197201/52‑1 (MQ=255)
                                              tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                       >  1:874430/1‑70 (MQ=255)
                                              tCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCataa                       >  1:2383178/1‑70 (MQ=255)
                                               cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaa                                        >  1:3682362/1‑52 (MQ=255)
                                               cACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaa                                        >  1:3636070/1‑52 (MQ=255)
                                                aCGTGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                    >  1:1721183/1‑71 (MQ=255)
                                                aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTa                                             >  1:801775/1‑46 (MQ=255)
                                                aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                       >  1:314147/1‑52 (MQ=255)
                                                aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                    <  1:3568748/71‑1 (MQ=255)
                                                aCGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATac                    <  1:2946840/71‑1 (MQ=255)
                                                 cttgCAACGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATac                    >  1:1342083/4‑70 (MQ=255)
                                                  gggCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCata                        >  1:2664461/1‑65 (MQ=255)
                                                  gggCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATTATacga                  <  1:3212337/71‑1 (MQ=255)
                                                    gCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAAc                                           >  1:373295/1‑44 (MQ=255)
                                                    gCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgacg                <  1:1001202/71‑1 (MQ=255)
                                                     cAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGaaa                                       <  1:1176573/47‑1 (MQ=255)
                                                        cGTACAGTTGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGtc            <  1:3945102/71‑1 (MQ=255)
                                                        cGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCt                                 >  1:2175408/1‑50 (MQ=255)
                                                           aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGc                                  >  1:1685103/1‑46 (MQ=255)
                                                           aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat          >  1:3091605/1‑70 (MQ=255)
                                                           aCAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcat          >  1:2562071/1‑70 (MQ=255)
                                                            cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc         <  1:1719053/70‑1 (MQ=255)
                                                            cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc         <  1:1020139/70‑1 (MQ=255)
                                                            cAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGtcatc         <  1:146099/70‑1 (MQ=255)
                                                             aGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCTATAATACGACGGGTcatcat       <  1:1842042/71‑1 (MQ=255)
                                                               ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcattc     >  1:2098882/2‑69 (MQ=255)
                                                               ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcattc     >  1:1474238/2‑69 (MQ=255)
                                                               ttGAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTcatcat       <  1:3442194/68‑1 (MQ=255)
                                                               cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACggg              <  1:2773660/62‑1 (MQ=255)
                                                               cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATg      <  1:3422759/70‑1 (MQ=255)
                                                               cTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATg      <  1:3424822/70‑1 (MQ=255)
                                                                tGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATacgac                 <  1:2419705/58‑1 (MQ=255)
                                                                 gAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCgg   >  1:613890/1‑71 (MQ=255)
                                                                 gAATACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg   >  1:2435959/1‑71 (MQ=255)
                                                                  aaTGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCgg   <  1:1966525/70‑1 (MQ=255)
                                                                  aaTACGGCCACCGACATCCTGCAACGTTACAAAGGAGGCTTTCCCCATGATACGACGGGTCATCATTCgg   <  1:1549006/70‑1 (MQ=255)
                                                                   aTGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGc  <  1:3845564/70‑1 (MQ=255)
                                                                     |                                                                   
GTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGCGACCGCCAACGTCCTGCAGGGTAACGAAAGACGCTTTACCCATAATACGACGGGTCATCATGCGGC  >  minE/1934092‑1934228

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: