Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,677,646 C→A 24.9% G126G (GGG→GGT yphF ← predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,677,6460CA24.9% 48.8 / 6.7 24G126G (GGG→GGTyphFpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base A (4/2);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGAC  >  W3110S.gb/2677579‑2677707
                                                                   |                                                             
cGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGAccccc                                                            <  1:483283/71‑1 (MQ=255)
cGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGAccccc                                                            <  1:2005060/71‑1 (MQ=255)
     cccAGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCtttt                                                        <  1:730347/70‑1 (MQ=255)
       cAGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGAt                                                     >  1:948760/1‑71 (MQ=255)
        aGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCtttt                                                        <  1:2085504/67‑1 (MQ=255)
          ttttttACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTaa                                                  <  1:381231/71‑1 (MQ=255)
                    aaTTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGtgtt                                        >  1:818382/1‑71 (MQ=255)
                        ccAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGa                                                      >  1:3787499/1‑53 (MQ=255)
                                   ccGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCgg                         <  1:3823130/71‑1 (MQ=255)
                                      aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                      >  1:1329517/1‑71 (MQ=255)
                                      aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                      >  1:259059/1‑71 (MQ=255)
                                      aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                      >  1:1817159/1‑71 (MQ=255)
                                      aCCAGATATGCCGAGACATATTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAAt                      >  1:1353477/1‑71 (MQ=255)
                                          gATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAg                                             <  1:1638914/44‑1 (MQ=255)
                                           aTACGCCGAGACATATTTATCGACCCCCTTTTGATTAATAc                                               <  1:2928694/41‑1 (MQ=255)
                                           aTACGCCGAGACATATTTATCGACCCCCTTTTGATTAATAc                                               <  1:3207217/41‑1 (MQ=255)
                                             aCGCCGAGACATATTTATCGACCCCCTTTTGATTAATa                                                >  1:3680779/1‑38 (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:3152600/71‑1 (MQ=255)
                                                       atatTTATCGACACCCTTTTGATTAATACAGGTGTTGTAGCAAATCACTGGAATGCCCGCCTCACTGgcgc     <  1:2385326/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)
                                                                   |                                                             
CGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGACCAGATACGCCGAGACATATTTATCGACCCCCTTTTGATTAATACAGGTGTTGTAGCAAATCACCGGAATGCCCGCTTCACTGGCGCGAC  >  W3110S.gb/2677579‑2677707

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: