Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,526,495 T→C 12.9% M70V (ATG→GTG)  yjhF ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,526,4950TC12.9% 158.0 / 15.5 62M70V (ATG→GTG) yjhFpredicted transporter
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base C (5/3);  total (36/26)
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.81e-01

CCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGAC  >  W3110S.gb/4526425‑4526564
                                                                      |                                                                     
ccGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCAt                                                                       <  1:3601402/71‑1 (MQ=255)
  gAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATaa                                                                     >  1:135042/1‑71 (MQ=255)
     ataatTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAacac                                                                  >  1:37467/1‑71 (MQ=255)
      taatTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACAcc                                                                 <  1:3817296/71‑1 (MQ=255)
             gCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCa                                                          <  1:4595878/71‑1 (MQ=255)
               aTGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCaaa                                                        <  1:2970111/71‑1 (MQ=255)
               aTGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCaaa                                                        <  1:4200379/71‑1 (MQ=255)
                   tCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAatat                                                     >  1:3457055/1‑70 (MQ=255)
                     aGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCaa                                                  <  1:1181027/71‑1 (MQ=255)
                     aGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCa                                                   >  1:2197773/1‑70 (MQ=255)
                       cgcTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCa                                                          <  1:1268506/61‑1 (MQ=255)
                       cgcTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATAACAAAc                                                >  1:3733546/1‑71 (MQ=255)
                        gcTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAAcc                                               <  1:245188/71‑1 (MQ=255)
                         cTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAatat                                                     >  1:3997924/1‑64 (MQ=255)
                          tATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCa                                                          <  1:1856986/58‑1 (MQ=255)
                           aTTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCaa                                             >  1:4075835/1‑70 (MQ=255)
                           aTTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCaa                                             >  1:2189080/1‑70 (MQ=255)
                            ttCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATaa                                                                     >  1:1107847/1‑45 (MQ=255)
                            ttCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAata                                           >  1:2073974/1‑71 (MQ=255)
                            ttCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCaa                                             <  1:4430597/69‑1 (MQ=255)
                            ttCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCAata                                           >  1:2012660/1‑71 (MQ=255)
                              cGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAata                                           >  1:452165/1‑69 (MQ=255)
                              cGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATg                                         >  1:3872718/1‑71 (MQ=255)
                                       ccaccCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAAcc                                               <  1:1348787/56‑1 (MQ=255)
                                        caccCGCATCGGCCAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTa                               <  1:4191474/71‑1 (MQ=255)
                                         accCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCaa                                             >  1:264432/1‑56 (MQ=255)
                                         accCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATgg                                        >  1:1275817/1‑61 (MQ=255)
                                           ccGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAatat                                                     <  1:3275726/46‑1 (MQ=255)
                                            cGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCaa                                             >  1:620740/1‑53 (MQ=255)
                                             gCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCaaa                                                        >  1:222019/1‑41 (MQ=255)
                                             gCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGc                           >  1:2882465/1‑70 (MQ=255)
                                               aTCGGCCAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCaa                                             <  1:4295961/50‑1 (MQ=255)
                                                tCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAAc                                                <  1:2506942/46‑1 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAata                                           >  1:4548632/1‑49 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAat                                            >  1:2368344/1‑48 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTc                     >  1:3266640/1‑71 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTc                     >  1:512554/1‑71 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTc                     >  1:3794715/1‑71 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTc                     >  1:2244344/1‑71 (MQ=255)
                                                  ggCCAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTc                     >  1:3785260/1‑71 (MQ=255)
                                                    ccAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCg                   >  1:3044239/1‑71 (MQ=255)
                                                    ccAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCg                   <  1:2762853/71‑1 (MQ=255)
                                                     cAGCAAGCGACCGAGCACAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCg                   >  1:2439148/1‑70 (MQ=255)
                                                      aGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCgca                 <  1:3459172/71‑3 (MQ=255)
                                                        cAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGAtt               >  1:2216215/1‑71 (MQ=255)
                                                        cAAGCGACAGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGt                                <  1:1478290/54‑1 (MQ=255)
                                                         aaGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTc              >  1:1144977/1‑71 (MQ=255)
                                                         aaGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTc              <  1:3699725/71‑1 (MQ=255)
                                                           gCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGa            >  1:1474551/1‑71 (MQ=255)
                                                             gACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACt          >  1:3826562/1‑71 (MQ=255)
                                                             gACCGAGCACAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGa                 >  1:1873227/1‑64 (MQ=255)
                                                              aCCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTg         <  1:4564670/71‑1 (MQ=255)
                                                               ccGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCg                                 <  1:2732015/46‑1 (MQ=255)
                                                               ccGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCg                                 <  1:1292266/46‑1 (MQ=255)
                                                               ccGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCaa                        >  1:3400077/1‑55 (MQ=255)
                                                                 gAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAg       >  1:1201387/1‑70 (MQ=255)
                                                                  aGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAg                                   <  1:804091/41‑1 (MQ=255)
                                                                   gCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGt     >  1:4045110/1‑70 (MQ=255)
                                                                   gCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTg    <  1:3194584/71‑1 (MQ=255)
                                                                   gCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTg    <  1:2875534/71‑1 (MQ=255)
                                                                    cATAACACCGAAGCCAAATATCAAACCAATATGGc                                       >  1:1287086/1‑35 (MQ=255)
                                                                     aTAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCg                                 <  1:3004494/40‑1 (MQ=255)
                                                                     aTAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCg                                 <  1:671136/40‑1 (MQ=255)
                                                                      tAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCg             <  1:2317658/59‑1 (MQ=255)
                                                                      tAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTgac  >  1:4472988/1‑70 (MQ=255)
                                                                      tAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTgac  >  1:1364855/1‑70 (MQ=255)
                                                                      |                                                                     
CCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGAC  >  W3110S.gb/4526425‑4526564

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: