Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,569,965 T→C 100% Q191Q (CAA→CAG yddW ← predicted liprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,569,9650TC81.6% 88.4 / 12.8 38Q191Q (CAA→CAGyddWpredicted liprotein
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base C (14/17);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 3.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGA  >  W3110S.gb/1569903‑1570017
                                                              |                                                    
aCCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATa                                              <  1:1280574/71‑1 (MQ=255)
aCCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATa                                              <  1:946773/71‑1 (MQ=255)
aCCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATa                                              <  1:357621/71‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                  <  1:1075915/57‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                  <  1:975693/57‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                  <  1:95618/57‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                  <  1:899500/57‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                  <  1:1972898/57‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGca                                                  <  1:1695376/57‑1 (MQ=255)
          tGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCAAGTCGCGGTGTTGca                                                  <  1:467048/57‑1 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCTGTGCTGCACATAGACGCTCGCCgg                                  >  1:437025/1‑71 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCgg                                  >  1:1297116/1‑71 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCgg                                  >  1:1525235/1‑71 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCgg                                  >  1:844573/1‑71 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCgg                                  >  1:532977/1‑71 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCgg                                  >  1:1402280/1‑71 (MQ=255)
            cccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCgg                                  >  1:356917/1‑71 (MQ=255)
             ccGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgt                                 <  1:738173/71‑1 (MQ=255)
                  tCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttgtt                             >  1:995902/1‑70 (MQ=255)
                  tCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttgt                              >  1:861938/1‑69 (MQ=255)
                       gACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttg                               >  1:914720/1‑63 (MQ=255)
                       gACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgttg                               >  1:1459487/1‑63 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:381970/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:1201470/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:1495113/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:1351964/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:1656184/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:1471943/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:49996/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:386759/71‑1 (MQ=255)
                        aCAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACag                      <  1:34807/71‑1 (MQ=255)
                              cGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGgt                                 <  1:1532213/54‑1 (MQ=255)
                                           ttCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTgag                          <  1:893737/48‑1 (MQ=255)
                                           ttCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTgag                          <  1:1034817/48‑1 (MQ=255)
                                             cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa  >  1:726433/1‑70 (MQ=255)
                                             cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa  >  1:52363/1‑70 (MQ=255)
                                             cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa  >  1:299876/1‑70 (MQ=255)
                                             cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa  <  1:1478447/70‑1 (MQ=255)
                                                              |                                                    
ACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGA  >  W3110S.gb/1569903‑1570017

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: