Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,198 C→T 9.5% E415E (GAG→GAA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,1980CT 9.5% 191.2 / 8.5 63E415E (GAG→GAAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (24/33);  new base T (3/3);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGA  >  minE/1933131‑1933264
                                                                   |                                                                  
gcgcCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTc                                                                  <  1:226877/70‑1 (MQ=255)
 cgcCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCa                                                                 <  1:1674580/70‑1 (MQ=255)
   ccTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAg                                                              <  1:2905683/71‑1 (MQ=255)
     tGATCTTCTGCGTCGTTTAATTCGCTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACg                                                            >  1:3551218/1‑71 (MQ=255)
       aTCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTc                                                                  <  1:1361521/63‑1 (MQ=255)
       aTCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTc                                                                  <  1:405186/63‑1 (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: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)
                      cAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAAc                                                                    >  1:1774919/1‑46 (MQ=255)
                      cAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAAc                                                                    >  1:59231/1‑46 (MQ=255)
                         cTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCgg                                                           <  1:3787812/52‑1 (MQ=255)
                         cTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCgg                                                           <  1:1341404/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)
                                                                   cTCAAAGCGGTCTTTGATTTCCGGGTTAACGTCGTTac                               <  1:366091/38‑1 (MQ=255)
                                                                   cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                               <  1:1159705/38‑1 (MQ=255)
                                                                   cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa  <  1:3607278/67‑1 (MQ=255)
                                                                   |                                                                  
GCGCCTGATCTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGA  >  minE/1933131‑1933264

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: