Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,622,878 T→C 9.4% G393G (GGT→GGC ppk → polyphosphate kinase, component of RNA degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,622,8780TC 9.4% 174.3 / 9.6 64G393G (GGT→GGCppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base C (3/3);  total (30/34)
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.22e-01

AATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGA  >  W3110S.gb/2622809‑2622945
                                                                     |                                                                   
aaTTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGt                                                                     <  1:2195466/70‑1 (MQ=255)
          tACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                           <  1:4861754/70‑1 (MQ=255)
           aCCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGCAAGAAAGTCAc                                                          <  1:496232/70‑1 (MQ=255)
             cgcgTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCAc                                                          <  1:2059632/68‑1 (MQ=255)
             cgcgTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgt                                                       >  1:2005954/1‑71 (MQ=255)
               cgTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtgg                                                     <  1:231293/71‑1 (MQ=255)
                gTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtggt                                                    <  1:4665979/71‑1 (MQ=255)
                 tGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAg                                                                  <  1:490434/56‑1 (MQ=255)
                 tGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                           <  1:5118077/63‑1 (MQ=255)
                 tGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:4715710/68‑1 (MQ=255)
                 tGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:461334/68‑1 (MQ=255)
                    cGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGtt                                                <  1:4455884/71‑1 (MQ=255)
                        aGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGaaa                                                               >  1:2414101/1‑52 (MQ=255)
                        aGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAg                                                              <  1:13771/53‑1 (MQ=255)
                         gATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGtt                                           <  1:312507/71‑1 (MQ=255)
                           ttCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGt                                            >  1:930472/1‑68 (MQ=255)
                            tCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTa                                          >  1:4496536/1‑69 (MQ=255)
                             cACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTggtggt                                                 >  1:2612982/1‑61 (MQ=255)
                                gCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGgcg                                    <  1:4539970/71‑1 (MQ=255)
                                   tcatcGACTCGATGATCCACGCCGCACATAACGGt                                                                     <  1:1570239/35‑1 (MQ=255)
                                   tcatcGACTCGATGATCCACGCCGCACATAACGGTAAGa                                                                 <  1:593511/39‑1 (MQ=255)
                                   tcatcGACTCGATGATCCACGCCGCACATAACGGTAAGa                                                                 <  1:3046351/39‑1 (MQ=255)
                                   tcatcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:3689432/50‑1 (MQ=255)
                                   tcatcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:3518405/50‑1 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTaa                                                                   >  1:2901876/1‑35 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAg                                                              >  1:4232991/1‑40 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGt                                                             >  1:370154/1‑41 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCg                                                        >  1:4611720/1‑46 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCg                                                        >  1:2043623/1‑46 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGttt                               >  1:1668428/1‑71 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGttt                               >  1:1525256/1‑71 (MQ=255)
                                     atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGtt                                <  1:4864666/70‑1 (MQ=255)
                                      tcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:4343761/47‑1 (MQ=255)
                                      tcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:1429259/47‑1 (MQ=255)
                                      tcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTc                              >  1:1926343/1‑71 (MQ=255)
                                       cGACTCGATGATCCACTCCGCACATAACGGTAAGAAAGTCACCgtg                                                      <  1:2758247/46‑1 (MQ=255)
                                       cGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCAc                                                          >  1:3440245/1‑42 (MQ=255)
                                        gACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGtt                                >  1:4790100/1‑67 (MQ=255)
                                        gACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcga                            >  1:2395546/1‑71 (MQ=255)
                                           tCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgt                                                       <  1:1927919/41‑1 (MQ=255)
                                           tCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgt                                                       <  1:2102739/41‑1 (MQ=255)
                                           tCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgt                                                       <  1:1501366/41‑1 (MQ=255)
                                            cGATGATCCACGCCGCACATAACGGTAAGAAAGTCAc                                                          >  1:3664938/1‑37 (MQ=255)
                                            cGATGATCCACGCCGCACATAACGGTAAGAAAGTCAc                                                          >  1:138806/1‑37 (MQ=255)
                                            cGATGATCCACGCCGCACATAACGGTAAGAAAGTCAc                                                          >  1:4847571/1‑37 (MQ=255)
                                             gatgatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTa                                          >  1:820751/1‑52 (MQ=255)
                                                gatCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTAc                                         >  1:3276129/1‑50 (MQ=255)
                                                  tCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGcc                  <  1:1382205/71‑1 (MQ=255)
                                                    cACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCa                 >  1:4945467/1‑70 (MQ=255)
                                                     aCGCCGCACATAACGGTACGAAAGTCACCGTGGTGGTTGAGTTACAGgcg                                    <  1:1032955/50‑1 (MQ=255)
                                                      cgccgcACATAACGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACa              <  1:41086/71‑1 (MQ=255)
                                                          gcACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCa          <  1:2125389/71‑1 (MQ=255)
                                                           cacaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                          <  1:3258426/54‑1 (MQ=255)
                                                            acaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAgg                                      <  1:3974335/41‑1 (MQ=255)
                                                            acaTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAgg                                      <  1:23785/41‑1 (MQ=255)
                                                             caTAACGGTAAGAAAGTCACCGTGGTTGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTg       <  1:1796936/71‑1 (MQ=255)
                                                              aTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGgcg                                    >  1:2716849/1‑41 (MQ=255)
                                                              aTAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGgcg                                    >  1:1154204/1‑41 (MQ=255)
                                                               tAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcga                            >  1:1062260/1‑48 (MQ=255)
                                                               tAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACAtt            <  1:1671134/64‑1 (MQ=255)
                                                               tAACGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                          >  1:4371150/1‑50 (MQ=255)
                                                               tAACGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCa          >  1:1110987/1‑66 (MQ=255)
                                                               tAACGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCa          >  1:502231/1‑66 (MQ=255)
                                                                aaCGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTggg     >  1:4684899/1‑70 (MQ=255)
                                                                aaCGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGc    >  1:612276/1‑71 (MQ=255)
                                                                  cGGCAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGa  <  1:29755/71‑1 (MQ=255)
                                                                     tAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCaa                >  1:2311902/1‑54 (MQ=255)
                                                                     |                                                                   
AATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGA  >  W3110S.gb/2622809‑2622945

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: