Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,231 G→A 18.6% N404N (AAC→AAT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2310GA18.6% 141.0 / 23.6 59N404N (AAC→AATlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (19/29);  new base A (7/4);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTT  >  minE/1933163‑1933301
                                                                    |                                                                      
aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                                                       <  1:310737/70‑1 (MQ=255)
aaGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                                                       <  1:382135/70‑1 (MQ=255)
   ccGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTacg                                                                   <  1:4146368/71‑1 (MQ=255)
        aCCGATTTCACGACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTa                                                                     >  1:393359/1‑64 (MQ=255)
        aCCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgcg                                                              >  1:1446853/1‑71 (MQ=255)
        aCCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTACGTCGTTACGACgcgc                                                             >  1:2849895/1‑71 (MQ=255)
           gATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgc                                                               <  1:1462510/67‑1 (MQ=255)
           gATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCa                                                           <  1:1090304/71‑1 (MQ=255)
            aTTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAg                                                          >  1:835765/1‑71 (MQ=255)
              ttCACGACCACCAATGAAGACCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                                        >  1:277583/1‑71 (MQ=255)
               tCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                                        <  1:1292465/70‑1 (MQ=255)
               tCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                                        <  1:2875756/70‑1 (MQ=255)
                  cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                    <  1:2515502/71‑1 (MQ=255)
                   gACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                                        >  1:3538469/1‑66 (MQ=255)
                    accaccGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGgggaca                                                  >  1:3168460/1‑69 (MQ=255)
                     ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                 <  1:2176779/71‑1 (MQ=255)
                      caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                <  1:3398231/71‑1 (MQ=255)
                      caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                <  1:188691/71‑1 (MQ=255)
                       accaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                               >  1:1417476/1‑71 (MQ=255)
                          aaTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTg                                            >  1:3490247/1‑71 (MQ=255)
                           aTGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAAcgtcgt                                                                       <  1:1295768/43‑1 (MQ=255)
                            tGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                <  1:517299/65‑1 (MQ=255)
                              aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTc                                              >  1:821439/1‑65 (MQ=255)
                              aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                                         <  1:3500476/70‑1 (MQ=255)
                               agaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                       <  1:900427/71‑1 (MQ=255)
                               agaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGtt                                                                      <  1:791869/40‑1 (MQ=255)
                                 aaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGAcgcg                                                              <  1:478069/46‑1 (MQ=255)
                                 aaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCg                                                        <  1:2243539/52‑1 (MQ=255)
                                  attCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAgggg                                                       >  1:1521751/3‑52 (MQ=255)
                                  attCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCgg                                                       >  1:1986668/3‑52 (MQ=255)
                                   cTCAAAGCGGTCTTTGATTTCCGGGTTAACGTCGTTac                                                                    <  1:366091/38‑1 (MQ=255)
                                   cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                                                    <  1:1159705/38‑1 (MQ=255)
                                   cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                       <  1:3607278/67‑1 (MQ=255)
                                    tCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                                                    <  1:768298/37‑1 (MQ=255)
                                    tCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTac                                                                    <  1:3785117/37‑1 (MQ=255)
                                      aaaaCGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCg                                                        >  1:1037387/1‑47 (MQ=19)
                                      aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGata                                     <  1:2608841/66‑1 (MQ=255)
                                      aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAg                                >  1:592188/1‑71 (MQ=255)
                                      aaaGCGGTCTGTGATTTCCGGGATAACGTCGTTacga                                                                  >  1:3417208/1‑37 (MQ=255)
                                       aaacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGt                               <  1:2930847/68‑1 (MQ=255)
                                           ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                           <  1:3276287/71‑1 (MQ=255)
                                            gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                          <  1:804056/71‑1 (MQ=255)
                                               tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                    <  1:4118789/42‑1 (MQ=255)
                                               tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                       <  1:1018364/55‑1 (MQ=255)
                                               tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAgg                       >  1:946646/1‑71 (MQ=255)
                                                 tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAg                        <  1:1782671/68‑1 (MQ=255)
                                                 tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                     >  1:3032730/1‑71 (MQ=255)
                                                 tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                     >  1:597976/1‑71 (MQ=255)
                                                 tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGacac                                                 <  1:1363723/43‑1 (MQ=16)
                                                 tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGt                      <  1:1036006/70‑1 (MQ=255)
                                                  gATTTCCGGTTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCg                    >  1:2912796/1‑71 (MQ=255)
                                                   aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                      >  1:441733/1‑68 (MQ=255)
                                                   aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCgg                   <  1:694078/71‑1 (MQ=255)
                                                   aTTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACtt                                               >  1:1314136/1‑43 (MQ=16)
                                                    tttCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt                                   <  1:2840527/54‑1 (MQ=255)
                                                      tCCGGGTTAACGTCGTTACGACGCGCCAGCGgaga                                                    <  1:2955099/35‑1 (MQ=255)
                                                      tCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCt                 >  1:812172/1‑70 (MQ=255)
                                                         gggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt             <  1:548156/71‑1 (MQ=255)
                                                             tAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt             >  1:2835017/1‑67 (MQ=255)
                                                             tAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGat                                      >  1:3854578/1‑42 (MQ=16)
                                                              aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCa                                 >  1:2515140/1‑46 (MQ=255)
                                                              aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAt        >  1:1902277/1‑71 (MQ=255)
                                                                cggcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGAtgt      <  1:595333/68‑1 (MQ=255)
                                                                  tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                       <  1:1600906/36‑1 (MQ=255)
                                                                  tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGt                               <  1:1739025/44‑1 (MQ=255)
                                                                  tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                      >  1:2519386/1‑53 (MQ=255)
                                                                    gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGaa                         >  1:2960271/1‑48 (MQ=255)
                                                                    gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa           <  1:114563/62‑1 (MQ=255)
                                                                    gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt  <  1:3072780/71‑1 (MQ=255)
                                                                    gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt  >  1:3825365/1‑71 (MQ=255)
                                                                    |                                                                      
AAGCCGTTACCGATTTCACGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTT  >  minE/1933163‑1933301

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: