Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,210 T→G 15.1% T411T (ACA→ACC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2100TG15.1% 102.2 / 21.5 53T411T (ACA→ACClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (17/28);  new base G (5/3);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGG  >  minE/1933140‑1933280
                                                                      |                                                                      
cTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCt                                                                        >  1:3302330/1‑71 (MQ=255)
  tCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                                      >  1:1475227/1‑71 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCTATGAAGAACTCAAAGCGGTCtg                                                                       >  1:315032/1‑69 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                                      >  1:3386400/1‑70 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                                      >  1:3032118/1‑70 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtg                                                                       >  1:3916144/1‑69 (MQ=255)
          aTTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGa                                                                    >  1:559190/1‑65 (MQ=255)
          aTTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGa                                                                    >  1:241460/1‑65 (MQ=255)
          aTTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCg                                                              >  1:3869722/1‑71 (MQ=255)
                 tCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                                      <  1:3596602/56‑1 (MQ=255)
                 tCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAc                                                       >  1:210115/1‑71 (MQ=255)
                   gCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAc                                                       <  1:3019216/69‑1 (MQ=255)
                    cTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGtt                                                          <  1:2369899/65‑1 (MQ=255)
                    cTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAcatc                                                    <  1:2614188/71‑1 (MQ=255)
                       aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtta                                                                    >  1:386701/1‑50 (MQ=255)
                       aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGa                                                                    >  1:1012146/1‑52 (MQ=255)
                       aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                                  <  1:310737/70‑1 (MQ=255)
                       aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                                  <  1:382135/70‑1 (MQ=255)
                          ccGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTacg                                              <  1:4146368/71‑1 (MQ=255)
                             ttACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtg                                                                     <  1:2148487/45‑1 (MQ=255)
                               aCCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTa                                                >  1:393359/1‑64 (MQ=255)
                               aCCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtg                                                                     <  1:2853619/43‑1 (MQ=255)
                               aCCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgcg                                         >  1:1446853/1‑71 (MQ=255)
                               aCCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTACGTCGTTACGACgcgc                                        >  1:2849895/1‑71 (MQ=255)
                                 cGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtg                                                                     <  1:3422081/41‑1 (MQ=255)
                                  gATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgc                                          <  1:1462510/67‑1 (MQ=255)
                                  gATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCa                                      <  1:1090304/71‑1 (MQ=255)
                                   aTTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAg                                     >  1:835765/1‑71 (MQ=255)
                                     ttCACGACCACCAATGAAGACCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                   >  1:277583/1‑71 (MQ=255)
                                     ttCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGtt                                                          <  1:663721/48‑1 (MQ=255)
                                      tCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                   <  1:1292465/70‑1 (MQ=255)
                                      tCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                   <  1:2875756/70‑1 (MQ=255)
                                         cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                               <  1:2515502/71‑1 (MQ=255)
                                          gACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                   >  1:3538469/1‑66 (MQ=255)
                                           accaccGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGgggaca                             >  1:3168460/1‑69 (MQ=255)
                                            ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                            <  1:2176779/71‑1 (MQ=255)
                                             caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                           <  1:3398231/71‑1 (MQ=255)
                                             caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                           <  1:188691/71‑1 (MQ=255)
                                              accaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                          >  1:1417476/1‑71 (MQ=255)
                                               ccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTaa                                                        >  1:3805603/1‑40 (MQ=255)
                                                caATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTaa                                                        >  1:749294/1‑39 (MQ=255)
                                                 aaTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTg                       >  1:3490247/1‑71 (MQ=255)
                                                  aTGAAGAATTCAAAACGGTCGGTGATTTCCGGGtt                                                          <  1:3956865/35‑1 (MQ=25)
                                                  aTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                                  <  1:1295768/43‑1 (MQ=255)
                                                   tGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                           <  1:517299/65‑1 (MQ=255)
                                                     aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTc                         >  1:821439/1‑65 (MQ=255)
                                                     aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                    <  1:3500476/70‑1 (MQ=255)
                                                      agaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                  <  1:900427/71‑1 (MQ=255)
                                                      agaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGtt                                                 <  1:791869/40‑1 (MQ=255)
                                                        aaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgcg                                         <  1:478069/46‑1 (MQ=255)
                                                        aaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                   <  1:2243539/52‑1 (MQ=255)
                                                         attCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAgggg                                  >  1:1521751/3‑52 (MQ=255)
                                                         attCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCgg                                  >  1:1986668/3‑52 (MQ=255)
                                                          cTCAAAGCGGTCTTTGATTTCCGGGTTAACGTCGTTac                                               <  1:366091/38‑1 (MQ=255)
                                                          cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                               <  1:1159705/38‑1 (MQ=255)
                                                          cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                  <  1:3607278/67‑1 (MQ=255)
                                                           tCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                               <  1:768298/37‑1 (MQ=255)
                                                           tCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                               <  1:3785117/37‑1 (MQ=255)
                                                             aaaaCGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCg                                   >  1:1037387/1‑47 (MQ=19)
                                                             aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGata                <  1:2608841/66‑1 (MQ=255)
                                                             aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAg           >  1:592188/1‑71 (MQ=255)
                                                             aaaGCGGTCTGTGATTTCCGGGATAACGTCGTTacga                                             >  1:3417208/1‑37 (MQ=255)
                                                              aaacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGt          <  1:2930847/68‑1 (MQ=255)
                                                                  ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg      <  1:3276287/71‑1 (MQ=255)
                                                                   gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa     <  1:804056/71‑1 (MQ=255)
                                                                      tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                               <  1:4118789/42‑1 (MQ=255)
                                                                      tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                  <  1:1018364/55‑1 (MQ=255)
                                                                      tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAgg  >  1:946646/1‑71 (MQ=255)
                                                                      |                                                                      
CTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGG  >  minE/1933140‑1933280

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: