Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,801 G→A 16.9% I214I (ATC→ATT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,8010GA16.9% 134.0 / 17.7 59I214I (ATC→ATTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (20/29);  new base A (8/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 3.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGA  >  minE/1933731‑1933869
                                                                      |                                                                    
tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAggg                                                                      <  1:182418/71‑1 (MQ=255)
  cGCAGGTACATGTCGAGATCCAGCGCGTTATTGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                                                     <  1:827530/70‑1 (MQ=255)
  cGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                                                     <  1:3943651/70‑1 (MQ=255)
    cAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCa                                                                  <  1:2229378/71‑1 (MQ=255)
     aGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAc                                                                 <  1:1902217/71‑1 (MQ=255)
       gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                                               <  1:2411358/71‑1 (MQ=255)
        tACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTg                                                              >  1:3909127/1‑71 (MQ=255)
           aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcat                                                           >  1:3289528/1‑71 (MQ=255)
           aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGca                                                            >  1:5083439/1‑70 (MQ=255)
            tGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcatc                                                          <  1:2953586/71‑1 (MQ=255)
              tCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcatca                                                         <  1:3846328/70‑1 (MQ=255)
                gagaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                                               >  1:2118789/1‑62 (MQ=255)
                gagaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCg                                                      >  1:1910079/1‑71 (MQ=255)
                   aTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGc                                                    >  1:218407/1‑70 (MQ=255)
                     ccAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAcc                                                                >  1:4806698/1‑56 (MQ=255)
                       aGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTc                                               <  1:299525/71‑1 (MQ=255)
                        gcgcGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTc                                               >  1:1235064/1‑70 (MQ=255)
                        gcgcGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCa                                              <  1:2873253/71‑1 (MQ=255)
                          gcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCaa                                             <  1:4238244/70‑1 (MQ=255)
                           cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt                                           >  1:3994413/1‑71 (MQ=255)
                           cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt                                           <  1:4157284/71‑1 (MQ=255)
                            gTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACtt                                          >  1:4633813/1‑71 (MQ=255)
                             ttATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTc                                         >  1:1266663/1‑71 (MQ=255)
                             ttATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTc                                         >  1:4341644/1‑71 (MQ=255)
                              tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc                                         >  1:2907183/1‑70 (MQ=255)
                              tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc                                         >  1:2900545/1‑70 (MQ=255)
                               atgatgGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCa                                       >  1:5466831/1‑71 (MQ=255)
                               atgatgGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCa                                       >  1:1629431/1‑71 (MQ=255)
                               aTGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCa                                       >  1:2968749/1‑71 (MQ=255)
                                tggtggGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCa                                       <  1:768147/70‑1 (MQ=255)
                                tggtggGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCAt                                      <  1:5063368/71‑1 (MQ=255)
                                 gatgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTACCAt                                      >  1:5275157/3‑70 (MQ=255)
                                  gtgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATaa                                    <  1:2920968/71‑1 (MQ=255)
                                     ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc                                 <  1:2354082/71‑1 (MQ=255)
                                     ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc                                 <  1:1764973/71‑1 (MQ=255)
                                     ggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAg                                  <  1:1906686/70‑1 (MQ=255)
                                      gtaaTAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGcc                                <  1:4942852/68‑1 (MQ=255)
                                        gATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcgc                              >  1:3084545/1‑71 (MQ=255)
                                         aTAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcgcg                             >  1:1694715/1‑71 (MQ=255)
                                          tAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCgg                            <  1:4098648/71‑1 (MQ=255)
                                           aaaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATaaa                                   >  1:810487/1‑63 (MQ=12)
                                           aaaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgcgc                            >  1:980354/1‑69 (MQ=255)
                                            aaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgcgcgc                          >  1:2829337/1‑68 (MQ=255)
                                               ggACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcg                               <  1:1974073/63‑1 (MQ=255)
                                               ggACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAc                       <  1:3778846/71‑1 (MQ=255)
                                                gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAc                       >  1:5549767/1‑70 (MQ=255)
                                                   gcgcAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATg                   <  1:832179/71‑1 (MQ=255)
                                                    cgcAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGttt                                                <  1:1535552/41‑1 (MQ=255)
                                                       aGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACt               <  1:910343/71‑1 (MQ=255)
                                                           gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg           <  1:1211218/71‑1 (MQ=255)
                                                           gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg           <  1:1625086/71‑1 (MQ=255)
                                                           gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg           <  1:5306629/71‑1 (MQ=255)
                                                           gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg           <  1:4079909/71‑1 (MQ=255)
                                                             aCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCAt                    <  1:4141221/60‑1 (MQ=255)
                                                             aCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGaa         >  1:3102022/1‑71 (MQ=255)
                                                             aCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGaa         >  1:4014552/1‑71 (MQ=255)
                                                             aCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGaa         >  1:4413358/1‑71 (MQ=255)
                                                              ccGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCAGCGGTTCACCATg                   <  1:4670574/60‑1 (MQ=255)
                                                                gCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATAc      >  1:5546473/1‑71 (MQ=255)
                                                                  cAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCa    <  1:1764683/71‑1 (MQ=255)
                                                                    gggATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCaga  <  1:452627/71‑1 (MQ=255)
                                                                     ggATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg           <  1:4245741/61‑1 (MQ=255)
                                                                      gATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCa                        >  1:369778/1‑47 (MQ=255)
                                                                      gATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCaga  >  1:3739751/1‑69 (MQ=255)
                                                                      |                                                                    
TACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGA  >  minE/1933731‑1933869

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: