Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,099 G→A 15.7% D448D (GAC→GAT) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,0990GA15.7% 222.5 / 28.6 89D448D (GAC→GAT) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (32/43);  new base A (3/11);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGT  >  minE/1933029‑1933169
                                                                      |                                                                      
tCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCg                                                                        <  1:773795/71‑1 (MQ=255)
  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:3943209/71‑1 (MQ=255)
      gTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTc                                                                  <  1:371382/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:3235053/1‑57 (MQ=255)
                                 aaTCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGAt                                                     >  1:2283533/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: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)
                                          catAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                   <  1:261850/48‑1 (MQ=255)
                                          cGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTcgcg                                                               <  1:2219100/38‑1 (MQ=255)
                                           gTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc                                     >  1:1102500/1‑63 (MQ=255)
                                           gTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTgcgc                                     >  1:2320695/1‑63 (MQ=255)
                                            tAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                   <  1:3013944/48‑1 (MQ=12)
                                            tAGAACATGGCTTCGTCGTCACCTGCAGCTTTAGCATTAACCTGTTcc                                                   <  1:3324889/48‑1 (MQ=12)
                                            tAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                                                        <  1:488882/43‑1 (MQ=255)
                                            tAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCt                                                        <  1:489204/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)
                                                                      |                                                                      
TCGGCGGTAAGCCATGTTCCAGTGCGGTGACGTAATCTTCATCGTAGAACATCGCTTCGTCGTCACCTGCGTCTTTCGCGGCAACCTGATCCAGGAAGCGTTGCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGT  >  minE/1933029‑1933169

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: