Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,252 A→C 33.3% V397V (GTT→GTG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2520AC33.3% 37.5 / 53.1 54V397V (GTT→GTGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (9/27);  new base C (7/11);  total (16/38)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCT  >  minE/1933183‑1933322
                                                                     |                                                                      
accaccAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGaaa                                                                       <  1:489697/71‑1 (MQ=255)
  caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                                      <  1:2885904/70‑1 (MQ=255)
  caccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACt                                                                     <  1:986997/71‑1 (MQ=255)
     caATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCt                                                                  <  1:1541563/71‑1 (MQ=255)
      gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTg                                                                 <  1:990147/70‑1 (MQ=255)
      gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTg                                                                 <  1:2143574/70‑1 (MQ=255)
         gaagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                                                              >  1:767057/1‑71 (MQ=255)
            gaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                                            >  1:3235564/1‑70 (MQ=255)
              attCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGatat                                                         <  1:1320466/69‑1 (MQ=255)
               cTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt                                                        <  1:3411527/71‑1 (MQ=255)
                 cAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCgg                                                             <  1:722270/64‑1 (MQ=255)
                 cAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATAtt                                                        <  1:704346/69‑1 (MQ=255)
                  aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTc                                                       <  1:1405310/69‑1 (MQ=255)
                      cggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGata                                                          >  1:245831/1‑62 (MQ=255)
                       ggTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                                <  1:147228/71‑1 (MQ=255)
                        gTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGa                                               <  1:1878731/71‑1 (MQ=255)
                        gTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAAtaa                                               <  1:3803418/71‑3 (MQ=255)
                         tCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                                                    <  1:2256603/49‑1 (MQ=255)
                          cggtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAg                                             >  1:2340843/3‑71 (MQ=255)
                           tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                            <  1:1264421/55‑1 (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:1288964/71‑1 (MQ=255)
                                                            cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt           >  1:1575206/1‑71 (MQ=255)
                                                            cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt           <  1:2490657/71‑1 (MQ=255)
                                                            cAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtctt           >  1:3056262/1‑71 (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)
                                                                     |                                                                      
ACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCGAAGATCT  >  minE/1933183‑1933322

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: