Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,559,087 T→C 13.9% V222V (GTA→GTG ddpF ← D‑Ala‑D‑Ala transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,559,0870TC13.9% 184.4 / 20.7 72V222V (GTA→GTGddpFD‑Ala‑D‑Ala transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (34/28);  new base C (5/5);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGA  >  W3110S.gb/1559018‑1559155
                                                                     |                                                                    
ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTa                                                                     >  1:2456439/1‑71 (MQ=255)
 gaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTAc                                                                    <  1:722692/71‑1 (MQ=255)
 gaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTAc                                                                    <  1:3476058/71‑1 (MQ=255)
  aaTCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACa                                                                   <  1:3550900/71‑1 (MQ=255)
          aataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGccc                                                           >  1:2704397/1‑71 (MQ=255)
            taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACa                                                                   <  1:3385541/61‑1 (MQ=255)
            taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATc                                                                <  1:74311/64‑1 (MQ=255)
              atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTAc                                                                    >  1:2828265/1‑58 (MQ=255)
                      tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                     <  1:1298509/65‑1 (MQ=255)
                        tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                             >  1:1666441/1‑71 (MQ=255)
                        tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                             >  1:628287/1‑71 (MQ=255)
                          gATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGc                                            >  1:4003454/1‑70 (MQ=255)
                          gATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGc                                            >  1:3838371/1‑70 (MQ=255)
                              tgCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAccc                                       >  1:3273402/1‑71 (MQ=255)
                                cAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGa                                     >  1:1022840/1‑71 (MQ=255)
                                cAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGa                                     <  1:1270310/71‑1 (MQ=255)
                                cAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGa                                     >  1:3243831/1‑71 (MQ=255)
                                  ggTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATAc                                                    >  1:1812627/1‑54 (MQ=255)
                                   gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCga                                                         >  1:3598842/1‑48 (MQ=255)
                                   gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                      <  1:1477793/51‑1 (MQ=255)
                                   gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                  >  1:4398754/1‑71 (MQ=255)
                                   gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                  >  1:924600/1‑71 (MQ=255)
                                   gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCCACAATCTGCCCGAGATACATCAcc                                              <  1:1283928/59‑1 (MQ=255)
                                   gTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCCACAATCTGCCCGAGATACATCAcc                                              <  1:3506755/59‑1 (MQ=255)
                                    tGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                   >  1:3801243/1‑69 (MQ=255)
                                      ccGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCg                                      <  1:1788286/64‑1 (MQ=255)
                                        gTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                       <  1:1304903/45‑1 (MQ=255)
                                         tCAGCACCTGCTGCGCGTCGCCCAGTTCCACAATCTGcc                                                            >  1:2677535/1‑39 (MQ=255)
                                         tCAGCACCTGCTGCGCGTCGCCCAGTTCCACAATCTGcc                                                            >  1:834002/1‑39 (MQ=255)
                                            tcACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCt                                <  1:3405183/63‑1 (MQ=255)
                                             cACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATc                                                 <  1:3411901/46‑1 (MQ=255)
                                             cACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                        <  1:2194596/71‑1 (MQ=255)
                                              aCCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                              >  1:4102020/1‑64 (MQ=255)
                                               ccTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCata                            <  1:1315642/65‑1 (MQ=255)
                                               ccTGCTGCGCGTCGCCCAGTTCCACAATCTGCCCGAGATACATCACCGCCACCCGATCGCt                                <  1:1958736/61‑1 (MQ=255)
                                                 tgctgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTAt                    >  1:859183/1‑71 (MQ=255)
                                                 tgctgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTAt                    >  1:4245036/1‑71 (MQ=255)
                                                 tgctgcGCGTCGCCCAGTTCCACAATCTGCCCGAGATa                                                     <  1:2325808/38‑1 (MQ=255)
                                                  gctgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTAt                    >  1:4057831/1‑70 (MQ=255)
                                                   ctgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATc                   <  1:3959011/70‑1 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAccc                                       <  1:3512348/49‑1 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAccc                                       <  1:352124/49‑1 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAccc                                       <  1:4609218/49‑1 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAccc                                       <  1:12167/49‑1 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAccc                                       <  1:155336/49‑1 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCAc                 >  1:1901369/1‑71 (MQ=255)
                                                    tgcGCGTCGCCCAGTTCCACAATCTGCCCGAGATACATc                                                 <  1:3310125/39‑1 (MQ=255)
                                                      cgcgTCGCCCAGTTCCACAATCTGCCCGAGATACATCAc                                               >  1:270133/1‑39 (MQ=255)
                                                          tCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                         >  1:2124137/1‑41 (MQ=255)
                                                           cGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCa                                          >  1:1863695/1‑39 (MQ=255)
                                                            gCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACa          <  1:2015622/70‑1 (MQ=255)
                                                            gCCCAGTTCCACAATCTGCCCGAGATAGATCACCGCCACCCGATCGCTCatat                           >  1:686375/1‑53 (MQ=25)
                                                               cAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACg                       <  1:1959327/54‑1 (MQ=255)
                                                               cAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACg                       <  1:55748/54‑1 (MQ=255)
                                                               cAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACg                       <  1:2017409/54‑1 (MQ=255)
                                                                 gTTCCACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                              >  1:4304456/1‑45 (MQ=255)
                                                                  ttCTACAATCTGCCCGAGATACATCACCGCCACCCg                                      <  1:2002965/36‑1 (MQ=255)
                                                                  ttCTACAATCTGCCCGAGATACATCACCGCCACCCg                                      <  1:125124/36‑1 (MQ=255)
                                                                  ttCTACAATCTGCCCGAGATACATCACCGCCACCCg                                      <  1:1057732/36‑1 (MQ=255)
                                                                  ttCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTAt                    >  1:1797559/1‑54 (MQ=255)
                                                                  ttCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTAt                    >  1:1340044/1‑54 (MQ=255)
                                                                  ttCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgtg   <  1:3444061/71‑1 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCata                            >  1:300314/1‑44 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCata                            >  1:2698089/1‑44 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                        >  1:2633415/1‑48 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACg                       >  1:976086/1‑49 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgaga            >  1:3020957/1‑60 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgaga            >  1:2546000/1‑60 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgaga            >  1:156080/1‑60 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgaga            >  1:1211046/1‑60 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa  >  1:3893085/1‑70 (MQ=255)
                                                                    cTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa  >  1:3360136/1‑70 (MQ=255)
                                                                     |                                                                    
GGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGA  >  W3110S.gb/1559018‑1559155

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: