Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,044 G→A 6.0% H467Y (CAT→TAT)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,0440GA 6.0% 331.9 / 7.3 100H467Y (CAT→TAT) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (45/49);  new base A (4/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 103 total aligned reads displayed.

GCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAC  >  minE/1932976‑1933113
                                                                    |                                                                     
gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGtt                                                                     >  1:3394044/1‑71 (MQ=255)
  tGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc                                                                    >  1:2915820/1‑70 (MQ=255)
        cgAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg                                                             <  1:4064093/70‑1 (MQ=255)
         gAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg                                                             >  1:2354297/1‑70 (MQ=255)
          aaCAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt                                                           >  1:2655096/1‑71 (MQ=255)
          aaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCATACCATATTCCAgcgcag                                                            >  1:2150261/1‑68 (MQ=255)
          aaCAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGt                                                           >  1:127449/1‑71 (MQ=255)
            cAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGa                                                         >  1:3129245/1‑71 (MQ=255)
            cAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGa                                                         >  1:3147988/1‑71 (MQ=255)
             aGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGAc                                                        >  1:1200501/1‑71 (MQ=255)
               cATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACg                                                       >  1:770972/1‑70 (MQ=255)
                 ttACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGt                                                      >  1:684961/1‑69 (MQ=255)
                 ttACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTaa                                                    <  1:1201579/71‑1 (MQ=255)
                 ttACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTaa                                                    <  1:2139306/71‑1 (MQ=255)
                 ttACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTaa                                                    >  1:4100059/1‑71 (MQ=255)
                   aCCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                  >  1:153411/1‑71 (MQ=255)
                   aCCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                  >  1:1759097/1‑71 (MQ=255)
                    ccATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATg                                                                       >  1:2943435/1‑49 (MQ=255)
                    ccATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATg                                                                       >  1:1646375/1‑49 (MQ=255)
                    ccATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCt                                                 >  1:174225/1‑71 (MQ=255)
                     cATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCtt                                                <  1:96404/71‑1 (MQ=255)
                      attcGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTc                                               >  1:3705486/4‑71 (MQ=255)
                      aTACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                  <  1:727059/68‑1 (MQ=255)
                      aTACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTc                                               <  1:1332598/71‑1 (MQ=255)
                       tACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATg                                                                       <  1:2368491/46‑1 (MQ=255)
                        aCGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGc                                                              <  1:811034/54‑1 (MQ=255)
                         cGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCAt                                             <  1:3354042/70‑1 (MQ=255)
                         cGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATc                                            >  1:2674786/1‑71 (MQ=255)
                              gATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                        <  1:3042518/70‑1 (MQ=255)
                               aTAGCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGaa                                      <  1:2580252/71‑1 (MQ=255)
                               aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                      <  1:2808897/71‑1 (MQ=255)
                               aTACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATTCCAGCGCAGTCACGTAATCTTCGTCATAGaa                                      <  1:2234781/71‑1 (MQ=255)
                                 aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGt                                                           >  1:1462464/1‑48 (MQ=255)
                                 aCCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                    >  1:1505129/1‑71 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                    >  1:546528/1‑70 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                    <  1:3797704/70‑1 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                   <  1:3956211/71‑1 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                   <  1:4037924/71‑1 (MQ=255)
                                  ccAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                   <  1:377206/71‑1 (MQ=255)
                                     aTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                    >  1:2628690/1‑67 (MQ=255)
                                     aTTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                    >  1:2290176/1‑67 (MQ=255)
                                      ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCg                                           <  1:1210598/59‑1 (MQ=255)
                                      ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                  <  1:436537/68‑1 (MQ=255)
                                      ttCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                               >  1:358568/1‑71 (MQ=255)
                                       tCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCtt                                                >  1:2958831/1‑53 (MQ=255)
                                       tCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCtt                              >  1:1759730/1‑71 (MQ=255)
                                        cccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg                                                             >  1:2142245/1‑39 (MQ=255)
                                        cccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGAc                                                        >  1:2129808/1‑44 (MQ=255)
                                        cccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                             >  1:2614259/1‑71 (MQ=255)
                                         ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                             <  1:3279410/70‑1 (MQ=255)
                                         ccAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                             <  1:3062278/70‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTa                                                     <  1:1052931/45‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTa                                                     <  1:2728139/45‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTa                                                     <  1:3883741/45‑1 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                           >  1:2620179/1‑71 (MQ=255)
                                          cAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgt                           >  1:3439755/1‑71 (MQ=255)
                                           aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                        <  1:1768908/57‑1 (MQ=255)
                                           aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                          <  1:1029968/71‑1 (MQ=255)
                                           aGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtc                          <  1:3570969/71‑1 (MQ=255)
                                             aCCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGt                                                      <  1:2571890/41‑1 (MQ=255)
                                             aCCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                        <  1:2796468/71‑1 (MQ=255)
                                              ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtcgtcgg                        >  1:2149095/1‑69 (MQ=255)
                                              ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                        >  1:1793098/1‑70 (MQ=255)
                                              ccTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCgtcgtc                       <  1:2487568/71‑1 (MQ=255)
                                               ctgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATc                                                  >  1:70646/1‑43 (MQ=255)
                                               ctgctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCgtcgtc                       <  1:109452/70‑1 (MQ=255)
                                                 gctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAcc                    >  1:2749653/1‑71 (MQ=255)
                                                 gctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCAcc                    >  1:3228973/1‑71 (MQ=255)
                                                  ctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                        >  1:487836/1‑50 (MQ=255)
                                                  ctgTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGaa                                      >  1:2553085/1‑52 (MQ=255)
                                                     tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACa                                    >  1:925193/1‑51 (MQ=255)
                                                     tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                               <  1:1387151/56‑1 (MQ=255)
                                                     tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                               <  1:1998831/56‑1 (MQ=255)
                                                     tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCt                               <  1:224303/56‑1 (MQ=255)
                                                       ggcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc              >  1:3454360/1‑71 (MQ=255)
                                                       ggcggTAAGCCATGTTACAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc              >  1:3760169/1‑71 (MQ=255)
                                                           gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAg                                        <  1:3008451/41‑1 (MQ=255)
                                                           gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                  <  1:1312081/47‑1 (MQ=255)
                                                           gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCg                <  1:3195856/65‑1 (MQ=255)
                                                           gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCg                <  1:1282271/65‑1 (MQ=255)
                                                            tAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATc                                  <  1:516229/46‑1 (MQ=255)
                                                             aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACAt                                   >  1:3995179/1‑44 (MQ=255)
                                                             aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgc        >  1:3583245/1‑71 (MQ=255)
                                                              aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                        >  1:2092251/1‑54 (MQ=255)
                                                              aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTcgtcgt                        >  1:371267/1‑54 (MQ=255)
                                                              aGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg       >  1:311101/1‑71 (MQ=255)
                                                                ccATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc     <  1:2123592/71‑1 (MQ=255)
                                                                 cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTtc                             <  1:996445/46‑1 (MQ=255)
                                                                 cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa    <  1:3625449/71‑1 (MQ=255)
                                                                 cATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa    <  1:152277/71‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg         <  1:1349057/64‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg         <  1:2522804/64‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcg         <  1:994455/64‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCaa   <  1:1067070/70‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCa    <  1:1929935/69‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc  <  1:363002/71‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc  <  1:1771479/71‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc  >  1:1174749/1‑71 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc  <  1:3492382/71‑1 (MQ=255)
                                                                   tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc  <  1:3954621/71‑1 (MQ=255)
                                                                    |                                                                     
GCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAC  >  minE/1932976‑1933113

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: