Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,688,036 T→C 71.8% I43V (ATT→GTT)  yfhG ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,688,0360TC71.8% 74.0 / 30.2 46I43V (ATT→GTT) yfhGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base C (16/17);  total (23/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 = 1.00e+00

CAGCGCCCAGATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGA  >  W3110S.gb/2687977‑2688103
                                                           |                                                                   
cAGCGCCCAGATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTctg                                                          <  1:2135852/71‑1 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAActtttt                                                                >  1:2969164/1‑56 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAActtttt                                                                >  1:2597143/1‑56 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAActtttt                                                                >  1:2453538/1‑56 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAActtttt                                                                >  1:1943055/1‑56 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAActtttt                                                                >  1:1822019/1‑56 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTAt                                                 <  1:3045875/71‑1 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTAt                                                 <  1:2357200/71‑1 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTAt                                                 <  1:2244010/71‑1 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTAt                                                 <  1:2006952/71‑1 (MQ=255)
         gATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTAt                                                 <  1:616276/71‑1 (MQ=255)
          aTGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTAt                                                 >  1:410653/1‑70 (MQ=255)
          aTGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATc                                                <  1:353243/71‑1 (MQ=255)
          aTGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATc                                                <  1:2098183/71‑1 (MQ=255)
          aTGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATc                                                <  1:930363/71‑1 (MQ=255)
            gTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCaa                                              >  1:2100574/1‑71 (MQ=255)
            gTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCaa                                              >  1:569077/1‑71 (MQ=255)
            gTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCaa                                              >  1:3022696/1‑71 (MQ=255)
               gCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAAtg                                            <  1:3227831/70‑1 (MQ=255)
                   cATTCAGTTGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAAtggctg                                        >  1:2878707/1‑70 (MQ=255)
                   cATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTctg                                                          >  1:1104467/1‑52 (MQ=255)
                   cATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAAtggctgg                                       >  1:93326/1‑71 (MQ=255)
                   cATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAAtggctgg                                       >  1:2876018/1‑71 (MQ=255)
                   cATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAAtggctgg                                       >  1:2620401/1‑71 (MQ=255)
                   cATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAAtggctg                                        >  1:1053986/1‑70 (MQ=255)
                     ttCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATggctggct                                     >  1:371018/1‑71 (MQ=255)
                     ttCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATggctggct                                     >  1:3381444/1‑71 (MQ=255)
                          tAGAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATtg                                >  1:3237909/1‑71 (MQ=255)
                            gAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATtgt                               <  1:501489/70‑1 (MQ=255)
                            gAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTttctt                                                             >  1:3051270/1‑40 (MQ=255)
                            gAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTttctt                                                             >  1:3359309/1‑40 (MQ=255)
                            gAAAGATAGTCAGCCAGTTGATAAACCGGAACTTTttctt                                                             >  1:2489316/1‑40 (MQ=255)
                               aGATAGTCAGCCAGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGAtt                           >  1:969008/1‑71 (MQ=255)
                                           aGTTGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGa                             >  1:2630752/1‑57 (MQ=255)
                                              tGATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAg            >  1:1131576/1‑71 (MQ=255)
                                               gATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAg            >  1:2974332/1‑70 (MQ=255)
                                               gATAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAg            >  1:548565/1‑70 (MQ=255)
                                                aTAAACCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGc           <  1:1810123/70‑1 (MQ=255)
                                                aTAAACCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGc           <  1:1154260/70‑1 (MQ=255)
                                                aTAAACCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGc           <  1:1726717/70‑1 (MQ=255)
                                                aTAAACCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGc           <  1:2850532/70‑1 (MQ=255)
                                                  aaaCCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCaa         >  1:2536137/1‑70 (MQ=255)
                                                   aaCCGGAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATg       >  1:1989967/1‑71 (MQ=255)
                                                    aCCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGc      <  1:1884875/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:1077421/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:860998/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:1359535/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:1722547/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:434903/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:3542677/71‑1 (MQ=255)
                                                        gAACTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGa  <  1:2262270/71‑1 (MQ=255)
                                                           |                                                                   
CAGCGCCCAGATGTCGCTACATTCAGTAGAAAGATAGTCAGCCAGTTGATAAACCGGAATTTTTTCTTCTGCTGGCGTATCAATGGCTGGCTTATTGTGATTCTGCACGCAACCCAGCAATGCCAGA  >  W3110S.gb/2687977‑2688103

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: