Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,270 A→C 32.6% T391T (ACT→ACG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2700AC32.6% 57.5 / 79.6 80T391T (ACT→ACGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (21/32);  new base C (12/14);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCC  >  minE/1933206‑1933338
                                                                |                                                                    
ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                                                <  1:147228/71‑1 (MQ=255)
 gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                                               <  1:1878731/71‑1 (MQ=255)
 gTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAAtaa                                                               <  1:3803418/71‑3 (MQ=255)
   cggtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAg                                                             >  1:2340843/3‑71 (MQ=255)
     gcgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTaa                                                                  <  1:3846345/62‑1 (MQ=255)
        aTTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTgg                                                        >  1:3352935/1‑71 (MQ=255)
            ccGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTg                                                     >  1:2317851/1‑70 (MQ=255)
               ggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                                           <  1:892255/61‑1 (MQ=255)
               ggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                                                  >  1:3920323/1‑70 (MQ=255)
                gTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                                                  >  1:1595899/1‑69 (MQ=255)
                 ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCg                                                         <  1:3165632/61‑1 (MQ=255)
                 ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                                               <  1:1366373/71‑1 (MQ=255)
                 ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                                               <  1:867695/71‑1 (MQ=255)
                     cgtcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgt                                           <  1:971606/71‑1 (MQ=255)
                      gtcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgtg                                          >  1:3734586/1‑71 (MQ=255)
                       tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGAATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGc                                         <  1:2590603/71‑1 (MQ=255)
                         gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGc                                         >  1:1879522/1‑69 (MQ=255)
                          ttACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGAtg                                            <  1:544849/65‑1 (MQ=255)
                           tACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt                                     >  1:1157883/1‑71 (MQ=255)
                            acggcgTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt                                     <  1:1208896/70‑1 (MQ=255)
                              gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGcc                                  <  1:2321248/71‑1 (MQ=255)
                               acgCGCCATCGGAGAAACTTCTGCCGGATATTCAg                                                                     <  1:405096/35‑1 (MQ=255)
                                 gtgcCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCAct                               <  1:3935828/69‑1 (MQ=255)
                                     cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt                           >  1:1575206/1‑71 (MQ=255)
                                     cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt                           >  1:3056262/1‑71 (MQ=255)
                                     cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt                           <  1:1288964/71‑1 (MQ=255)
                                     cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt                           <  1:2490657/71‑1 (MQ=255)
                                      aGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACttcttc                          <  1:2505579/71‑1 (MQ=255)
                                       gCGGGGACACTTCTGCCGTATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatca                         <  1:1798835/71‑2 (MQ=255)
                                       gCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatca                         <  1:1455906/71‑2 (MQ=255)
                                       gCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatc                          >  1:2304842/1‑70 (MQ=255)
                                       gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                                               <  1:3688637/49‑1 (MQ=255)
                                       gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACttcttc                          >  1:3209953/1‑70 (MQ=255)
                                       gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                         <  1:1282024/71‑1 (MQ=255)
                                       gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                         <  1:2743481/71‑1 (MQ=255)
                                        cGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCActt                              <  1:1325994/65‑1 (MQ=255)
                                        cGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGa                        <  1:3332325/71‑1 (MQ=255)
                                         ggAGAAACTTCTGCCGTATATTCTTTAATGAAGGTGGGCTGAATCAGAtg                                            <  1:1280261/50‑1 (MQ=255)
                                            gaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCa                                 >  1:2147195/1‑58 (MQ=255)
                                             aaaCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTc                                      <  1:3325177/52‑1 (MQ=255)
                                              aaCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGAtct                  <  1:3191140/71‑1 (MQ=255)
                                               aCTTCTGCCGGATATTCCGTAATAAAGGGTGGCTGAATCAGATGTGCTTCTGCCActtca                            <  1:1423155/60‑2 (MQ=255)
                                               aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                                                <  1:3497932/40‑1 (MQ=255)
                                               aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                                                <  1:1739470/40‑1 (MQ=255)
                                                cTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                         <  1:1180128/62‑1 (MQ=255)
                                                cTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                <  1:677601/71‑1 (MQ=255)
                                                 ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt                                     <  1:837524/49‑1 (MQ=38)
                                                 ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatc                          >  1:3214535/1‑60 (MQ=255)
                                                 ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGAtctctg               >  1:3427040/1‑69 (MQ=255)
                                                 ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGAtctctg               >  1:1462838/1‑69 (MQ=255)
                                                 ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGAtctctg               >  1:1793457/1‑69 (MQ=255)
                                                 ttCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCgg               >  1:546764/1‑71 (MQ=255)
                                                  tCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                                                <  1:3571730/37‑1 (MQ=255)
                                                  tCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGa                                              >  1:2275243/1‑39 (MQ=255)
                                                   cTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt              >  1:1528363/1‑70 (MQ=255)
                                                   cTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt              <  1:3167784/70‑1 (MQ=255)
                                                   cTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGt              >  1:324878/1‑70 (MQ=255)
                                                    tGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTaa            >  1:1404275/1‑71 (MQ=255)
                                                     gCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCAc                                >  1:3494584/1‑50 (MQ=255)
                                                     gCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtct                            >  1:1932959/1‑54 (MQ=255)
                                                      ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTa             <  1:2387211/68‑1 (MQ=255)
                                                      ccGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTa             <  1:1810719/68‑1 (MQ=255)
                                                       cGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGAca          <  1:2562765/70‑2 (MQ=255)
                                                       cGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGAc           <  1:3765587/69‑1 (MQ=255)
                                                       cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCt                                     <  1:2843124/43‑1 (MQ=255)
                                                       cGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGACACTTCTTCGAAGATCTCGGTAACGa         >  1:831199/1‑71 (MQ=255)
                                                        ggATATTCCGTAATAAAGGTTGGCTGAATCAGAtg                                            <  1:3517823/35‑1 (MQ=25)
                                                        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGc                                         <  1:1248171/38‑1 (MQ=255)
                                                        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGa                     <  1:2126397/58‑1 (MQ=255)
                                                        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa         >  1:1592644/1‑70 (MQ=255)
                                                        ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa         >  1:3495522/1‑70 (MQ=255)
                                                         gATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATa       <  1:3900767/71‑1 (MQ=255)
                                                           tatTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCt                                     <  1:3785542/39‑1 (MQ=255)
                                                           tatTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACg     >  1:436776/1‑71 (MQ=255)
                                                              tCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGAtctc                 <  1:2358898/56‑1 (MQ=255)
                                                              tCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGcc  >  1:3371050/1‑71 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                              >  1:3545151/1‑42 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                              >  1:3913153/1‑42 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                              >  1:3497586/1‑42 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa         >  1:1451068/1‑63 (MQ=255)
                                                                aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                <  1:3691115/39‑1 (MQ=255)
                                                                |                                                                    
GGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCC  >  minE/1933206‑1933338

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: