Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 859,171 A→G 7.3% S105P (TCG→CCG)  ybiU ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb859,1710AG 7.3% 148.9 / 4.3 55S105P (TCG→CCG) ybiUhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

AGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAT  >  W3110S.gb/859105‑859240
                                                                  |                                                                     
agaaATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTg                                                                    <  1:559223/70‑1 (MQ=255)
agaaATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGa                                                                   <  1:3596858/71‑1 (MQ=255)
      tATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAg                                                             >  1:679483/1‑71 (MQ=255)
                gCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCt                                                    >  1:436769/1‑70 (MQ=255)
                  cTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCt                                                    <  1:805431/68‑1 (MQ=255)
                  cTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCt                                                    <  1:3075284/68‑1 (MQ=255)
                  cTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 <  1:4495054/71‑1 (MQ=255)
                  cTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 <  1:3494082/71‑1 (MQ=255)
                   tttGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 >  1:915032/1‑70 (MQ=255)
                    ttGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCtt                                                   >  1:3460696/1‑67 (MQ=255)
                    ttGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTc                                               >  1:1189817/1‑71 (MQ=255)
                       tAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTcgc                                             <  1:4134195/70‑1 (MQ=255)
                       tAGACCTCGTCAAAGAGGTTGCGGTCCAGATAGTCCAGCATCGa                                                                       >  1:3905736/1‑44 (MQ=255)
                            cTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTg                                                  <  1:3485705/60‑1 (MQ=255)
                            cTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGa                                       <  1:316466/71‑1 (MQ=255)
                            cTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGa                                       <  1:3103113/71‑1 (MQ=255)
                            cTCGTCAAAGCGGTTGCAGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGa                                       <  1:4753893/71‑1 (MQ=255)
                             tCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTg                                                                    <  1:2231348/41‑1 (MQ=255)
                              cGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                     >  1:3735167/1‑71 (MQ=255)
                               gTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCtt                                                   >  1:772140/1‑56 (MQ=255)
                               gTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaa                                      <  1:1848033/69‑1 (MQ=255)
                               gTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                                    >  1:997855/1‑71 (MQ=255)
                               gTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                                    <  1:2903477/71‑1 (MQ=255)
                                  aaaGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCa                                                              <  1:669429/42‑1 (MQ=255)
                                  aaaGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTg                                                  >  1:1289816/1‑54 (MQ=255)
                                  aaaGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTg                                                  >  1:2821181/1‑54 (MQ=255)
                                    aGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCtt                                                   <  1:4791281/51‑1 (MQ=255)
                                     gCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggggg                                        >  1:3055435/1‑61 (MQ=255)
                                         ttGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTAt                          <  1:3351670/71‑1 (MQ=255)
                                           gCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTcgcg                                            >  1:1800439/1‑51 (MQ=255)
                                               tCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGa                                       <  1:1386314/52‑1 (MQ=255)
                                               tCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGa                                       <  1:1995131/52‑1 (MQ=255)
                                                ccAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc                       >  1:4179081/1‑67 (MQ=255)
                                                  aGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca                 <  1:2593587/71‑1 (MQ=255)
                                                  aGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca                 <  1:262881/71‑1 (MQ=255)
                                                    aTAGTCCAGCATCGGCTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                     <  1:3785764/49‑1 (MQ=255)
                                                    aTAGTCCAGCATCGGCTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                     <  1:1205055/49‑1 (MQ=255)
                                                    aTAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATc                         <  1:1458081/61‑1 (MQ=255)
                                                      aGTCCAGCATCGGCTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATg                                   >  1:2404956/1‑49 (MQ=255)
                                                      aGTCCAGCATCGGCTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATg                                   >  1:1982477/1‑49 (MQ=255)
                                                       gTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACcgc            <  1:2034273/71‑1 (MQ=255)
                                                       gTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACcgc            <  1:2654175/71‑1 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCCTAGCGCTTGTTcgcg                                            >  1:1602281/1‑36 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCCTAGCGCTTGTTcgcg                                            >  1:3100876/1‑36 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCCTAGCGCTTGTTcgc                                             >  1:634905/1‑35 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCTGGGGAAATGGCCTTTTATCACCGCACAACCGCGAc         >  1:1187159/1‑71 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                     >  1:4262002/1‑43 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaaa                                     >  1:2753890/1‑43 (MQ=255)
                                                          cAGCATCGACTGATCCCAGCATAGCGCTTGTTCGCGGGGGAAATGGCCTTTTAt                          >  1:3761084/1‑54 (MQ=255)
                                                          cAGCATCGACTGATCCCAGACTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGa          >  1:3440287/1‑70 (MQ=255)
                                                           aGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCttt                             <  1:3123951/50‑1 (MQ=255)
                                                               tCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTAt                          <  1:4612047/49‑1 (MQ=255)
                                                               tCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTAt                          <  1:2061182/49‑1 (MQ=255)
                                                                 gACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTaa   >  1:2437527/1‑70 (MQ=255)
                                                                 gACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  >  1:450636/1‑71 (MQ=255)
                                                                 gACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  >  1:1539773/1‑71 (MQ=255)
                                                                  |                                                                     
AGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAT  >  W3110S.gb/859105‑859240

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: