Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,275 T→C 6.6% K56K (AAA→AAG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2750TC 6.6% 208.6 / 5.2 75K56K (AAA→AAGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base C (3/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAAT  >  minE/1934205‑1934343
                                                                      |                                                                    
taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                      <  1:2329692/71‑1 (MQ=255)
taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                      <  1:170121/71‑1 (MQ=255)
taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                      <  1:425590/71‑1 (MQ=255)
 aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                     >  1:3586116/1‑71 (MQ=255)
 aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                     >  1:37839/1‑71 (MQ=255)
 aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                     >  1:3785940/1‑71 (MQ=255)
   taCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTTATCctt                                                                    <  1:2698314/70‑1 (MQ=255)
   taCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                     <  1:190178/69‑1 (MQ=255)
       acgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                     <  1:760319/65‑1 (MQ=255)
       acgGGTCATCATGCGGCCAGCAACGGAGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                     <  1:2024128/65‑1 (MQ=255)
        cgGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                                                                    >  1:3508625/1‑65 (MQ=255)
         gggTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                                                             <  1:3623730/71‑1 (MQ=255)
           gTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                           <  1:1904250/71‑1 (MQ=255)
            tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                          >  1:2079312/1‑71 (MQ=255)
             catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                          >  1:2449958/1‑70 (MQ=255)
              atcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                      <  1:2672363/57‑1 (MQ=255)
                catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCttt                                                                    >  1:1086959/1‑57 (MQ=255)
                catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt                                                       <  1:3895473/70‑1 (MQ=255)
                catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg                                                      >  1:1004369/1‑71 (MQ=255)
                 atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc                                                     >  1:3580720/1‑71 (MQ=255)
                      gCCAGCAACGGCGACTTCTATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa                                                <  1:3632756/71‑1 (MQ=255)
                      gCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                           >  1:430752/1‑60 (MQ=255)
                       ccAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                          <  1:590094/60‑1 (MQ=255)
                        cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAt                                              >  1:797984/1‑71 (MQ=255)
                         aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                         >  1:989762/1‑59 (MQ=255)
                         aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt                                             >  1:56651/1‑71 (MQ=255)
                         aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt                                             <  1:2212748/71‑1 (MQ=255)
                          gCAACGGCGACTTCGATGTTCAGTGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                            <  1:3438780/71‑1 (MQ=255)
                          gCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                            <  1:994345/71‑1 (MQ=255)
                           cAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                           >  1:1360918/1‑55 (MQ=255)
                             aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa                                                <  1:3577816/64‑1 (MQ=255)
                             aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                            <  1:2617371/68‑1 (MQ=255)
                             aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTc                                         <  1:2290633/71‑1 (MQ=255)
                                  gACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGcc                                                                 <  1:238993/42‑1 (MQ=255)
                                   aCTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAgg                                    <  1:1001321/70‑1 (MQ=255)
                                    cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGt                                                               >  1:413024/1‑42 (MQ=255)
                                    cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                                   >  1:3077641/1‑54 (MQ=255)
                                    cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                                   >  1:1893917/1‑54 (MQ=255)
                                    cTTCGATGTTCAGCGCTTACAGTTCTTCGTTCTCTTTGCCg                                                                >  1:2315693/1‑41 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                          <  1:1798553/44‑1 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTg                                                  >  1:1663257/1‑52 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGt                                          >  1:3374137/1‑60 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTAt                                 >  1:2489960/1‑69 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATg                                <  1:19321/70‑1 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                               >  1:874366/1‑71 (MQ=255)
                                       cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                               >  1:2930972/1‑71 (MQ=255)
                                        gATGTTCAGCGCTTCCAGTTCTTCGGTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                               <  1:1010110/70‑1 (MQ=255)
                                         aTGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCaga                                      >  1:1201739/1‑62 (MQ=255)
                                          tGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                           >  1:3409677/1‑40 (MQ=255)
                                            ttCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTc                                                        <  1:2001492/41‑1 (MQ=255)
                                             tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                                                             <  1:810349/35‑1 (MQ=255)
                                             tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATcgc                           >  1:177321/1‑69 (MQ=255)
                                              cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                         >  1:2138374/1‑38 (MQ=255)
                                              cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                        >  1:3699312/1‑71 (MQ=255)
                                              cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                        >  1:3419140/1‑71 (MQ=255)
                                              cAGAGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                         >  1:1711027/1‑38 (MQ=255)
                                               aGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACg                       >  1:2465454/1‑71 (MQ=255)
                                                gcgcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACgg                      <  1:249032/71‑1 (MQ=255)
                                                 cgcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGa                     <  1:640171/71‑1 (MQ=255)
                                                  gcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGaa                    <  1:95399/71‑1 (MQ=255)
                                                  gattCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGaaa                    <  1:3352812/69‑1 (MQ=255)
                                                   cTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTTGTCAGAGGTATGATCGCGACGGaa                    <  1:1306342/70‑1 (MQ=255)
                                                    ttCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGaaaa                   >  1:3717359/1‑70 (MQ=255)
                                                    ttCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAAt                  >  1:1094387/1‑71 (MQ=255)
                                                     tCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATc                 >  1:1628646/1‑71 (MQ=255)
                                                     tCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAAt                  >  1:1441068/1‑70 (MQ=255)
                                                      ccAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAAt                  <  1:914869/69‑1 (MQ=255)
                                                       cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGa                         >  1:135864/1‑61 (MQ=255)
                                                       cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCg                <  1:277073/70‑1 (MQ=255)
                                                       cAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAATca                <  1:2108073/70‑2 (MQ=255)
                                                       cAGTTATTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                            >  1:1832902/1‑42 (MQ=255)
                                                        aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt              >  1:2157153/1‑71 (MQ=255)
                                                        aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt              >  1:1236295/1‑71 (MQ=255)
                                                          ttcttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCg            <  1:3722214/71‑1 (MQ=255)
                                                            cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa                                        <  1:514162/41‑1 (MQ=255)
                                                            cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa                                        <  1:3893408/41‑1 (MQ=255)
                                                            cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg           >  1:1519415/1‑70 (MQ=255)
                                                            cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg           >  1:1925090/1‑70 (MQ=255)
                                                             ttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCagag                                     <  1:1117770/43‑1 (MQ=255)
                                                               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa        <  1:2734720/70‑1 (MQ=255)
                                                               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa        <  1:2298148/70‑1 (MQ=255)
                                                               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGa         <  1:161816/69‑1 (MQ=255)
                                                               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGa         <  1:2677430/69‑1 (MQ=255)
                                                                gTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATcgc                           >  1:1207895/1‑50 (MQ=255)
                                                                gTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAgg      <  1:2219731/71‑1 (MQ=255)
                                                                 ttCTCTTTGCCGTCGAATTCTGCGTGCATTTGGTCAGAGGTATGa                               <  1:2630468/45‑1 (MQ=255)
                                                                 ttCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                               <  1:3478289/45‑1 (MQ=255)
                                                                 ttCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGc     <  1:699005/71‑1 (MQ=255)
                                                                   ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa   >  1:3139375/1‑71 (MQ=255)
                                                                   ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa   >  1:3811096/1‑71 (MQ=255)
                                                                   atctTTGACGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGAt                              >  1:2815884/2‑44 (MQ=255)
                                                                      tttGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa   <  1:2659247/68‑1 (MQ=255)
                                                                      tttGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAAt  <  1:1168111/69‑1 (MQ=255)
                                                                      |                                                                    
TAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAAT  >  minE/1934205‑1934343

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: