Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,204 G→A 12.1% R413R (CGC→CGT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2040GA12.1% 185.3 / 14.1 66R413R (CGC→CGTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base A (5/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

CCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAG  >  minE/1933134‑1933271
                                                                      |                                                                   
ccTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAg                                                                     <  1:2905683/71‑1 (MQ=255)
  tGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACg                                                                   >  1:3551218/1‑71 (MQ=255)
    aTCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCgg                                                                  >  1:778883/1‑70 (MQ=255)
      cTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCt                                                               >  1:3302330/1‑71 (MQ=255)
        tCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                             >  1:1475227/1‑71 (MQ=255)
         ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCTATGAAGAACTCAAAGCGGTCtg                                                              >  1:315032/1‑69 (MQ=255)
         ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                             >  1:3032118/1‑70 (MQ=255)
         ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                             >  1:3386400/1‑70 (MQ=255)
         ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtg                                                              >  1:3916144/1‑69 (MQ=255)
                aTTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGa                                                           >  1:241460/1‑65 (MQ=255)
                aTTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGa                                                           >  1:559190/1‑65 (MQ=255)
                aTTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCg                                                     >  1:3869722/1‑71 (MQ=255)
                      cTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCgg                                                                  <  1:1341404/52‑1 (MQ=255)
                      cTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCgg                                                                  <  1:3787812/52‑1 (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)
                            aaaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCg                                                                   >  1:2196216/1‑45 (MQ=255)
                             aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtta                                                           >  1:386701/1‑50 (MQ=255)
                             aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGa                                                           >  1:1012146/1‑52 (MQ=255)
                             aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                         <  1:382135/70‑1 (MQ=255)
                             aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                         <  1:310737/70‑1 (MQ=255)
                                ccGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTacg                                     <  1:4146368/71‑1 (MQ=255)
                                 cGTTACCGATATCACGACCACCAATGAAGAACTCAAAg                                                                     <  1:1265064/38‑1 (MQ=255)
                                   ttACCGATTTCACGACCACCAATGAAGAACTCAAAg                                                                     <  1:2501875/36‑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:2875756/70‑1 (MQ=255)
                                            tCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                          <  1:1292465/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:3785117/37‑1 (MQ=255)
                                                                 tCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                      <  1:768298/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)
                                                                      |                                                                   
CCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAG  >  minE/1933134‑1933271

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: