Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,164,775 C→T 63.2% A534T (GCA→ACA)  ygiS ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,164,7750CT63.2% 45.9 / 60.0 57A534T (GCA→ACA) ygiSpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base T (14/22);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

TCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA  >  W3110S.gb/3164710‑3164841
                                                                 |                                                                  
tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCtt                                                                 <  1:1015260/69‑1 (MQ=255)
tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCtt                                                                 <  1:1890891/69‑1 (MQ=255)
tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCtt                                                                 <  1:1682553/69‑1 (MQ=255)
 cacaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTg                                                                <  1:1025088/69‑1 (MQ=255)
  acaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGa                                                               >  1:1016204/1‑69 (MQ=255)
  acaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGa                                                               <  1:1525563/69‑1 (MQ=255)
  acaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGa                                                               >  1:75943/1‑69 (MQ=255)
  aaaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGa                                                               >  1:822905/3‑69 (MQ=255)
       agaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATAc                                                          <  1:2447668/69‑1 (MQ=255)
        gaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATAc                                                          <  1:1843966/68‑1 (MQ=255)
                     aGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTg                                              <  1:1177539/67‑1 (MQ=255)
                      gcgcCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAg                                           <  1:215735/69‑1 (MQ=255)
                            ccgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATaa                                     <  1:299057/69‑1 (MQ=255)
                            ccgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATa                                      <  1:3095095/68‑1 (MQ=255)
                             cgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATaa                                     <  1:2476280/68‑1 (MQ=255)
                                 cGTGCTCTGACGGCGATAAAGGCATCAATGTGTCtt                                                                 >  1:1807905/1‑36 (MQ=37)
                                 cGTGCTCCGACGGCGATAAAGGCATCAATGTGTCtt                                                                 >  1:2582802/1‑36 (MQ=255)
                                   tGCTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTggg                              <  1:2805813/69‑1 (MQ=255)
                                     cTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                             >  1:942670/1‑68 (MQ=255)
                                     cTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                             >  1:2858721/1‑68 (MQ=255)
                                     cTCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                             >  1:1280526/1‑68 (MQ=255)
                                      tCCGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                             >  1:2908682/1‑67 (MQ=255)
                                       ccGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAt                           <  1:3299304/68‑1 (MQ=255)
                                       ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                          >  1:1830724/1‑69 (MQ=255)
                                       ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                          >  1:1750155/1‑69 (MQ=255)
                                       ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                          >  1:1745798/1‑69 (MQ=255)
                                        cGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACa                                        <  1:1568186/54‑1 (MQ=255)
                                        cGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACa                                        <  1:3511916/54‑1 (MQ=255)
                                        cGACGGCGATAAAGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACa                                        <  1:1733602/54‑1 (MQ=255)
                                                 tAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTg                                              <  1:3444932/39‑1 (MQ=255)
                                                  aaaGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACa                                        <  1:1770215/44‑1 (MQ=255)
                                                  aaaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaa               <  1:2113437/69‑1 (MQ=255)
                                                  aaaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaa               <  1:199050/69‑1 (MQ=255)
                                                   aaGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa              >  1:3275267/1‑69 (MQ=255)
                                                   aaGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa              >  1:1118764/1‑69 (MQ=255)
                                                   aaGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa              <  1:499123/69‑1 (MQ=255)
                                                   aaGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa              <  1:1995966/69‑1 (MQ=255)
                                                   aaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa              >  1:980253/1‑69 (MQ=255)
                                                    aGGCATCAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaa              <  1:236996/68‑1 (MQ=255)
                                                          cAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGcc        >  1:3396582/1‑68 (MQ=255)
                                                          cAATGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGcc        >  1:3400777/1‑68 (MQ=255)
                                                           aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa      >  1:2219935/1‑69 (MQ=255)
                                                           aaTGTGTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCa       >  1:3719812/1‑68 (MQ=255)
                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa      >  1:3300251/1‑69 (MQ=255)
                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa      >  1:11999/1‑69 (MQ=255)
                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa      >  1:3461785/1‑69 (MQ=255)
                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa      >  1:2035514/1‑69 (MQ=255)
                                                           aaTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCaa      >  1:433219/1‑69 (MQ=255)
                                                             tgtgTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATg                       <  1:1563431/50‑1 (MQ=255)
                                                             tgtgTCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATg                       <  1:1733564/50‑1 (MQ=255)
                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                  <  1:1130798/38‑1 (MQ=255)
                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                  <  1:2753114/38‑1 (MQ=255)
                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATcc                                  <  1:3304839/38‑1 (MQ=255)
                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGa                            >  1:2521399/1‑44 (MQ=255)
                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAt                           >  1:3440666/1‑45 (MQ=255)
                                                              gtgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAt                           >  1:3436375/1‑45 (MQ=255)
                                                               tgCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa  <  1:1432060/69‑1 (MQ=255)
                                                                 |                                                                  
TCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA  >  W3110S.gb/3164710‑3164841

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: