Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,210 T→G 52.6% 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,2100TG52.6% 39.2 / 112.8 99T411T (ACA→ACClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (30/17);  new base G (21/31);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 2.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAAT  >  minE/1933140‑1933275
                                                                      |                                                                 
cTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCt                                                                   >  1:83250/1‑71 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgt                                                                 <  1:148358/70‑1 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCt                                                                   <  1:181014/68‑1 (MQ=255)
   ccGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCt                                                                   <  1:3784642/68‑1 (MQ=255)
    cgcgTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtg                                                                <  1:5053157/70‑1 (MQ=255)
      cgTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGAtt                                                             <  1:1958619/71‑1 (MQ=255)
        tCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGAttt                                                            <  1:2735161/70‑1 (MQ=255)
            tCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCggg                                                       >  1:4314297/1‑71 (MQ=255)
              aGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGt                                                      >  1:5368628/1‑70 (MQ=255)
              aGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGt                                                      >  1:3826387/1‑70 (MQ=255)
               attCGCTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGtt                                                     <  1:4205209/68‑1 (MQ=255)
                ttCGCTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTc                                                           <  1:3011458/62‑1 (MQ=255)
                ttCGCTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTaa                                                   >  1:1234982/2‑71 (MQ=255)
                 tCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAc                                                  >  1:3613736/1‑71 (MQ=255)
                  cGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCg                                                         >  1:2196321/1‑63 (MQ=255)
                   gCTAAAACCATTACCGATTTTACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAcat                                                <  1:4925232/71‑3 (MQ=255)
                   gCTAAAACCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAcat                                                <  1:321381/71‑3 (MQ=255)
                      aaaaCCATTACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAcatcat                                             >  1:2156035/1‑69 (MQ=255)
                      aaaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                             >  1:5640703/1‑71 (MQ=255)
                         gCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                          >  1:2357625/1‑71 (MQ=255)
                         gCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                          >  1:3459425/1‑71 (MQ=255)
                         gCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                          >  1:3650974/1‑71 (MQ=255)
                         gCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                          >  1:1998655/1‑71 (MQ=255)
                           cattACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacgg                                        <  1:217585/69‑2 (MQ=255)
                            gTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTacgac                                       >  1:169675/1‑71 (MQ=255)
                             ttACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAcatcat                                             <  1:4769723/64‑3 (MQ=255)
                             ttACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGAttt                                                            >  1:1390382/1‑49 (MQ=255)
                              tACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGa                                                               >  1:5626423/1‑45 (MQ=255)
                              tACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGa                                                               >  1:4416406/1‑45 (MQ=255)
                              tACCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGa                                                               >  1:3099723/1‑45 (MQ=255)
                              gaCCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGa                                                               >  1:5177432/2‑45 (MQ=255)
                               aCCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgcg                                    >  1:1279968/1‑71 (MQ=255)
                                   aTTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAcatcat                                             >  1:1087368/1‑56 (MQ=255)
                                   aTTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                >  1:610317/1‑71 (MQ=255)
                                   aTTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:4415099/71‑1 (MQ=255)
                                   aTTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:5034452/71‑1 (MQ=255)
                                   aTTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:784231/71‑1 (MQ=255)
                                   aTTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:5182046/71‑1 (MQ=255)
                                   aTTTCACGACCACCAATGAAGAACTCAAAGCGGTCtgtg                                                                <  1:5181708/39‑1 (MQ=255)
                                    tttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacggcgtg                                    <  1:1903804/66‑3 (MQ=255)
                                    tttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacggcgtg                                    <  1:4017068/66‑3 (MQ=255)
                                    tttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacggcgtg                                    <  1:4283486/66‑3 (MQ=255)
                                    tttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGc                               <  1:4345069/71‑1 (MQ=255)
                                    tttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGc                               <  1:114508/71‑1 (MQ=255)
                                    tttCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGc                               <  1:1902831/71‑1 (MQ=255)
                                    tttCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGc                               <  1:2494489/71‑1 (MQ=255)
                                     ttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAc                                                  >  1:1764189/1‑51 (MQ=255)
                                     ttCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAAc                                                  >  1:170271/1‑51 (MQ=255)
                                      tCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGcc                                  >  1:3604451/1‑66 (MQ=255)
                                      tCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCg                              <  1:5217669/70‑1 (MQ=255)
                                      tCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCg                              <  1:2327593/70‑1 (MQ=255)
                                       cACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                              <  1:1954216/69‑1 (MQ=255)
                                        aCGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                             >  1:4191302/1‑53 (MQ=255)
                                         cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGGCAGCGgag                           >  1:4861846/1‑70 (MQ=255)
                                         cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgag                           >  1:2485440/1‑70 (MQ=255)
                                           accaccAATGAAGAACTCTAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                          <  1:5526757/53‑1 (MQ=255)
                                            ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                       >  1:1928819/1‑71 (MQ=255)
                                            ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                       >  1:1245625/1‑71 (MQ=255)
                                              accaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGaa                         <  1:5295457/67‑1 (MQ=255)
                                                 gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTaa                                                   >  1:168162/2‑38 (MQ=25)
                                                 gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTg                  <  1:5022301/70‑1 (MQ=255)
                                                 actGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCt                   <  1:5183383/68‑1 (MQ=255)
                                                 aaTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGcc                                  >  1:3819135/1‑55 (MQ=255)
                                                 aaTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGcc                                  >  1:5497285/1‑55 (MQ=255)
                                                  aTGAAGAATTCAAAACGGTCGGTGATTTCCGGGtt                                                     <  1:1059066/35‑1 (MQ=25)
                                                  aTGAAGAATTCAAAACGGGCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGc                 <  1:162148/71‑1 (MQ=255)
                                                   tGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCt                   <  1:3702050/68‑1 (MQ=255)
                                                    gaagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                       <  1:1030113/63‑1 (MQ=255)
                                                     aagaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCg               <  1:4965131/70‑1 (MQ=255)
                                                     aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAAc                                                  >  1:3454373/1‑35 (MQ=255)
                                                     aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg               >  1:2914009/1‑70 (MQ=255)
                                                      agaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                       <  1:4604216/61‑1 (MQ=255)
                                                       gaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTac                                          <  1:3750045/41‑1 (MQ=19)
                                                       gaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa             >  1:2674907/1‑70 (MQ=255)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacggcg                                      >  1:573003/2‑42 (MQ=16)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacg                                         <  1:5389028/38‑1 (MQ=21)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacg                                         <  1:188270/38‑1 (MQ=21)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTacg                                         <  1:3141335/38‑1 (MQ=21)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:5411983/47‑1 (MQ=19)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:2235573/47‑1 (MQ=19)
                                                          ttCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:1416983/47‑1 (MQ=19)
                                                           tCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCg               >  1:2266241/1‑64 (MQ=255)
                                                           tCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCg               >  1:4740757/1‑64 (MQ=255)
                                                           tCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTc        >  1:5697636/1‑71 (MQ=255)
                                                            cAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt         >  1:4524972/1‑69 (MQ=255)
                                                             aaaaCGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTccg      >  1:4684276/1‑69 (MQ=255)
                                                             aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa             >  1:1058383/1‑64 (MQ=255)
                                                               aacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTa    <  1:5613074/69‑1 (MQ=255)
                                                               aacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTa    <  1:2695751/69‑1 (MQ=255)
                                                               aGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTa    <  1:5057228/71‑1 (MQ=255)
                                                                gCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACgcgc                                   >  1:2978820/1‑39 (MQ=255)
                                                                gCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACgcgc                                   >  1:766922/1‑39 (MQ=255)
                                                                gCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTaa   >  1:2623724/1‑71 (MQ=255)
                                                                 cggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAg                                <  1:3392334/41‑1 (MQ=19)
                                                                 cggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGggga                          >  1:3525716/1‑47 (MQ=19)
                                                                 cGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAAt  <  1:2902075/71‑1 (MQ=255)
                                                                  ggTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGata           >  1:3379239/1‑61 (MQ=255)
                                                                      tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                              <  1:5655235/38‑1 (MQ=255)
                                                                      tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                              <  1:1272112/38‑1 (MQ=255)
                                                                      tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                     <  1:3115275/47‑1 (MQ=255)
                                                                      |                                                                 
CTTCCGCGTCATTCAGCTCGCTAAAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAAT  >  minE/1933140‑1933275

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: