Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 704,335 A→T 6.3% S104T (TCA→ACA) ‡ ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE704,3350AT 6.3% 177.8 / 3.6 63S104T (TCA→ACA) ‡ymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATG  >  minE/704268‑704399
                                                                   |                                                                
aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg                                                               <  1:298368/71‑1 (MQ=255)
 cTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa                                                              <  1:4845245/71‑1 (MQ=255)
   tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                            >  1:3675203/1‑71 (MQ=255)
   tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                            >  1:1750856/1‑71 (MQ=255)
    gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg                                                           >  1:2796996/1‑71 (MQ=255)
     cGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg                                                           <  1:2704464/70‑1 (MQ=255)
     cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                                          <  1:4096640/71‑1 (MQ=255)
     cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                                          <  1:3651759/71‑1 (MQ=255)
             aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                                  <  1:2746868/71‑1 (MQ=255)
             aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                                  <  1:2124752/71‑1 (MQ=255)
              gTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                                 <  1:3557014/71‑1 (MQ=255)
              gTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCAGAAAt                                                 <  1:2269246/71‑1 (MQ=255)
               tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                                   >  1:1528892/1‑68 (MQ=255)
               tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                                   >  1:130526/1‑68 (MQ=255)
                aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACt                                                                >  1:3500166/1‑54 (MQ=255)
                aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACt                                                                >  1:1236879/1‑54 (MQ=255)
                aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt                                               >  1:1163262/1‑71 (MQ=255)
                 aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt                                              >  1:1807894/1‑71 (MQ=255)
                 aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt                                              >  1:3100519/1‑71 (MQ=255)
                 aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt                                               >  1:4132175/1‑70 (MQ=255)
                   cttcttCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTcc                                            >  1:4450791/1‑71 (MQ=255)
                    ttcttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                                 >  1:3540141/1‑65 (MQ=255)
                    ttcttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCAGCTGAAATCTTcc                                            >  1:2706530/1‑70 (MQ=255)
                      cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                                         >  1:4357608/1‑71 (MQ=255)
                      cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                                         >  1:3531913/1‑71 (MQ=255)
                            acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                                   <  1:4293144/55‑1 (MQ=255)
                            acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCta                                   >  1:428365/1‑71 (MQ=255)
                             caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGGTGAAATCTTCCCACTGGc                                     <  1:4277656/68‑1 (MQ=255)
                             caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                                 >  1:875528/1‑56 (MQ=255)
                             caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCtat                                  <  1:1734292/71‑1 (MQ=255)
                                 cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt                                                      >  1:1814948/1‑47 (MQ=255)
                                 cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt                                                      >  1:2227201/1‑47 (MQ=255)
                                  gAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa                                                              >  1:3783593/1‑38 (MQ=255)
                                   aaGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                                 <  1:5651640/50‑1 (MQ=255)
                                    aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                                         >  1:4380159/1‑57 (MQ=255)
                                    aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                           <  1:3769852/71‑1 (MQ=255)
                                    aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                           <  1:5130217/71‑1 (MQ=255)
                                     gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                                     <  1:4611536/44‑1 (MQ=255)
                                     gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                                     <  1:4574332/44‑1 (MQ=255)
                                      tttAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGc                            <  1:3889830/68‑1 (MQ=255)
                                        tAAAATCGAATGATTTGTCTTTGTATGACTGATcgcc                                                         >  1:1557522/1‑37 (MQ=255)
                                        tAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCa                                          >  1:2202457/1‑52 (MQ=255)
                                        tAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCgtgt                        <  1:2930932/70‑1 (MQ=255)
                                          aaaTCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTcc                                            >  1:3072299/1‑48 (MQ=255)
                                          aaaTCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTgcg                     <  1:4934951/71‑1 (MQ=255)
                                             tCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGcac                  <  1:3229673/71‑1 (MQ=255)
                                               gAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGcaca                 <  1:4553275/70‑1 (MQ=255)
                                                 aTGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt                                               >  1:2240798/1‑38 (MQ=255)
                                                  tGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg             >  1:5656822/1‑71 (MQ=255)
                                                  tGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg             >  1:3307530/1‑71 (MQ=255)
                                                   gATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg             >  1:4534918/1‑70 (MQ=255)
                                                    aTTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGta           >  1:877800/1‑71 (MQ=255)
                                                    aTTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGta           >  1:3034772/1‑71 (MQ=255)
                                                     tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:4540870/71‑1 (MQ=255)
                                                     tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:1405437/71‑1 (MQ=255)
                                                     tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:3776153/71‑1 (MQ=255)
                                                     tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:1142498/71‑1 (MQ=255)
                                                         tCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACg      <  1:3424822/71‑1 (MQ=255)
                                                         tCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACg      <  1:3245202/71‑1 (MQ=255)
                                                             tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGaa    <  1:4487875/69‑1 (MQ=255)
                                                             tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGa     <  1:5520817/68‑1 (MQ=255)
                                                             tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGa     <  1:843602/68‑1 (MQ=255)
                                                             tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATg  >  1:184676/1‑71 (MQ=255)
                                                                   aCTGATGGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                           >  1:715/1‑40 (MQ=255)
                                                                   aCTGATCGCCGCTGAAATCTTCCCACTGGCTATAttt                              >  1:996110/1‑37 (MQ=255)
                                                                   |                                                                
ACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATG  >  minE/704268‑704399

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: