Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,494,403 G→A 42.6% L115L (CTG→TTG)  yjgR ← predicted ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,494,4030GA42.6% 21.2 / 56.9 47L115L (CTG→TTG) yjgRpredicted ATPase
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base A (8/12);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

ATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCA  >  W3110S.gb/4494339‑4494473
                                                                |                                                                      
atatTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCg                                                                   <  1:380451/70‑1 (MQ=255)
atatTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCg                                                                   <  1:3620452/70‑1 (MQ=255)
    tCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAAc                                                               <  1:4446016/70‑1 (MQ=255)
     cAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCg                                                             >  1:4000559/1‑71 (MQ=255)
      aGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGt                                                            >  1:3669465/1‑71 (MQ=255)
        cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTc                                                           >  1:736585/1‑70 (MQ=255)
        cACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCg                                                          >  1:3195794/1‑71 (MQ=255)
         aCGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGc                                                         <  1:3511209/71‑1 (MQ=255)
           gCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGa                                                                  >  1:5129760/1‑60 (MQ=255)
           gCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGa                                                                  >  1:481690/1‑60 (MQ=255)
             cAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCg                                                      >  1:2070293/1‑70 (MQ=255)
                aTTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAg                                                                        >  1:3166038/1‑49 (MQ=255)
                 ttGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCg                                                 <  1:4744095/71‑1 (MQ=255)
                       aTCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGc                                           <  1:1029902/71‑1 (MQ=255)
                       aTCGTTGAGATTCAACAGTAGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGc                                           <  1:3821571/71‑1 (MQ=255)
                          gTTGAGATTCAACAGTCGTGCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCttt                                       <  1:1694043/71‑1 (MQ=255)
                           ttGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCg                                                 <  1:1190188/61‑1 (MQ=255)
                            tGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGa                                               >  1:3023886/1‑62 (MQ=255)
                             gagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTtc                                      <  1:4643915/70‑1 (MQ=255)
                              agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAAc                                                               >  1:4756593/1‑44 (MQ=255)
                              agaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAAc                                                               >  1:4401117/1‑44 (MQ=255)
                                aTTCAACAGTCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCt                                         <  1:2246007/64‑1 (MQ=255)
                                   cAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCa                                <  1:1941851/70‑1 (MQ=255)
                                    aaCAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCaaa                              <  1:1185980/71‑1 (MQ=255)
                                     aCAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAg                             <  1:4591706/71‑1 (MQ=255)
                                         tCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGc                                                         <  1:4475361/39‑1 (MQ=255)
                                         tCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCa                                                    >  1:586089/1‑44 (MQ=255)
                                         tCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGatat                         <  1:1006398/71‑1 (MQ=255)
                                         tCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGatat                         >  1:4816913/1‑71 (MQ=255)
                                         tCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGatat                         >  1:1179730/1‑71 (MQ=255)
                                         tCGTGCCAGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGata                          >  1:2409994/1‑70 (MQ=255)
                                           gTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGccc                                                       >  1:4296609/1‑39 (MQ=255)
                                           gTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGAt                                              >  1:1316085/1‑48 (MQ=255)
                                                aGCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCacac                  >  1:5067202/1‑71 (MQ=255)
                                                 gCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCAcacc                 <  1:1030031/71‑1 (MQ=255)
                                                 gCAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGa                            <  1:3479258/60‑1 (MQ=255)
                                                  cAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAg                             <  1:3213516/58‑1 (MQ=255)
                                                  cAACAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATcc                       >  1:1165759/1‑64 (MQ=255)
                                                   aaCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGa                                               >  1:2396783/1‑39 (MQ=255)
                                                   aaCAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGata                          >  1:77903/1‑60 (MQ=255)
                                                   aaCAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGata                          >  1:1407572/1‑60 (MQ=255)
                                                    aCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGcc                                 <  1:2591561/52‑1 (MQ=255)
                                                    aCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCAcaccac               >  1:2291936/1‑70 (MQ=255)
                                                        cGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCg          <  1:4838539/71‑1 (MQ=255)
                                                         ggCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATccc                      >  1:4464597/1‑58 (MQ=255)
                                                             ccAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCca                     <  1:8320/55‑1 (MQ=255)
                                                             ccAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCg          <  1:713564/66‑1 (MQ=255)
                                                             ccAAATCCGAAACCGTCGCCCGCACCGGATGACATTTCTCGCCAAAGATATCCCACACCACCACCg          <  1:625665/66‑1 (MQ=255)
                                                                gATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCca                     >  1:4406663/1‑52 (MQ=255)
                                                                gATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCa  >  1:4883319/1‑71 (MQ=255)
                                                                |                                                                      
ATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCA  >  W3110S.gb/4494339‑4494473

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: