Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,411,479 C→T 12.7% G48G (GGC→GGT rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,411,4790CT12.7% 129.8 / 6.4 47G48G (GGC→GGTrnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base T (2/4);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTT  >  W3110S.gb/4411411‑4411544
                                                                    |                                                                 
gAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGTGa                                                                 <  1:3641253/71‑1 (MQ=255)
 aaCATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATTAGCTGGCGGTAGAACTGCACATTGAAGGCg                                                                  <  1:895322/69‑1 (MQ=255)
  aCATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAag                                                               <  1:38678/71‑1 (MQ=255)
  aCATTTAACCAAACGTGAAAAACCGCCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCg                                                                  >  1:3829710/1‑68 (MQ=255)
       tAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGGACTGCACATTGAAGGCGAAGagca                                                           <  1:816628/70‑1 (MQ=255)
           cAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTg                                                      <  1:1036729/71‑1 (MQ=255)
              aCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAg                                                   <  1:68731/71‑1 (MQ=255)
                   aaaaaCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGTGAAGAGCAGCTTGAAGGCCTg                                              >  1:1176186/1‑71 (MQ=255)
                   aaaaaCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTg                                              >  1:3860033/1‑71 (MQ=255)
                   aaaaaCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGAAGCTTGAAGGCCTg                                              >  1:503836/1‑71 (MQ=255)
                     aaaCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGTGAAGAGCAGCTTGAAGGCCTGCg                                            >  1:1180220/1‑71 (MQ=255)
                     aaaCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGaa                                                                >  1:1934104/1‑51 (MQ=255)
                     aaaCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCt                                                        >  1:3942839/1‑59 (MQ=255)
                         cGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgc                                        >  1:632215/1‑71 (MQ=255)
                         cGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgc                                        >  1:2760516/1‑71 (MQ=255)
                            ccAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgccgc                                     <  1:3393024/71‑1 (MQ=255)
                                ccGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGTGa                                                                 <  1:4125506/39‑1 (MQ=255)
                                 cGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGaa                                                    >  1:164906/1‑51 (MQ=255)
                                   tgatgaGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTg                                                      <  1:1932086/47‑1 (MQ=255)
                                     atgaGCTGGCGGTAGAACTGCACATTGAAGGTGAAGAGCAGCTTAAGGCCTGCGTCGCCGCCTgcgcgc                             <  1:2570500/69‑1 (MQ=255)
                                     atgaGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcgcgc                             >  1:3255048/1‑70 (MQ=255)
                                     atgaGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcgcgc                             >  1:3878654/1‑70 (MQ=255)
                                     atgaGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcgcgc                             <  1:3252354/70‑1 (MQ=255)
                                      tgaGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCt                                               <  1:2275864/51‑1 (MQ=255)
                                        aGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGAt                          <  1:1396546/70‑1 (MQ=255)
                                        aGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATg                         <  1:2324779/71‑1 (MQ=255)
                                         gCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATg                         >  1:801007/1‑70 (MQ=255)
                                              cGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                   <  1:443285/71‑1 (MQ=255)
                                              cGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                   <  1:4440353/71‑1 (MQ=255)
                                                gTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGAt                          <  1:2993675/62‑1 (MQ=255)
                                                gTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGAt                          <  1:92006/62‑1 (MQ=255)
                                                gTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                   <  1:1067050/69‑1 (MQ=255)
                                                  aGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAg                                                   <  1:1305725/35‑1 (MQ=255)
                                                  aGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgc                                        >  1:552742/1‑46 (MQ=255)
                                                  aGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgc                                        >  1:963983/1‑46 (MQ=255)
                                                  aGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATg               >  1:1267541/1‑71 (MQ=255)
                                                   gAACTGCACATTGAAGGTGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                   <  1:1771286/66‑1 (MQ=255)
                                                   gAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATgg              >  1:2407797/1‑71 (MQ=255)
                                                           caTTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcgcgc                             <  1:3166358/48‑1 (MQ=255)
                                                           caTTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTgg      >  1:4018700/1‑71 (MQ=255)
                                                           caTTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTgg      >  1:4023998/1‑71 (MQ=255)
                                                             ttGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcgcgc                             <  1:489365/46‑1 (MQ=255)
                                                               gAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATg                         >  1:1636918/1‑48 (MQ=255)
                                                               gAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCtt  >  1:2788506/1‑71 (MQ=255)
                                                               gAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCtt  >  1:1104657/1‑71 (MQ=255)
                                                                aaGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCtt  <  1:1091823/70‑1 (MQ=255)
                                                                   gCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTg       <  1:360606/62‑1 (MQ=255)
                                                                    |                                                                 
GAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTT  >  W3110S.gb/4411411‑4411544

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: