Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,282 C→T 33.7% P387P (CCG→CCA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2820CT33.7% 97.8 / 75.5 87P387P (CCG→CCAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base T (14/15);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCC  >  minE/1933214‑1933346
                                                                    |                                                                
aTTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTgg                                                                >  1:3352935/1‑71 (MQ=255)
    ccGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTg                                                             >  1:2317851/1‑70 (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:867695/71‑1 (MQ=255)
         ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                                                       <  1:1366373/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)
                         gtgcCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCAct                                       <  1:3935828/69‑1 (MQ=255)
                             cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt                                   >  1:3056262/1‑71 (MQ=255)
                             cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt                                   >  1:1575206/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:2743481/71‑1 (MQ=255)
                               gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                                 <  1:1282024/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:3495522/1‑70 (MQ=255)
                                                ggATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                 >  1:1592644/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)
                                                       ccgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATa               >  1:489336/3‑65 (MQ=255)
                                                       cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                                      >  1:3497586/1‑42 (MQ=255)
                                                       cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                                      >  1:3913153/1‑42 (MQ=255)
                                                       cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                                      >  1:3545151/1‑42 (MQ=255)
                                                       cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGa                 >  1:1451068/1‑63 (MQ=255)
                                                        cgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACttcatc                                  <  1:2671391/44‑1 (MQ=255)
                                                        aGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc                                        <  1:3691115/39‑1 (MQ=255)
                                                          tAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg                                 <  1:3327095/44‑1 (MQ=255)
                                                           aaTAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCaacc     >  1:522311/1‑68 (MQ=255)
                                                             tGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc   <  1:3098484/71‑1 (MQ=255)
                                                              gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc   >  1:2327411/1‑70 (MQ=255)
                                                              gAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAccc   >  1:2167486/1‑70 (MQ=255)
                                                               aaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGtga                    <  1:1770045/52‑3 (MQ=255)
                                                               aaGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATACGTCCCAAcccc  >  1:815524/1‑69 (MQ=255)
                                                               aaGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGAcccc  <  1:3802313/70‑1 (MQ=255)
                                                                aGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                        >  1:1087968/1‑47 (MQ=255)
                                                                aGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCg                        >  1:783357/1‑47 (MQ=255)
                                                                  gTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGAtct                          <  1:2675862/43‑1 (MQ=255)
                                                                    |                                                                
ATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCCC  >  minE/1933214‑1933346

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: