Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,677,682 C→T 15.1% P114P (CCG→CCA yphF ← predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,677,6820CT15.1% 92.0 / 8.3 40P114P (CCG→CCAyphFpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (4/2);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGACACTGCCGATA  >  W3110S.gb/2677614‑2677746
                                                                    |                                                                
ccGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCgg                                                                <  1:3823130/71‑1 (MQ=255)
   aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                                                             >  1:259059/1‑71 (MQ=255)
   aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                                                             >  1:1329517/1‑71 (MQ=255)
   aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                                                             >  1:1353477/1‑71 (MQ=255)
   aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                                                             >  1:1817159/1‑71 (MQ=255)
          aCGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATg                                                            >  1:3558463/1‑65 (MQ=255)
           cGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCtt                                                     <  1:122814/71‑1 (MQ=255)
               gagaCATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACt                                                 <  1:3130190/71‑1 (MQ=255)
                    atatTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAATGCCCGCCTCACTGgcgc                                            <  1:2385326/71‑1 (MQ=255)
                    atatTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAATGCCCGCCTCACTGgcgc                                            <  1:3152600/71‑1 (MQ=255)
                     tatTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGgcgc                                            <  1:305613/70‑1 (MQ=255)
                       tttATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCg                                                        >  1:2434068/1‑56 (MQ=255)
                       tttATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCg                                                        >  1:2238240/1‑56 (MQ=255)
                       tttATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGcgac                                         >  1:2759804/1‑71 (MQ=255)
                                   cTTTTGATTAATACAGGTGTTGTAGCAAATCACCGTAATGCCCGCTTCACTGGCGCGACGAACgg                                   <  1:2446128/65‑1 (MQ=255)
                                   cTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAGTGCCCGCTTCACTGGCGCGACGAACgg                                   <  1:1020862/65‑1 (MQ=255)
                                     tttGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTa                           >  1:2117767/1‑71 (MQ=255)
                                      ttGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCt                            <  1:2969339/69‑1 (MQ=255)
                                       tGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACt                         >  1:891118/1‑71 (MQ=255)
                                          ttAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCt                            <  1:965232/65‑1 (MQ=255)
                                             aTACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGcgacga                                       <  1:3068550/51‑1 (MQ=255)
                                             aTACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGcgacga                                       <  1:599450/51‑1 (MQ=255)
                                             aTACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGcgacga                                       <  1:1202533/51‑1 (MQ=255)
                                               aCAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCAtttt                 >  1:3230340/1‑71 (MQ=255)
                                                cAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACt                                                 >  1:621582/1‑38 (MQ=255)
                                                cAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACt                                                 >  1:1000928/1‑38 (MQ=255)
                                                cAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCAtttt                 >  1:1187565/1‑70 (MQ=255)
                                                    tgttgtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGac            <  1:839260/71‑1 (MQ=255)
                                                      ttgtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCt                            <  1:1004759/53‑1 (MQ=255)
                                                      ttgtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCAtttt                 <  1:21051/64‑1 (MQ=255)
                                                      ttgtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGaca           <  1:3226090/70‑1 (MQ=255)
                                                        gtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGt                                  >  1:1397556/1‑45 (MQ=255)
                                                        gtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTa                                 >  1:99379/1‑46 (MQ=255)
                                                        gtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCa               >  1:3355452/1‑64 (MQ=255)
                                                        gtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCa               >  1:2081229/1‑64 (MQ=255)
                                                        gtAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGACACTg        >  1:436027/1‑71 (MQ=255)
                                                        gtAGCAAATCACAGGAATGCCCGCTTCACTGGCGCGACGAACGGTa                                 >  1:2842479/1‑46 (MQ=255)
                                                           gCAAATGACCGGAATGCCAGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGACACTGcc      >  1:3626/1‑70 (MQ=255)
                                                               aTCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGACACTGCCGATa  >  1:2269192/1‑70 (MQ=255)
                                                                acACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCAt                    >  1:2597642/2‑51 (MQ=255)
                                                                    |                                                                
CCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGACGAACGGTACGGCTACTGCCATTTTCAGACACTGCCGATA  >  W3110S.gb/2677614‑2677746

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