Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,100 T→G 16.7% D448A (GAC→GCC) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1000TG16.7% 166.5 / 40.2 90D448A (GAC→GCC) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (32/43);  new base G (3/12);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTT  >  minE/1933031‑1933170
                                                                     |                                                                      
ggcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGt                                                                        >  1:2236883/1‑70 (MQ=255)
ggcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                                                       <  1:51534/71‑1 (MQ=255)
 gcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTc                                                                       <  1:1933875/70‑1 (MQ=255)
 gcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                                                      <  1:972765/71‑1 (MQ=255)
 gcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                                                      <  1:2166646/71‑1 (MQ=255)
 gcggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                                                      <  1:1826420/71‑1 (MQ=255)
  cggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCtt                                                                     >  1:2751356/1‑71 (MQ=255)
   ggTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCttt                                                                    <  1:1822936/71‑1 (MQ=255)
    gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTc                                                                   <  1:371382/71‑1 (MQ=255)
    gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTc                                                                   <  1:3943209/71‑1 (MQ=255)
      aaGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgc                                                                 >  1:465693/1‑71 (MQ=255)
       gaccATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagca                                                                <  1:3431786/69‑2 (MQ=255)
       gaccATATTCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTagca                                                                <  1:3303413/69‑2 (MQ=255)
           aTGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGc                                                              >  1:1187514/1‑69 (MQ=255)
            tGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAAc                                                           >  1:462003/1‑71 (MQ=255)
              ttCCAGCGCAGTCACGTAATCTTCGTCATAGAACATGGCTTCGTCGTCACATGCAGCTTTagcatt                                                              >  1:3003156/1‑63 (MQ=255)
                ccAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCt                                                                      <  1:2700032/56‑1 (MQ=255)
                ccAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGa                                                       >  1:1081700/1‑71 (MQ=255)
                 cAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGAt                                                      <  1:1971145/71‑1 (MQ=255)
                  aGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                     <  1:746536/71‑1 (MQ=255)
                    tGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATc                                                     >  1:563076/1‑69 (MQ=255)
                      cGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAg                                                  <  1:2388968/70‑1 (MQ=255)
                          gACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGc                                             >  1:429908/1‑71 (MQ=255)
                           aCGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCtt                                                                     <  1:192207/46‑1 (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)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGAt                                                      >  1:2021538/1‑59 (MQ=255)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGAt                                                      >  1:1103000/1‑59 (MQ=255)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGtt                                          <  1:2438942/71‑1 (MQ=255)
                             gTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCAGCAACCTGATCCAGGAAGCGtt                                          >  1:3385580/1‑71 (MQ=255)
                              tAATCTTCATCGTAGAACATCGCTTCTTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                            <  1:2984257/68‑1 (MQ=255)
                              tAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                                                <  1:906911/48‑1 (MQ=255)
                              tAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCg                                            <  1:1175115/68‑1 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGAt                                                      >  1:2283533/1‑57 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGAt                                                      >  1:3235053/1‑57 (MQ=255)
                               aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                        >  1:2933978/1‑71 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcg                                       >  1:3277849/1‑71 (MQ=255)
                                aTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgc                                        <  1:3570258/70‑1 (MQ=255)
                                   ttCGTCATAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                    >  1:3735240/1‑55 (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:3545562/1‑71 (MQ=255)
                                      aTCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGAt                                 >  1:384461/1‑71 (MQ=255)
                                       tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                                <  1:610835/71‑1 (MQ=255)
                                       tCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATc                                >  1:3267113/1‑71 (MQ=255)
                                        catAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                    <  1:261850/48‑1 (MQ=255)
                                        cGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                                                <  1:2219100/38‑1 (MQ=255)
                                         gTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc                                      >  1:2320695/1‑63 (MQ=255)
                                         gTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc                                      >  1:1102500/1‑63 (MQ=255)
                                          tAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                    <  1:3013944/48‑1 (MQ=12)
                                          tAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                    <  1:3324889/48‑1 (MQ=12)
                                          tAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                                                         <  1:489204/43‑1 (MQ=255)
                                          tAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                                                         <  1:488882/43‑1 (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:3332342/70‑1 (MQ=255)
                                              aCATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTCCTGGAAACGTTCAGCCTGATCTTCtgc                          <  1:139504/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)
                                                                     |                                                                      
GGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTT  >  minE/1933031‑1933170

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: