Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,225,819 A→T 8.8% intergenic (+335/+37) ycgI → / ← minE hypothetical protein/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,225,8190AT 8.8% 180.5 / 9.4 68intergenic (+335/+37)ycgI/minEhypothetical protein/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base T (4/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

GCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAC  >  W3110S.gb/1225755‑1225888
                                                                |                                                                     
gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                 <  1:1445424/71‑1 (MQ=255)
   aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                >  1:7360424/1‑69 (MQ=255)
   aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                              <  1:6113811/71‑1 (MQ=255)
    gAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                              >  1:1258811/1‑70 (MQ=255)
      aTACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                             >  1:3670320/1‑69 (MQ=255)
       tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                          <  1:1044483/71‑1 (MQ=255)
         cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGa                                                         <  1:2978467/70‑1 (MQ=255)
         cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                        <  1:4589808/71‑1 (MQ=255)
         cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                        <  1:743882/71‑1 (MQ=255)
            gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaaa                                                     >  1:2149230/1‑71 (MQ=255)
             ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAct                                                                     >  1:53215/1‑54 (MQ=255)
             ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAAct                                                                     >  1:638864/1‑54 (MQ=255)
             ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                    <  1:381470/71‑1 (MQ=255)
              cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                   >  1:1889594/1‑71 (MQ=255)
              cgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATg                                                   >  1:68299/1‑71 (MQ=255)
               gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActct                                                                   >  1:3973679/1‑54 (MQ=255)
               gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                  >  1:4348233/1‑71 (MQ=255)
               gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                  >  1:7106393/1‑71 (MQ=255)
                cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                                 >  1:3761247/1‑71 (MQ=255)
                 tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                       >  1:2266672/1‑48 (MQ=255)
                 tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                          <  1:125694/61‑1 (MQ=255)
                 tttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                >  1:4742371/1‑71 (MQ=255)
                  ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                               >  1:7038584/1‑71 (MQ=255)
                  ttACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                               >  1:4558819/1‑71 (MQ=255)
                      cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                           <  1:5500267/71‑1 (MQ=255)
                        caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                               >  1:6238972/1‑65 (MQ=255)
                        caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                           >  1:2550394/1‑69 (MQ=255)
                        caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                           >  1:6065406/1‑69 (MQ=255)
                         aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                        >  1:920132/1‑71 (MQ=255)
                           cGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                      >  1:1653794/1‑71 (MQ=255)
                            gTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                     <  1:6767975/71‑1 (MQ=255)
                            gTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                     <  1:7521284/71‑1 (MQ=255)
                               gATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTGTAAAAGCgcat                                   <  1:2121482/70‑4 (MQ=255)
                                aTGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                     <  1:5288757/67‑1 (MQ=255)
                                   agagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAt                               <  1:827758/70‑1 (MQ=255)
                                    gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                             <  1:5369527/71‑1 (MQ=255)
                                      ggtaTGAAAAACTCAAGGCAGTGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                            <  1:3412021/67‑1 (MQ=255)
                                         aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                       >  1:3203705/1‑56 (MQ=255)
                                            aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAttt                             <  1:613258/63‑1 (MQ=255)
                                            aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGctct                      <  1:5911357/70‑1 (MQ=255)
                                            aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtctt                     <  1:7809269/71‑1 (MQ=255)
                                             aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                    <  1:3069975/71‑1 (MQ=255)
                                               aaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttct                   >  1:3259003/1‑70 (MQ=255)
                                                aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                     >  1:4951468/1‑51 (MQ=255)
                                                aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                  >  1:1444189/1‑70 (MQ=255)
                                                    aaGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                                 <  1:5975/35‑1 (MQ=255)
                                                    aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTGTGCTTcc             >  1:6997066/1‑71 (MQ=255)
                                                    aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc             >  1:5124286/1‑71 (MQ=255)
                                                     aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTAc                                          <  1:273651/41‑1 (MQ=21)
                                                     aGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTAc                                          <  1:2004083/41‑1 (MQ=21)
                                                     aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg            >  1:662957/1‑71 (MQ=255)
                                                       gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAg                                        >  1:5641261/1‑41 (MQ=25)
                                                       gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAg                                        >  1:4130513/1‑41 (MQ=25)
                                                       gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCgg           >  1:4177815/1‑70 (MQ=255)
                                                       gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCgg           >  1:3557465/1‑70 (MQ=255)
                                                       gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                        <  1:2011373/41‑1 (MQ=255)
                                                       gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                        <  1:2119099/41‑1 (MQ=255)
                                                       gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAg                                        <  1:7240866/41‑1 (MQ=255)
                                                       gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt          <  1:237033/71‑1 (MQ=255)
                                                       gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt          <  1:809162/71‑1 (MQ=255)
                                                         agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                              >  1:968508/1‑49 (MQ=255)
                                                         agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTAtt                              >  1:4836855/1‑49 (MQ=255)
                                                         agagATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt          >  1:3179636/1‑69 (MQ=255)
                                                          gagaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg       <  1:361630/71‑1 (MQ=255)
                                                           agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg      >  1:1165348/1‑71 (MQ=255)
                                                             aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGt          <  1:2711735/65‑1 (MQ=255)
                                                             aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGt     <  1:1279513/70‑1 (MQ=255)
                                                              tAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTc    <  1:6095792/70‑1 (MQ=255)
                                                              tAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa   >  1:1228936/1‑71 (MQ=255)
                                                               aaCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                         >  1:1831687/1‑48 (MQ=255)
                                                               aaCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAc  >  1:7786626/1‑71 (MQ=255)
                                                               aaCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAc  >  1:897727/1‑71 (MQ=255)
                                                                |                                                                     
GCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAC  >  W3110S.gb/1225755‑1225888

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: