Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,276 G→A 34.1% F389F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2760GA34.1% 92.3 / 71.2 82F389F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (21/33);  new base A (14/14);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCC  >  minE/1933206‑1933345
                                                                      |                                                                     
ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                                                       <  1:147228/71‑1 (MQ=255)
 gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                                                      <  1:1878731/71‑1 (MQ=255)
   cggtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAg                                                                    >  1:2340843/3‑71 (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: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:1739470/40‑1 (MQ=255)
                                               aCTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                                                       <  1:3497932/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:1462838/1‑69 (MQ=255)
                                                 ttCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGAtctctg                      >  1:3427040/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:3167784/70‑1 (MQ=255)
                                                   cTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGt                     >  1:1528363/1‑70 (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)
                                                               ccgTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCACTTCATCAAAGATCTCTGTGACAATa              >  1:489336/3‑65 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                                     >  1:3497586/1‑42 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                                     >  1:3545151/1‑42 (MQ=255)
                                                               cAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActt                                     >  1:3913153/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)
                                                                      |                                                                     
GGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCTCGGTAACGATACGGCCCAGACCC  >  minE/1933206‑1933345

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: