Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,276,691 Δ1 bp 100% coding (994/1815 nt) yejA → predicted oligopeptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,276,6910T.83.3% 121.5 / 25.6 42coding (994/1815 nt)yejApredicted oligopeptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (7/0);  new base . (20/15);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 3.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGTTA  >  W3110S.gb/2276623‑2276752
                                                                    |                                                             
ggctggcGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCct                                                               >  1:3031700/1‑69 (MQ=255)
    ggcGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCcttgac                                                           >  1:2407974/1‑66 (MQ=255)
        tttAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      <  1:3357699/69‑1 (MQ=255)
        tttAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgactt                                                        >  1:679978/1‑68 (MQ=255)
         ttAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgactttg                                                      <  1:2555683/69‑1 (MQ=255)
         ttAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgactttg                                                      <  1:698727/69‑1 (MQ=255)
          tAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgactttga                                                     >  1:3110170/1‑69 (MQ=255)
                 cAACGTCCGGTATTCAGCGCTCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATg                                              >  1:3026049/1‑68 (MQ=255)
                 cAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      <  1:345840/60‑1 (MQ=255)
                 cAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      <  1:2171321/60‑1 (MQ=255)
                 cAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      <  1:507824/60‑1 (MQ=255)
                 cAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGa                                             <  1:1169708/69‑1 (MQ=255)
                     gTCGGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATg                                                  <  1:2049372/60‑1 (MQ=255)
                        cGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAgg                                       >  1:681056/1‑69 (MQ=255)
                             ttCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCgttg                                  >  1:3275206/1‑69 (MQ=255)
                             ttCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCgttg                                  >  1:3686479/1‑69 (MQ=255)
                              tCAGCGATCGCCGGGTCCGTGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCgttgt                                 <  1:310238/68‑1 (MQ=255)
                              tCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCgttgt                                 <  1:1865787/68‑1 (MQ=255)
                              tCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCgttgt                                 <  1:1484622/68‑1 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgact                                                         >  1:3630403/1‑44 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCctttgact                                                         >  1:2797118/1‑44 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACa                                          >  1:2519247/1‑59 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCgttgt                                 >  1:2241904/1‑68 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGTATGAACAAGGCGTTGttt                               <  1:488967/69‑1 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCgttgtt                                >  1:3536052/1‑68 (MQ=255)
                               cAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGttt                               >  1:3018183/1‑69 (MQ=255)
                                  cGATCGCCGGGTCCGTGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTAc                            <  1:360635/69‑1 (MQ=255)
                                  cGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGTATGAACAAGGCGTTGTTTTa                             <  1:2017163/68‑1 (MQ=255)
                                  cGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGGTGAACAAGGCGTTGTTTTAc                            <  1:718055/69‑1 (MQ=255)
                                  cGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTa                             <  1:244780/68‑1 (MQ=255)
                                  cGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTAc                            <  1:2233826/69‑1 (MQ=255)
                                  cGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTAc                            <  1:1051932/69‑1 (MQ=255)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      >  1:2092730/1‑39 (MQ=38)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      <  1:148395/39‑1 (MQ=38)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGccttgactttg                                                      <  1:876148/39‑1 (MQ=38)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAAt                         >  1:1880477/1‑68 (MQ=255)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAAt                         >  1:1567966/1‑68 (MQ=255)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATg                        >  1:1068319/1‑69 (MQ=255)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATg                        >  1:742372/1‑69 (MQ=255)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATg                        >  1:838375/1‑69 (MQ=255)
                                      cGCCGGGTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATg                        >  1:3337334/1‑69 (MQ=255)
                                          gggTCCGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGcc                      >  1:3606685/1‑67 (MQ=255)
                                              ccGGGAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGt                >  1:712166/1‑69 (MQ=255)
                                                gggAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCg              >  1:1154456/1‑69 (MQ=255)
                                                gggAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCg              >  1:2526450/1‑69 (MQ=255)
                                                gggAAGCTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCg              >  1:830565/1‑69 (MQ=255)
                                                      cTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGaa        <  1:3007162/69‑1 (MQ=255)
                                                      cTATCACTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGaa        <  1:1317528/69‑1 (MQ=255)
                                                           aCTCTCGCCTTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTg                    >  1:2167467/1‑53 (MQ=255)
                                                           aCTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGttt                               >  1:3659011/1‑41 (MQ=38)
                                                           aCTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGttt                               >  1:2128058/1‑41 (MQ=38)
                                                           aCTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGttt                               >  1:1632618/1‑41 (MQ=38)
                                                           aCTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGtt   >  1:977366/1‑69 (MQ=255)
                                                           aCTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGtt   >  1:1898150/1‑69 (MQ=255)
                                                           aCTCTCGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGtt   >  1:3207543/1‑69 (MQ=255)
                                                            ctctCGCCTTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGtt   >  1:1387241/1‑69 (MQ=255)
                                                             tctcGCCTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGTTa  <  1:1140433/68‑1 (MQ=255)
                                                                    tttgactttgaATGGATGAACAAGGCGTTGTTTTAc                            <  1:390546/36‑1 (MQ=255)
                                                                    |                                                             
GGCTGGCGTTTAATATCCAACGTCCGGTATTCAGCGATCGCCGGGTCCGGGAAGCTATCACTCTCGCCTTTGACTTTGAATGGATGAACAAGGCGTTGTTTTACAATGCCTGGAGTCGCACGAACAGTTA  >  W3110S.gb/2276623‑2276752

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: