Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,516,711 C→T 17.2% A79V (GCG→GTG)  ydcV → predicted spermidine/putrescine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,516,7110CT17.2% 69.3 / 7.8 29A79V (GCG→GTG) ydcVpredicted spermidine/putrescine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (2/3);  total (12/17)
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.86e-01

GCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTC  >  W3110S.gb/1516641‑1516780
                                                                      |                                                                     
gcaCAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTgc                                                                       <  1:1906252/71‑1 (MQ=255)
  acaGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTgtgc                                                                     >  1:419479/1‑71 (MQ=255)
   caGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCt                                                                    <  1:1071891/71‑1 (MQ=255)
          gTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGtt                                                              <  1:430438/70‑1 (MQ=255)
          gTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGtt                                                              <  1:543415/70‑1 (MQ=255)
          gTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTa                                                             <  1:1629080/71‑1 (MQ=255)
             atatTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAgg                                                           <  1:1453457/70‑1 (MQ=255)
             atatTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAgg                                                           <  1:820439/70‑1 (MQ=255)
                  cTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCt                                                     >  1:1811435/1‑71 (MQ=255)
                   ttGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGTGCTGGTGTTAGGGACGCTg                                                    >  1:699098/1‑71 (MQ=255)
                   ttGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGa                                                         >  1:2582190/1‑66 (MQ=255)
                       tGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAg                                                <  1:3009589/71‑1 (MQ=255)
                           gTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGTGCTGGTGTTAg                                                            >  1:205200/1‑55 (MQ=255)
                              acacTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCcgc                                         >  1:2203147/1‑71 (MQ=255)
                                   gTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTgg                                                                  <  1:3643867/41‑1 (MQ=255)
                                    tCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCtgtg                                   >  1:4247292/1‑71 (MQ=255)
                                    tCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCtgt                                    <  1:3050661/70‑1 (MQ=255)
                                    tCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCtg                                     <  1:1405450/69‑1 (MQ=255)
                                     cACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGtt                                                              >  1:2951884/1‑43 (MQ=255)
                                     cACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGtt                                                              >  1:3904905/1‑43 (MQ=255)
                                       cTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAgg                                                           <  1:634340/44‑1 (MQ=255)
                                        ttAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGcga                               >  1:2216273/1‑71 (MQ=255)
                                          aaaGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGcgac                              >  1:2291734/1‑70 (MQ=255)
                                                 cggcgCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACttt                      <  1:2062134/71‑1 (MQ=255)
                                                   gcgcTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACttttt                    <  1:911651/71‑1 (MQ=255)
                                                       tGGCGACATTAATTGTGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCt                                      <  1:1846531/49‑1 (MQ=255)
                                                         gCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCa               >  1:345967/1‑70 (MQ=255)
                                                          cGACATTAATTGTGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCaaa             <  1:3957644/71‑1 (MQ=255)
                                                           gACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGCCaaaa            <  1:3587695/71‑1 (MQ=255)
                                                             cATTAATTGTGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCt                                      <  1:785822/43‑1 (MQ=255)
                                                                 aaTTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCaaaaa           >  1:3878335/1‑66 (MQ=255)
                                                                     gcgcTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCAttt   <  1:3186863/70‑1 (MQ=255)
                                                                     gcgcTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTc  >  1:3978408/1‑71 (MQ=255)
                                                                      |                                                                     
GCACAGCGTAGTGATATTCTTGATGCCGTGACACTGTCACTTAAAGTGGCGGCGCTGGCGACATTAATTGCGCTGGTGTTAGGGACGCTGGCAGCTGCCGCGCTGTGGCGACGAGACTTTTTCGGCAAAAACGCCATTTC  >  W3110S.gb/1516641‑1516780

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: