Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,724,748 T→C 26.8% V305A (GTC→GCC)  ydhK → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,724,7480TC26.8% 62.8 / 26.9 41V305A (GTC→GCC) ydhKconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (17/13);  new base C (5/6);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGT  >  W3110S.gb/1724686‑1724818
                                                              |                                                                      
cAGTTCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGCCGCCGACt                                                                <  1:1628882/71‑1 (MQ=255)
cAGTTCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGCCGCCGACt                                                                <  1:4985565/71‑1 (MQ=255)
 aGTTCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTa                                                               <  1:824971/71‑1 (MQ=255)
    tCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTAt                                                              >  1:4427037/1‑69 (MQ=255)
    tCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGCCGCCGACTATCg                                                            >  1:748814/1‑71 (MQ=255)
       cAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAAcgccgc                                                                     >  1:4464669/1‑59 (MQ=255)
       cAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAAcgccgc                                                                     >  1:256931/1‑59 (MQ=255)
       cAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAAcgccgc                                                                     >  1:250801/1‑59 (MQ=255)
        aaaCAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCAc                                                        <  1:3546296/71‑1 (MQ=255)
          aCAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGt                                                      >  1:1323931/1‑71 (MQ=255)
                gTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCtt                                                >  1:5098088/1‑71 (MQ=255)
                  ttACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCt                                                 <  1:1885869/68‑1 (MQ=255)
                  ttACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCt                                                 <  1:3736853/68‑1 (MQ=255)
                   tACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGCCGCCGACTATCGGCACGTCGCCTTCTg                                             >  1:2345969/1‑71 (MQ=255)
                        cGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGc                                        <  1:3226147/71‑1 (MQ=255)
                        cGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGc                                        <  1:2155287/71‑1 (MQ=255)
                            gCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGAc                                     >  1:5001419/1‑70 (MQ=255)
                                 tattatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGt                                >  1:1099309/1‑70 (MQ=255)
                                  attatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTa                              >  1:397323/1‑71 (MQ=255)
                                  attatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTa                              >  1:4492415/1‑71 (MQ=255)
                                   ttatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTa                              <  1:3251137/70‑1 (MQ=255)
                                   ttatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAt                             <  1:3201870/71‑1 (MQ=255)
                                   ttatCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAt                             <  1:3665747/71‑1 (MQ=255)
                                   ttatCGCCCCGCTACGCCCGACCAACGCCGCCGACTATCgg                                                           <  1:2907746/41‑1 (MQ=255)
                                   ttatCGCCCCGCTACGCCCGACCAACGCCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAt                             >  1:2356172/1‑71 (MQ=255)
                                     atCGCCCCGCTGCGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAttt                           >  1:3738196/1‑71 (MQ=255)
                                     atCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAttt                           >  1:886915/1‑71 (MQ=255)
                                     atCGCCCCGCTACGCCCGACCAACGCCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAttt                           >  1:3126299/1‑71 (MQ=255)
                                      tCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAtttt                          <  1:288772/71‑1 (MQ=255)
                                        gCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGCCTACGTTAttttt                         <  1:5013714/70‑1 (MQ=255)
                                        gCCCCGCTACGCCCGACCAACGCCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAttttt                         <  1:634363/70‑1 (MQ=255)
                                        gCCCCGCTACGCCCGACCAACGCCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTAttttt                         <  1:2405494/70‑1 (MQ=255)
                                           ccGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTgcc                      >  1:3118554/1‑70 (MQ=255)
                                            cGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGa                                      >  1:1426729/1‑53 (MQ=255)
                                                aCGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCtt                 >  1:10999/1‑70 (MQ=255)
                                                aCGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCtt                 >  1:2764212/1‑70 (MQ=255)
                                                 cGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCttt                <  1:4955413/70‑1 (MQ=255)
                                                   cccGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATc             <  1:3193417/71‑1 (MQ=255)
                                                      gACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTgc                       >  1:52292/1‑58 (MQ=255)
                                                       aCCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGc          >  1:4325007/1‑70 (MQ=255)
                                                        ccAACGCCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTg                        >  1:4521006/1‑55 (MQ=255)
                                                          aaCGCCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCAAag      <  1:2135074/71‑1 (MQ=255)
                                                           aCGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTgcc                      >  1:4241372/1‑54 (MQ=255)
                                                           aCGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATc             >  1:2686055/1‑63 (MQ=255)
                                                              tCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCa         <  1:2617967/64‑1 (MQ=255)
                                                              tCGCCGACTATCGGCACGTCGCCTTCTGGCAGCAACTACGTTATTTTTGCCGCCTTTATCTGCAAagtagt  >  1:2190085/1‑71 (MQ=255)
                                                              |                                                                      
CAGTTCGCAAACAGATGTTTACACCGTCGCACGTATTATCGCCCCGCTACGCCCGACCAACGTCGCCGACTATCGGCACGTCGCCTTCTGGCAGCGACTACGTTATTTTTGCCGCCTTTATCTGCAAAGTAGT  >  W3110S.gb/1724686‑1724818

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: