Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,272 C→A 7.0% E57D (GAG→GAT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2720CA 7.0% 276.3 / 3.3 84E57D (GAG→GATlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base A (3/3);  total (42/42)
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.93e-01

TAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAA  >  minE/1934205‑1934342
                                                                   |                                                                      
tgataCGACGGGTCATCATTCGGCCAGCAACCGAGACTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcc                                                                     <  1:2818886/69‑2 (MQ=255)
taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                     <  1:170121/71‑1 (MQ=255)
taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                     <  1:2329692/71‑1 (MQ=255)
taataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTctct                                                                     <  1:425590/71‑1 (MQ=255)
 aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                    >  1:3785940/1‑71 (MQ=255)
 aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                    >  1:37839/1‑71 (MQ=255)
 aataCGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCtt                                                                    >  1:3586116/1‑71 (MQ=255)
  ataCGACGGGTCATCATTCGGCCAGCAACCGAGCCTTCAATGTTTAAGGATTCCAGTTCCTGGTtatcct                                                                    <  1:860711/70‑3 (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)
            tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTtctc                                                                      >  1:509796/1‑58 (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:1893917/1‑54 (MQ=255)
                                    cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                                  >  1:3077641/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:3419140/1‑71 (MQ=255)
                                              cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                       >  1:3699312/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:3811096/1‑71 (MQ=255)
                                                                   ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa  >  1:3139375/1‑71 (MQ=255)
                                                                   |                                                                      
TAATACGACGGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAA  >  minE/1934205‑1934342

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: