Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,221,135 T→C 15.4% T253A (ACC→GCC)  yehY ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,221,1350TC15.4% 93.0 / 10.9 39T253A (ACC→GCC) yehYpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (18/15);  new base C (4/2);  total (22/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

ACGCCGCGCACCAGCGGCAGCAAGGCATAGAGCACCAGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCAACA  >  W3110S.gb/2221071‑2221203
                                                                |                                                                    
aCGCCGCGCACCAGCGGCAGCAAGGCATAGAGCACCAGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGc                                                                >  1:3020660/1‑71 (MQ=255)
                         cATAGAGCACCAGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGa                                       <  1:3444429/71‑1 (MQ=255)
                           tAGAGCACCAGCGCAATCAGTGCGGGTGTCATTCCGGCTCCTGCTATg                                                            >  1:3507465/1‑48 (MQ=255)
                           tAGAGCACCAGCGCAATCAGTGCGGGTGTCATTCCGGCTCCTGCTATGCCGAGCGTCCCCAGCCACGGa                                       >  1:3030026/1‑69 (MQ=255)
                             gagCACCAGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGaaa                                     >  1:2135355/1‑69 (MQ=255)
                             gagCACCAGCGCAATCAGTGCGGGTGTCATTCCGGTTACTGCTATGCCGAGCGTCCCCAGc                                             >  1:3768465/1‑61 (MQ=255)
                                    aGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAAc                            >  1:926032/1‑71 (MQ=255)
                                    aGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAAc                            <  1:2135785/71‑1 (MQ=255)
                                    aGCGCAATCAGTGCGGGTGTCATTACGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTaa                             >  1:948171/1‑70 (MQ=255)
                                        cAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAgc                        >  1:2227521/1‑71 (MQ=255)
                                          aTCAGTGCGGGTGTCATTCCGGCTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAgcgc                      >  1:2778766/1‑71 (MQ=255)
                                            cAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCgg                    >  1:1552876/1‑71 (MQ=255)
                                              gTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCa                  <  1:199851/71‑1 (MQ=255)
                                              gTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCa                  <  1:627266/71‑1 (MQ=255)
                                                 cGGGTGTCATTCCGGCTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAgc               <  1:468328/71‑1 (MQ=255)
                                                     tgtCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGa                                       <  1:830280/43‑1 (MQ=255)
                                                     tgtCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGa                                       <  1:1204837/43‑1 (MQ=255)
                                                     tgtCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGa                                       <  1:1404247/43‑1 (MQ=255)
                                                     tgtCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAgcggcg           >  1:1737332/1‑71 (MQ=255)
                                                      gtCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAgcggcg           >  1:1472557/1‑70 (MQ=255)
                                                       tCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGgcgc          >  1:2758707/1‑70 (MQ=255)
                                                       tCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGgcgc          >  1:968351/1‑70 (MQ=255)
                                                         aTTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAAt       >  1:2688989/1‑71 (MQ=255)
                                                         aTTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAAt       >  1:864448/1‑71 (MQ=255)
                                                         aTTCCGGCTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAAt       >  1:932733/1‑71 (MQ=255)
                                                          ttCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATc      <  1:102757/71‑1 (MQ=255)
                                                          ttCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATc      <  1:2432417/71‑1 (MQ=255)
                                                          ttCCGGCTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATc      <  1:3075690/71‑1 (MQ=255)
                                                            ccGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    >  1:3527770/1‑71 (MQ=255)
                                                            ccGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    >  1:3892076/1‑71 (MQ=255)
                                                            ccGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    >  1:3034363/1‑71 (MQ=255)
                                                            ccGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    >  1:961721/1‑71 (MQ=255)
                                                            ccGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCa     <  1:719058/70‑1 (MQ=255)
                                                            ccGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCACGGCAAGCGGCGCAATCaa    >  1:3417761/1‑71 (MQ=255)
                                                             cGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    <  1:1891156/70‑1 (MQ=255)
                                                             cGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    <  1:3583840/70‑1 (MQ=255)
                                                             cGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    <  1:1304762/70‑1 (MQ=255)
                                                             cGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCaa    <  1:1300065/70‑1 (MQ=255)
                                                              ggTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCAACa  <  1:1456027/71‑1 (MQ=255)
                                                                |                                                                    
ACGCCGCGCACCAGCGGCAGCAAGGCATAGAGCACCAGCGCAATCAGTGCGGGTGTCATTCCGGTTCCTGCTATGCCGAGCGTCCCCAGCCACGGAAAGGCCGTAACCAGCGCGGCAAGCGGCGCAATCAACA  >  W3110S.gb/2221071‑2221203

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