Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,039,152 T→C 9.8% Y274Y (TAT→TAC ftsX → predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,039,1520TC 9.8% 131.8 / 5.5 51Y274Y (TAT→TACftsXpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (25/21);  new base C (3/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCCGTGACTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGTGCTGCGAT  >  W3110S.gb/4039082‑4039222
                                                                      |                                                                      
gCCGTGACTCCTTTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTAt                                                                        <  1:3458469/71‑1 (MQ=255)
    tGACTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTAtggtg                                                                    <  1:864604/71‑1 (MQ=255)
    tGACTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTAtggtg                                                                    <  1:1265403/71‑1 (MQ=255)
    tGACTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTAtggtg                                                                    <  1:374838/71‑1 (MQ=255)
      aCTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGgc                                                                  >  1:1662345/1‑71 (MQ=255)
      aCTCCATTAACGTACAGAAACTGATTGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGgcg                                                                 >  1:578703/1‑71 (MQ=255)
       cTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGgcg                                                                 <  1:1141541/71‑1 (MQ=255)
       cTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGgcg                                                                 <  1:2514/71‑1 (MQ=255)
          cATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGgcgc                                                                <  1:2671813/69‑1 (MQ=255)
          cATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGgcg                                                                 <  1:3364462/68‑1 (MQ=255)
          cATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCAc                                                              >  1:575578/1‑71 (MQ=255)
          cATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCAc                                                              >  1:3248283/1‑71 (MQ=255)
             tAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgc                                                           >  1:175414/1‑71 (MQ=255)
              aaCGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgct                                                          >  1:3694884/1‑71 (MQ=255)
              aaCGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgct                                                          >  1:107488/1‑71 (MQ=255)
              aaCGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgct                                                          >  1:3367413/1‑71 (MQ=255)
              aaCGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgct                                                          >  1:2408800/1‑71 (MQ=255)
              aaCGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgct                                                          >  1:1066933/1‑71 (MQ=255)
              aaCGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCActgct                                                          >  1:115484/1‑71 (MQ=255)
                    cAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGAtt                                                    >  1:3630459/1‑71 (MQ=255)
                    cAGAAAATGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGAt                                                     <  1:453344/70‑1 (MQ=255)
                      gAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGAtttt                                                  >  1:3035174/1‑71 (MQ=255)
                       aaaCTGATTGGTGCGACGGATGGTTTCATCCTGCGCCCGTTCCTGTACGGTGGCGCACTGCTGGGATTTTc                                                 <  1:2573403/71‑1 (MQ=255)
                       aaaCTGATTGGTGCGACGGATGGTTTCATCCTGCGCCCGTTCCTGTACGGTGGCGCACTGCTGGGATTTTc                                                 <  1:2577846/71‑1 (MQ=255)
                       aaaCTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTc                                                 <  1:2623051/71‑1 (MQ=255)
                             aTTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGgcg                                            <  1:123055/70‑1 (MQ=255)
                                ggTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCAt                                         >  1:3626312/1‑70 (MQ=255)
                                 gTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCAtt                                        <  1:3664037/70‑1 (MQ=255)
                                 gTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCAtt                                        <  1:2459387/70‑1 (MQ=255)
                                 gTGCGACAGATGGATACATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCAttg                                       <  1:1923201/71‑1 (MQ=255)
                                       cAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCa                                                               >  1:2154641/1‑41 (MQ=255)
                                       cAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCa                                                               >  1:1771871/1‑41 (MQ=255)
                                       cAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCa                                                               >  1:3613855/1‑41 (MQ=255)
                                          aTGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCAtt                               <  1:3524740/70‑1 (MQ=255)
                                             gATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAAt                            <  1:1569887/70‑1 (MQ=255)
                                              aTTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCAtt                                        <  1:1907044/57‑1 (MQ=255)
                                               ttCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAAttt                          >  1:2344117/1‑70 (MQ=255)
                                                tCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTg                        >  1:1181389/1‑71 (MQ=255)
                                                 cATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAAtttt                         <  1:2129626/69‑1 (MQ=255)
                                                        cgcCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGaaa                 <  1:2609034/70‑1 (MQ=255)
                                                        cgcCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGaaa                 <  1:1320196/70‑1 (MQ=255)
                                                         gcCCGTTCCTGTATGTTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAAtt               <  1:2052421/71‑1 (MQ=255)
                                                         gcCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCa                                          >  1:2026032/1‑44 (MQ=255)
                                                           ccGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCt             <  1:339963/71‑1 (MQ=255)
                                                           ccGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCt             <  1:2323583/71‑1 (MQ=255)
                                                                ccTTTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGtgc        >  1:2649410/1‑71 (MQ=255)
                                                                ccTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGaaa                 >  1:3192790/1‑62 (MQ=255)
                                                                ccTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGtgc        >  1:1842012/1‑71 (MQ=255)
                                                                ccTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGtg         <  1:3447823/70‑1 (MQ=255)
                                                                ccTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGtg         <  1:2661785/70‑1 (MQ=255)
                                                                 cTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGtgct       >  1:1293824/1‑71 (MQ=255)
                                                                 cTGTACGGTGGCGCACTGCTGGGATTTTCTGgcgc                                           >  1:2607435/1‑35 (MQ=255)
                                                                 cTGTACGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAAt                            >  1:3389047/1‑50 (MQ=255)
                                                                 cTGTACGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAAt                            >  1:3483311/1‑50 (MQ=255)
                                                                   gTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGtgc        >  1:639639/1‑68 (MQ=255)
                                                                      tggtggCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTc              <  1:223668/59‑1 (MQ=255)
                                                                      tggtggCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTc              <  1:2676326/59‑1 (MQ=255)
                                                                      tggtggCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTc              <  1:3359678/59‑1 (MQ=255)
                                                                      tggtggCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTg            >  1:2661859/1‑61 (MQ=255)
                                                                      tggtggCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGTGCTGCGa   <  1:1524397/70‑1 (MQ=255)
                                                                      tggtggCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGTGCTGCGAt  >  1:1214305/1‑71 (MQ=255)
                                                                      |                                                                      
GCCGTGACTCCATTAACGTACAGAAACTGATTGGTGCGACAGATGGATTCATCCTGCGCCCGTTCCTGTATGGTGGCGCACTGCTGGGATTTTCTGGCGCATTGTTGTCATTAATTTTGTCAGAAATTCTGGTGCTGCGAT  >  W3110S.gb/4039082‑4039222

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