Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,255,918 C→T 41.7% G538S (GGC→AGC)  ycgV ← predicted adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,255,9180CT41.7% 10.7 / 25.6 24G538S (GGC→AGC) ycgVpredicted adhesin
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (3/7);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGCACTACTCCCGCTGATATTCAATAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACGTTTAATGTGGCAAA  >  W3110S.gb/1255859‑1255984
                                                           |                                                                  
aGCACTACTCCCGCTGATATTCAATAAATCCCCTTTATTATTAACGCCATTACCCTCGCTAACAACATCAg                                                         <  1:556848/71‑1 (MQ=255)
       ctCCCGCTGATATTCAATAAATCCCCTTTATTATTAACGCCATTACCCTCGCTAACAACATCAGCACGCa                                                   >  1:1713845/1‑70 (MQ=255)
          ccGCTGATATTCAATAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCataa                                                >  1:668830/1‑70 (MQ=255)
                    tCAATAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTg                                     >  1:560504/1‑71 (MQ=255)
                      aaTAAATCCCCTTTATTATTAACGCCATTACCCTCGCTAACAACATCAGCACGCATAATAAAGGTACTGt                                    <  1:795790/70‑1 (MQ=255)
                       aTAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACAt                                                            <  1:141593/45‑1 (MQ=255)
                        tAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTAc                                 >  1:1643507/1‑71 (MQ=255)
                        tAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTAc                                 <  1:324269/71‑1 (MQ=255)
                        tAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATAAGCACGCATAATAAAGGTACTGTTAc                                 >  1:1156994/1‑71 (MQ=255)
                                cTTTATTATTAACGCCATTACCCTCGCTAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAgg                         <  1:1635134/71‑1 (MQ=255)
                                cTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCataa                                                <  1:395724/48‑1 (MQ=255)
                                       attaACGCCATTACCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTctct                   <  1:840020/70‑1 (MQ=255)
                                         taACGCCATTACCCTCGCTAACAACATCAGCACGCATAATAAAGGTACTg                                     >  1:1576491/1‑50 (MQ=255)
                                             gCCATTACCCTCGCCAACAACATCAGCACGCataat                                               >  1:358547/1‑36 (MQ=255)
                                             gCCATTACCCTCGCCAACAACATCAGCACGCataat                                               >  1:1741845/1‑36 (MQ=255)
                                             gCCATTACCCTCGCCAACAACATCAGCACGCataat                                               >  1:1617955/1‑36 (MQ=255)
                                                aTTACCCTCGCTAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACGTTTAAtg         <  1:1223456/71‑1 (MQ=255)
                                                 ttACCCTCGCTAACAACATCAGCACGCATAATAAAGGTa                                        <  1:456287/39‑1 (MQ=255)
                                                   aCCCTCGCTAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACGTTTAAtgtg       >  1:1387337/1‑70 (MQ=255)
                                                   aCCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACg                <  1:366408/61‑1 (MQ=255)
                                                   aCCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACg                <  1:1209761/61‑1 (MQ=255)
                                                       tCGCTAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAgg                         <  1:1749047/48‑1 (MQ=255)
                                                       tCGCTAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAgg                         <  1:416485/48‑1 (MQ=255)
                                                         gCCAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACGTTTAATGTGGCaaa  <  1:48876/69‑1 (MQ=255)
                                                           |                                                                  
AGCACTACTCCCGCTGATATTCAATAAATCCCCTTTATTATTAACGCCATTACCCTCGCCAACAACATCAGCACGCATAATAAAGGTACTGTTACCGCTCAGGTTCTCTACGTTTAATGTGGCAAA  >  W3110S.gb/1255859‑1255984

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: