Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,927 G→C 9.1% L172L (CTC→CTG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9270GC 9.1% 159.1 / 6.1 55L172L (CTC→CTGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base C (3/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GCGAATACCAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGT  >  minE/1933858‑1933996
                                                                     |                                                                     
gCGAATACCAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCg                                                                       <  1:2385275/70‑1 (MQ=255)
    aTACCAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATa                                                                  <  1:1754711/71‑1 (MQ=255)
     tACCAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATaa                                                                 <  1:2259530/71‑1 (MQ=255)
     tACCAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATaa                                                                 <  1:1600095/71‑1 (MQ=255)
       ccAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgaga                                                                    >  1:1411995/1‑66 (MQ=255)
       ccAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACg                                                               >  1:460010/1‑71 (MQ=255)
              ggATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGa                                                        <  1:3847390/71‑1 (MQ=255)
                 tCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCga                                                                      <  1:2194199/54‑1 (MQ=255)
                   tGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgaga                                                                    >  1:3227490/1‑54 (MQ=255)
                   tGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGa                                                           <  1:719903/63‑1 (MQ=255)
                    gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgc                                                  >  1:4098853/1‑71 (MQ=255)
                    gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgc                                                  >  1:2401052/1‑71 (MQ=255)
                    gCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcg                                                   >  1:2299453/1‑70 (MQ=255)
                           cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGAc                                                          <  1:3576545/56‑1 (MQ=255)
                           cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                           <  1:1039876/71‑1 (MQ=255)
                              tttAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATaa                                                                 <  1:3074317/46‑1 (MQ=255)
                              tttAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcgcgc                                                >  1:1055040/1‑63 (MQ=255)
                               ttAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGAt                                                                   <  1:751647/43‑1 (MQ=255)
                                tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgag                                                                     >  1:1707386/1‑40 (MQ=255)
                                tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgag                                                                     >  1:1172067/1‑40 (MQ=255)
                                tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGa                                                           >  1:2789213/1‑50 (MQ=255)
                                tAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                           >  1:656168/1‑66 (MQ=255)
                                 aaaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgag                                                                     >  1:1351284/1‑39 (MQ=255)
                                 aaaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgag                                                                     >  1:3440667/1‑39 (MQ=255)
                                  aaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCgag                                                                     >  1:2322251/1‑38 (MQ=255)
                                   aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                   <  1:3104870/71‑1 (MQ=255)
                                   aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                   <  1:3453558/71‑1 (MQ=255)
                                   aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGc                                   <  1:735565/71‑1 (MQ=255)
                                        ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACg                                                         <  1:1639053/44‑1 (MQ=255)
                                        ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACg                                                         <  1:860012/44‑1 (MQ=255)
                                        ttGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                           >  1:2727446/1‑58 (MQ=255)
                                          gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGAt                                                       <  1:1261431/44‑1 (MQ=255)
                                          gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAgcg                                                   >  1:2007165/1‑48 (MQ=255)
                                             ggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTgg                          <  1:3509917/70‑1 (MQ=255)
                                              gATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGaa                        <  1:1842658/71‑1 (MQ=255)
                                               aTTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAt                       >  1:352920/1‑71 (MQ=255)
                                                tttaTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAtt                      >  1:395179/4‑71 (MQ=255)
                                                ttCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                           >  1:3150550/1‑50 (MQ=255)
                                                 tCATCGTTCGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAtt                      >  1:3280268/1‑70 (MQ=255)
                                                  taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                     <  1:1860825/69‑1 (MQ=255)
                                                  taTCGTTAGCGATGAGGTCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAAttt                     <  1:2011318/69‑1 (MQ=255)
                                                  cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                     >  1:3001961/1‑70 (MQ=255)
                                                  cATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAAttt                     <  1:2690654/70‑1 (MQ=255)
                                                   aTCATTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCt                                           >  1:2223781/1‑47 (MQ=255)
                                                     cGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATcc                 >  1:571348/1‑71 (MQ=255)
                                                       ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGa                         <  1:878876/61‑1 (MQ=255)
                                                       ttGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCgg               <  1:3810504/71‑1 (MQ=255)
                                                         ggAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGa                                         >  1:307841/1‑43 (MQ=255)
                                                          gagatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATc                                       <  1:1448303/44‑1 (MQ=255)
                                                           agatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGc           >  1:2199005/1‑71 (MQ=255)
                                                            gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg          <  1:2532192/71‑1 (MQ=255)
                                                            gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg          <  1:616547/71‑1 (MQ=255)
                                                              tgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa        >  1:3823945/1‑71 (MQ=255)
                                                               gaggtCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCa                                  >  1:370217/1‑44 (MQ=255)
                                                               gaggtCCAGATAACGTTGACGATAACGGACTTCCTGATCCTGCa                                  >  1:365271/1‑44 (MQ=255)
                                                               gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa        >  1:3702366/1‑70 (MQ=255)
                                                               gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc       >  1:3210298/1‑71 (MQ=255)
                                                                agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTa                    <  1:3750292/57‑1 (MQ=255)
                                                                agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc       <  1:385157/70‑1 (MQ=255)
                                                                agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc       <  1:1913520/70‑1 (MQ=255)
                                                                  atCGAGATAACGCTGACGATTGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGc              <  1:656722/61‑1 (MQ=255)
                                                                  atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACg      >  1:3037702/1‑69 (MQ=255)
                                                                    cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt  >  1:1552525/1‑71 (MQ=255)
                                                                    cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt  >  1:2960645/1‑71 (MQ=255)
                                                                     |                                                                     
GCGAATACCAGAGAGGATCTGCGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGT  >  minE/1933858‑1933996

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: