Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,588,402 C→T 52.0% L1032L (CTG→CTA hsdR ← endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,588,4020CT52.0% ‑6.5 / 34.4 25L1032L (CTG→CTAhsdRendonuclease R
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base T (9/4);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

CGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACC  >  W3110S.gb/4588334‑4588471
                                                                    |                                                                     
cGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCa                                                                      <  1:1379227/70‑1 (MQ=255)
   gTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGc                                                                    >  1:1123373/1‑69 (MQ=255)
                       aataaCTGCCTGCAATGCCGGTTGCGCGTTCGGCGAACGTAGCACTAGCGAGTCAAACGCTTCGAGGa                                                 >  1:1570362/1‑68 (MQ=255)
                             tgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcg                                        >  1:264974/1‑71 (MQ=255)
                             tgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTgcg                                        >  1:1745297/1‑71 (MQ=255)
                                       gCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTc                                                      >  1:160357/1‑47 (MQ=255)
                                           gTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                                      <  1:529231/59‑1 (MQ=255)
                                           gTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTgcggc                                      <  1:18462/59‑1 (MQ=255)
                                               cgcgTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCg                       <  1:55158/70‑1 (MQ=255)
                                               cgcgTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCg                       <  1:1340823/70‑1 (MQ=255)
                                               cgcgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCg                       >  1:504522/1‑70 (MQ=255)
                                                gcgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                     <  1:615332/71‑1 (MQ=255)
                                                gcgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                     <  1:286384/71‑1 (MQ=255)
                                                gcgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                     <  1:1653109/71‑1 (MQ=255)
                                                 cgTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTc                                 <  1:887002/58‑1 (MQ=255)
                                                 cgTTCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTa                     >  1:1626001/1‑70 (MQ=255)
                                                   ttCGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGt                               >  1:458559/1‑58 (MQ=255)
                                                    tCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATc                                             >  1:852429/1‑43 (MQ=255)
                                                     cGGGGAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAAt                            >  1:1121062/1‑59 (MQ=255)
                                                      ggggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGa               <  1:186840/71‑1 (MQ=255)
                                                         gAACGTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGAttt            <  1:1089920/71‑1 (MQ=255)
                                                            cGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTac         <  1:1504112/71‑1 (MQ=255)
                                                             gTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTa                              >  1:655744/1‑49 (MQ=255)
                                                             gTTGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAAt                            >  1:1245459/1‑51 (MQ=255)
                                                               tGCACTAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACt      >  1:1391644/1‑71 (MQ=255)
                                                                   ccAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTcacc  <  1:1541033/71‑1 (MQ=255)
                                                                    |                                                                     
CGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACC  >  W3110S.gb/4588334‑4588471

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: