Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,121 A→T 18.1% L441Q (CTG→CAG)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1210AT18.1% 227.2 / 41.7 94L441Q (CTG→CAG) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (34/43);  new base T (11/6);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAAT  >  minE/1933053‑1933191
                                                                    |                                                                      
cGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAg                                                                       <  1:2388968/70‑1 (MQ=255)
    gACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGc                                                                  >  1:429908/1‑71 (MQ=255)
     aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGc                                                                  <  1:2875832/70‑1 (MQ=255)
     aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                 <  1:1330031/71‑1 (MQ=255)
      cGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGt                                                                <  1:2198952/71‑1 (MQ=255)
      cGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                                <  1:3680084/71‑1 (MQ=255)
       gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGtt                                                               <  1:2438942/71‑1 (MQ=255)
       gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCAGCAACCTGATCCAGGAAGCGtt                                                               >  1:3385580/1‑71 (MQ=255)
        tAATCTTCATCGTAGAACATCGCTTCTTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                 <  1:2984257/68‑1 (MQ=255)
        tAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                 <  1:1175115/68‑1 (MQ=255)
         aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                             >  1:2933978/1‑71 (MQ=255)
          aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                                            >  1:3277849/1‑71 (MQ=255)
          aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                                             <  1:3570258/70‑1 (MQ=255)
               cATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGa                                                       >  1:646527/1‑71 (MQ=255)
                aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                 >  1:379081/1‑60 (MQ=255)
                aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGa                                                       <  1:1915100/70‑1 (MQ=255)
                aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAt                                                      >  1:384461/1‑71 (MQ=255)
                aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAt                                                      >  1:3545562/1‑71 (MQ=255)
                 tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                                                     <  1:610835/71‑1 (MQ=255)
                 tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                                                     >  1:3267113/1‑71 (MQ=255)
                   gTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc                                                           >  1:2320695/1‑63 (MQ=255)
                   gTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc                                                           >  1:1102500/1‑63 (MQ=255)
                    tAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTc                                                  >  1:1165781/1‑71 (MQ=255)
                    tAGAACATCGCTTCGTCGTCACCTGCGGCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTc                                                  <  1:3782300/71‑1 (MQ=255)
                     aGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGa                                                       <  1:581029/65‑1 (MQ=255)
                     aGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCTTct                                                 >  1:1747131/1‑70 (MQ=255)
                     aGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTcc                                                 >  1:335065/1‑71 (MQ=255)
                      gAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTg                                                              >  1:1739927/1‑57 (MQ=255)
                       aaCATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCTTCtgc                                               <  1:1624122/71‑1 (MQ=255)
                        aCATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCTTCtgcg                                              <  1:1128480/71‑1 (MQ=255)
                        aCATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCTTCtgc                                               <  1:139504/70‑1 (MQ=255)
                        aCATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCTTCtgc                                               <  1:3332342/70‑1 (MQ=255)
                        aCATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCcgcg                                              <  1:2525895/71‑1 (MQ=255)
                          aTCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCcg                                                <  1:1610090/67‑1 (MQ=255)
                          aTCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCcg                                                <  1:1647274/67‑1 (MQ=255)
                           tCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCcg                                                <  1:1253598/66‑1 (MQ=255)
                            cGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCa                                                                        >  1:2153386/1‑41 (MQ=255)
                              cTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAg                                                                       <  1:34124/40‑1 (MQ=255)
                              cTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGa                                                       <  1:3738862/56‑1 (MQ=255)
                                 cgtcgtCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGt                                                                >  1:1621701/1‑44 (MQ=255)
                                 cgtcgtCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGc                                     >  1:2262898/1‑71 (MQ=255)
                                  gtcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGc                                     <  1:1384504/70‑1 (MQ=255)
                                  gtcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCt                                    <  1:2906310/71‑1 (MQ=255)
                                   tcgtcACCTGCGTCTTTCGCTGCAACCTGATCCAGGAAGCg                                                                 <  1:1332519/41‑1 (MQ=255)
                                   tcgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                 <  1:1229593/41‑1 (MQ=255)
                                    cgtcACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                                                 <  1:3714435/40‑1 (MQ=255)
                                        aCCTGCGTCTTTCTCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                               >  1:3859545/1‑70 (MQ=255)
                                        aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTaa                              <  1:3903421/71‑1 (MQ=255)
                                        aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                               <  1:2483047/70‑1 (MQ=255)
                                        aCCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGAGTCATTCAGCTCGCTa                               <  1:2476857/70‑1 (MQ=255)
                                         ccTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTaaa                             <  1:3922048/71‑1 (MQ=255)
                                          cTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                            >  1:1897727/1‑71 (MQ=255)
                                            tcGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTa                               <  1:2827202/65‑1 (MQ=255)
                                            gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGcc                          >  1:2650001/1‑71 (MQ=255)
                                            gCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGc                           >  1:1698278/1‑70 (MQ=255)
                                              gTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCa                                       >  1:1552691/1‑56 (MQ=255)
                                              gTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGt                        <  1:2317057/71‑1 (MQ=255)
                                               tCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGtt                       <  1:3289036/71‑1 (MQ=255)
                                                cTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGtt                       <  1:2592423/70‑1 (MQ=255)
                                                cTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTa                      >  1:242409/1‑71 (MQ=255)
                                                  ttCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                    >  1:313799/1‑71 (MQ=255)
                                                   tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAcc                    <  1:1573752/70‑1 (MQ=255)
                                                   tCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTAc                     <  1:3486576/69‑1 (MQ=255)
                                                    cgcgGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGa                  >  1:1126789/1‑71 (MQ=255)
                                                     gcgGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGAt                 >  1:2037557/1‑71 (MQ=255)
                                                       ggCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGAtt                <  1:403066/70‑1 (MQ=255)
                                                        ttaaCCTGTTCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGtt                                         >  1:2340006/3‑44 (MQ=12)
                                                        ttaaCCTGTTCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGtt                                         >  1:3865378/3‑44 (MQ=12)
                                                        gCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTc              <  1:2295148/71‑1 (MQ=255)
                                                         taaCCTGTTCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCa             >  1:1499791/2‑71 (MQ=255)
                                                         cAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCAtt                                         <  1:1707248/43‑1 (MQ=255)
                                                         cAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                            <  1:1790523/56‑1 (MQ=255)
                                                         cAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                            <  1:3034267/56‑1 (MQ=255)
                                                          aaCCTGTTCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAAcca                         >  1:1929449/1‑57 (MQ=255)
                                                           aCCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAg                            <  1:2360341/54‑1 (MQ=255)
                                                           aCCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACg           <  1:573730/71‑1 (MQ=255)
                                                            ccTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGa          >  1:3502104/1‑71 (MQ=255)
                                                             cTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGac         <  1:1682021/71‑1 (MQ=255)
                                                              tGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGAttt               <  1:458023/64‑1 (MQ=255)
                                                               gATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGacca       >  1:3153879/1‑71 (MQ=255)
                                                               gATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGacca       >  1:1219652/1‑71 (MQ=255)
                                                                ttCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGaccac      >  1:1321896/2‑71 (MQ=255)
                                                                ttCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGa          >  1:1094129/2‑67 (MQ=255)
                                                                ttCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGa          >  1:1303462/2‑67 (MQ=255)
                                                                aTCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGaccac      >  1:1386299/1‑71 (MQ=255)
                                                                aTCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGacc        <  1:3179727/69‑1 (MQ=255)
                                                                aTCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGacc        <  1:1249779/69‑1 (MQ=255)
                                                                 tCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGaccacc     >  1:433552/1‑71 (MQ=255)
                                                                 tCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGaccacc     >  1:438114/1‑71 (MQ=255)
                                                                 tCCTGGAAACGTTCAGCCTGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGaccac      >  1:1477415/1‑70 (MQ=255)
                                                                 tCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCg                   >  1:1966535/1‑57 (MQ=255)
                                                                 tCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGaccac      >  1:3856698/1‑70 (MQ=255)
                                                                  ccAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGaccacc     >  1:222669/1‑70 (MQ=255)
                                                                  ccAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGacca       >  1:2677720/1‑68 (MQ=255)
                                                                   cAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGa          >  1:2628038/1‑64 (MQ=255)
                                                                    aGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAAt  >  1:1174807/1‑71 (MQ=255)
                                                                    |                                                                      
CGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAAT  >  minE/1933053‑1933191

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