Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,970 G→A 11.9% T491T (ACC→ACT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,9700GA11.9% 223.5 / 16.5 84T491T (ACC→ACTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (43/31);  new base A (4/6);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGC  >  minE/1932900‑1933040
                                                                      |                                                                      
cGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg                                                                        >  1:2432956/1‑71 (MQ=255)
cGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg                                                                        >  1:1225247/1‑71 (MQ=255)
  gggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGt                                                                      <  1:41688/71‑1 (MQ=255)
   ggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTa                                                                     >  1:3629386/1‑71 (MQ=255)
   ggTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTa                                                                     >  1:1318654/1‑71 (MQ=255)
    gTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                                    <  1:2428158/71‑1 (MQ=255)
        tgAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGc                                                                 <  1:540989/70‑1 (MQ=255)
         gAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTg                                                               >  1:2583037/1‑71 (MQ=255)
           gCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATgg                                                                  >  1:1979411/1‑66 (MQ=255)
           gCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATgg                                                                  >  1:950560/1‑66 (MQ=255)
               aaCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGt                                                          >  1:2621741/1‑70 (MQ=255)
                aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                         <  1:2531095/70‑1 (MQ=255)
                 cGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGaa                                                       <  1:2635929/71‑1 (MQ=255)
                   tAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGt                                                              >  1:3772398/1‑62 (MQ=255)
                    aaTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATg                                                                   >  1:653092/1‑56 (MQ=255)
                    aaTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                    >  1:645167/1‑71 (MQ=255)
                     aTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                    >  1:3537682/1‑70 (MQ=255)
                     aTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                    >  1:3127157/1‑70 (MQ=255)
                         ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                               >  1:1149017/1‑71 (MQ=255)
                            ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                                             <  1:2293919/54‑1 (MQ=255)
                            ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAcc                                             <  1:1169249/70‑1 (MQ=255)
                            ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                            >  1:878366/1‑71 (MQ=255)
                            ttttACCGGACGCATCGCCGCGAACAGCATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                            >  1:2314682/1‑71 (MQ=255)
                              ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                                    <  1:2323056/45‑1 (MQ=255)
                              ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTAt                                                                    <  1:2990411/45‑1 (MQ=255)
                              ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                 >  1:3650037/1‑64 (MQ=255)
                               tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                          <  1:1280495/70‑1 (MQ=255)
                               tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                         >  1:2324401/1‑71 (MQ=255)
                                 ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCGGTTGGTGAACAGCATTACCATACg                                        <  1:3049851/70‑1 (MQ=255)
                                   ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                                             >  1:1446317/1‑47 (MQ=255)
                                   ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                     <  1:142934/71‑1 (MQ=255)
                                   ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                     <  1:1995375/71‑1 (MQ=255)
                                   ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                     <  1:3554281/71‑1 (MQ=255)
                                    gACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                     >  1:2499762/1‑70 (MQ=255)
                                       gCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                 >  1:3378994/1‑55 (MQ=255)
                                       gCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                 >  1:1704730/1‑55 (MQ=255)
                                        cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                            <  1:2212840/43‑1 (MQ=255)
                                        cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                            <  1:255271/43‑1 (MQ=255)
                                        cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                            <  1:2511229/43‑1 (MQ=255)
                                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGt                                                              >  1:3600651/1‑40 (MQ=255)
                                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                            >  1:342351/1‑58 (MQ=255)
                                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                >  1:399179/1‑70 (MQ=255)
                                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                >  1:1589153/1‑70 (MQ=255)
                                          tCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt                                                             <  1:998539/40‑1 (MQ=255)
                                           cGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                           >  1:1545448/1‑41 (MQ=255)
                                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                       >  1:2507471/1‑44 (MQ=21)
                                            gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                             >  1:389074/1‑70 (MQ=255)
                                             ccGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATa                                 <  1:2249842/65‑1 (MQ=255)
                                             ccGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAta                           <  1:1387639/71‑2 (MQ=255)
                                             ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt                            >  1:2309855/1‑70 (MQ=255)
                                              cGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAta                           <  1:3545116/70‑2 (MQ=255)
                                               gggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                         <  1:2036599/71‑1 (MQ=255)
                                                ggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                                           <  1:2573551/52‑1 (MQ=255)
                                                ggAACAGAATAACGTCGCGGATGGTATGGCTGTTTGTGat                                                       <  1:113343/40‑2 (MQ=39)
                                                ggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGa                                   <  1:957812/60‑1 (MQ=255)
                                                ggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                         >  1:2420740/1‑70 (MQ=255)
                                                 gAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                       <  1:3836490/71‑1 (MQ=255)
                                                 gAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                         >  1:51047/1‑69 (MQ=255)
                                                 gAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                         >  1:446491/1‑69 (MQ=255)
                                                  aaCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                         >  1:2017540/1‑68 (MQ=255)
                                                   aCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGa                     >  1:2874504/1‑71 (MQ=255)
                                                    gaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                    >  1:2195732/2‑71 (MQ=255)
                                                     aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt                                                <  1:3511196/42‑1 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                               >  1:1159958/1‑41 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                          >  1:2444327/1‑46 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                          >  1:2073444/1‑46 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                 >  1:3221528/1‑71 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                 >  1:729160/1‑71 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                 <  1:2702998/71‑1 (MQ=255)
                                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGAATTACCATa                                          >  1:2595537/1‑46 (MQ=255)
                                                       aaTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAcca                  >  1:439714/1‑69 (MQ=255)
                                                        aTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgc                >  1:265374/1‑71 (MQ=255)
                                                         tAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct               <  1:2236446/71‑1 (MQ=255)
                                                          aaCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgctg              >  1:881339/1‑71 (MQ=255)
                                                           aCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt             <  1:2107835/71‑1 (MQ=255)
                                                           aCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt             <  1:1883498/71‑1 (MQ=255)
                                                             gTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                     >  1:3812308/1‑45 (MQ=255)
                                                             gTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct               <  1:1849638/67‑1 (MQ=255)
                                                              tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                    <  1:1516799/61‑1 (MQ=255)
                                                              tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct               <  1:2579557/66‑1 (MQ=255)
                                                              tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcg           >  1:3709333/1‑70 (MQ=255)
                                                                 cgGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg       <  1:949646/71‑1 (MQ=255)
                                                                  ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct               <  1:596495/62‑1 (MQ=255)
                                                                  ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg       >  1:891360/1‑70 (MQ=255)
                                                                    aTGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt      >  1:622335/1‑69 (MQ=255)
                                                                     tGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAg   <  1:3250687/71‑1 (MQ=255)
                                                                      ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc  >  1:1002513/1‑71 (MQ=255)
                                                                      ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc  >  1:143374/1‑71 (MQ=255)
                                                                      ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc  >  1:2167885/1‑71 (MQ=255)
                                                                      |                                                                      
CGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGC  >  minE/1932900‑1933040

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: