Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 704,333 T→C 6.2% S104S (TCA→TCG) ‡ ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported704,3330TC 6.2% 184.3 / 2.3 64S104S (TCA→TCG) ‡ymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATG  >  Exported/704264‑704399
                                                                     |                                                                  
cGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGATTGATTTGTCTTTGTATg                                                                   <  1:2826470/71‑1 (MQ=255)
    aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg                                                               <  1:298368/71‑1 (MQ=255)
     cTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa                                                              <  1:4845245/71‑1 (MQ=255)
       tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                            >  1:1750856/1‑71 (MQ=255)
       tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                            >  1:3675203/1‑71 (MQ=255)
        gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg                                                           >  1:2796996/1‑71 (MQ=255)
         cGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg                                                           <  1:2704464/70‑1 (MQ=255)
         cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                                          <  1:3651759/71‑1 (MQ=255)
         cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                                          <  1:4096640/71‑1 (MQ=255)
                 aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                                  <  1:2124752/71‑1 (MQ=255)
                 aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                                  <  1:2746868/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:3531913/1‑71 (MQ=255)
                          cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                                         >  1:4357608/1‑71 (MQ=255)
                                acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt                                                                    <  1:2463480/38‑1 (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:5130217/71‑1 (MQ=255)
                                        aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                           <  1:3769852/71‑1 (MQ=255)
                                         gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                                     <  1:4574332/44‑1 (MQ=255)
                                         gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                                     <  1:4611536/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:1142498/71‑1 (MQ=255)
                                                         tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:4540870/71‑1 (MQ=255)
                                                         tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:3776153/71‑1 (MQ=255)
                                                         tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat          <  1:1405437/71‑1 (MQ=255)
                                                             tCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACg      <  1:3245202/71‑1 (MQ=255)
                                                             tCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACg      <  1:3424822/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)
                                                                     |                                                                  
CGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATG  >  Exported/704264‑704399

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: