Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,893,179 A→C I110L (ATT→CTT)  ygcP → predicted anti‑terminator regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,893,1790AC98.4% 211.0 / ‑5.2 62I110L (ATT→CTT) ygcPpredicted anti‑terminator regulatory protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (31/30);  total (32/30)
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.81e-01

AGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGA  >  W3110S.gb/2893113‑2893246
                                                                  |                                                                   
aGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCttgt                                                                 >  1:218168/1‑71 (MQ=255)
  tACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCttgtt                                                                >  1:902712/1‑70 (MQ=255)
    ccAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGAt                                                             >  1:2101925/1‑71 (MQ=255)
    ccAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGAt                                                             >  1:2162641/1‑71 (MQ=255)
    ccAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGAt                                                             >  1:2902625/1‑71 (MQ=255)
    ccAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGAt                                                             >  1:447992/1‑71 (MQ=255)
    ccAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGAt                                                             >  1:1622028/1‑71 (MQ=255)
    ccAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGAt                                                             >  1:4608995/1‑71 (MQ=255)
      aaaGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTc                                                           <  1:3929996/71‑1 (MQ=255)
      aaaGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTc                                                           <  1:2590397/71‑1 (MQ=255)
       aaGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCa                                                          <  1:1270504/71‑1 (MQ=255)
         gCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAAt                                                        >  1:3893913/1‑71 (MQ=255)
         gCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAAt                                                        >  1:2671138/1‑71 (MQ=255)
         gCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAAt                                                        >  1:257924/1‑71 (MQ=255)
                tGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGttt                                                 <  1:4210412/71‑1 (MQ=255)
                tGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGttt                                                 >  1:1185900/1‑71 (MQ=255)
                 gTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGttt                                                 >  1:3558421/1‑70 (MQ=255)
                  ttGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCaa                                                         >  1:1321932/1‑61 (MQ=255)
                  ttGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAAtt                                                       >  1:2998189/1‑63 (MQ=255)
                  ttGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTca                                               >  1:2444982/1‑71 (MQ=255)
                   tGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTca                                               >  1:2869320/1‑70 (MQ=255)
                     aaaGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaa                                            <  1:1107148/71‑1 (MQ=255)
                      aaGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaac                                           <  1:649923/71‑1 (MQ=255)
                       aGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaaca                                          >  1:1202466/1‑71 (MQ=255)
                       aGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaaca                                          >  1:3638626/1‑71 (MQ=255)
                       aGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaaca                                          >  1:825339/1‑71 (MQ=255)
                         cGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACAtt                                        <  1:832199/71‑1 (MQ=255)
                         cGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACAtt                                        <  1:1929773/71‑1 (MQ=255)
                         cGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACAtt                                        <  1:2598107/71‑1 (MQ=255)
                         cGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACAtt                                        <  1:1099224/71‑1 (MQ=255)
                            cAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCttgttc                                                               <  1:940002/45‑2 (MQ=255)
                            cAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaa                                            <  1:1980280/64‑1 (MQ=255)
                            cAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaa                                            <  1:1438530/64‑1 (MQ=255)
                                agCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTc                                                           >  1:435460/1‑45 (MQ=255)
                                   ggAAGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAg                              <  1:1588650/71‑1 (MQ=255)
                                      aGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTGAATTTCGttt                                                 >  1:3680222/1‑49 (MQ=255)
                                      aGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGttt                                                 >  1:4473031/1‑49 (MQ=255)
                                      aGGTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTg                           >  1:2765768/1‑71 (MQ=255)
                                       ggTTTTTTCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTgc                          >  1:1625672/1‑71 (MQ=255)
                                          tttttCTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCacaa                                            >  1:1867482/1‑50 (MQ=255)
                                               cTGTATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCg                   <  1:1829868/70‑1 (MQ=255)
                                                  tATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTcaca                                             <  1:396735/41‑1 (MQ=255)
                                                  tATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTcaca                                             <  1:2627017/41‑1 (MQ=255)
                                                  tATTCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGaa                 <  1:2130044/69‑1 (MQ=255)
                                                    ttCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTg                                       <  1:1448616/45‑1 (MQ=255)
                                                    ttCATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTg                                       <  1:637587/45‑1 (MQ=255)
                                                      cATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGa           >  1:2977291/1‑71 (MQ=255)
                                                      cATCGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGa           >  1:2518561/1‑71 (MQ=255)
                                                         cGCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        >  1:2438182/1‑71 (MQ=255)
                                                          gCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        <  1:2978350/70‑1 (MQ=255)
                                                          gCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        <  1:1878601/70‑1 (MQ=255)
                                                          gCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        <  1:1883911/70‑1 (MQ=255)
                                                          gCCTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        >  1:733391/1‑70 (MQ=255)
                                                          gCCTGTTTCTTGTTATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTa      >  1:1893537/1‑71 (MQ=255)
                                                            cTGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        <  1:821607/68‑1 (MQ=255)
                                                             tGTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        <  1:3543812/67‑1 (MQ=255)
                                                              gTTTCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTAtt    >  1:2209795/1‑70 (MQ=255)
                                                               tttCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGa                  <  1:2822559/55‑1 (MQ=255)
                                                               tttCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCg             <  1:2914138/60‑1 (MQ=255)
                                                               tttCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCg             <  1:2276543/60‑1 (MQ=255)
                                                               tttCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg        <  1:400369/65‑1 (MQ=255)
                                                               tttCTTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTga  <  1:3213817/71‑1 (MQ=255)
                                                                  |                                                                   
AGTACCAAAGCCTCAATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGA  >  W3110S.gb/2893113‑2893246

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: