Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,989,353 T→A 6.3% S98T (TCT→ACT)  yecI → predicted ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,989,3530TA 6.3% 171.4 / 4.0 63S98T (TCT→ACT) yecIpredicted ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base A (2/2);  total (34/29)
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.48e-01

GTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCT  >  W3110S.gb/1989284‑1989420
                                                                     |                                                                   
gTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTc                                                                    >  1:2152622/1‑71 (MQ=255)
 ttCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTc                                                                    <  1:1676218/70‑1 (MQ=255)
     cGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTa                                                               >  1:4070895/1‑71 (MQ=255)
       gTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACg                                                             >  1:3106574/1‑71 (MQ=255)
       gTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACg                                                             >  1:3101792/1‑71 (MQ=255)
           aaaaCTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTgg                                                         >  1:1207154/1‑71 (MQ=255)
              aCTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTgg                                                         >  1:1005486/1‑68 (MQ=255)
              aCTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcac                                                      <  1:4037170/71‑1 (MQ=255)
               cTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGca                                                       >  1:1724183/1‑69 (MQ=255)
                      ctctGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGcc                                              >  1:1411295/1‑71 (MQ=255)
                         tGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGAt                                           >  1:1484627/1‑71 (MQ=255)
                          ggAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcac                                                      <  1:2299111/59‑1 (MQ=255)
                           gaagaaCTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATg                                          <  1:2397629/70‑1 (MQ=255)
                            aagaaCTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATg                                          >  1:3668621/1‑69 (MQ=255)
                                aCTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGa                                    <  1:975676/71‑1 (MQ=255)
                                 cTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGaa                                   <  1:3583339/71‑1 (MQ=255)
                                     tCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGaa                               >  1:2712975/1‑71 (MQ=255)
                                      ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACg                                                             >  1:3045227/1‑40 (MQ=255)
                                      ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcaca                                                     >  1:2591661/1‑48 (MQ=255)
                                      ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGaa                               <  1:481235/70‑1 (MQ=255)
                                      ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAAc                              >  1:1975969/1‑71 (MQ=255)
                                        aaaaaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGa                                         <  1:1924712/58‑1 (MQ=255)
                                        aaaaaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                             >  1:2311721/1‑70 (MQ=255)
                                         aaaaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcac                                                      <  1:602557/44‑1 (MQ=255)
                                           aaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAat                         <  1:287349/71‑1 (MQ=255)
                                            aaCGATGGAAGAATACGAGCAACGTTCTAGTAGGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAat                         <  1:564947/70‑1 (MQ=255)
                                            aaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatg                        <  1:2378276/71‑1 (MQ=255)
                                              cGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgat                      <  1:788572/71‑1 (MQ=255)
                                               gATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                     <  1:3424932/71‑1 (MQ=255)
                                               gATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                     <  1:1526557/71‑1 (MQ=255)
                                                 tGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                             >  1:67219/1‑61 (MQ=255)
                                                 tGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                             >  1:2580830/1‑61 (MQ=255)
                                                 tGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                             >  1:2575722/1‑61 (MQ=255)
                                                    aagaaTACGATCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGaa                               >  1:2597295/1‑56 (MQ=255)
                                                     agaaTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                     >  1:1640498/1‑65 (MQ=255)
                                                       aaTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAAcc             <  1:2824126/71‑1 (MQ=255)
                                                        aTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAAcc             <  1:841721/70‑1 (MQ=255)
                                                        aTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCg            <  1:1975789/71‑1 (MQ=255)
                                                         tACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAg                                       >  1:2863778/1‑43 (MQ=255)
                                                         tACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt           >  1:638149/1‑71 (MQ=255)
                                                          aCGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCg                                             >  1:2374921/1‑36 (MQ=255)
                                                          aCGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt           >  1:1190295/1‑70 (MQ=255)
                                                           cGAGCAACGTTCTAGTACGTTGTCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTc          >  1:582570/1‑70 (MQ=255)
                                                           cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAg                                 >  1:3425003/1‑47 (MQ=255)
                                                           cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAg                                 >  1:1128228/1‑47 (MQ=255)
                                                           cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAAc              >  1:3879954/1‑66 (MQ=255)
                                                           cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTc          >  1:2209999/1‑70 (MQ=255)
                                                           cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCa         <  1:1887058/71‑1 (MQ=255)
                                                              gCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAAtt      <  1:341793/71‑1 (MQ=255)
                                                              ccAACGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt           <  1:1604204/65‑1 (MQ=39)
                                                              ccAACGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt           <  1:3826346/65‑1 (MQ=39)
                                                               cAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAAtgatga                    >  1:100485/1‑56 (MQ=255)
                                                                aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGa                                    <  1:1008349/39‑1 (MQ=255)
                                                                aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTg                            <  1:3066032/47‑1 (MQ=255)
                                                                aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAAttt     <  1:4075592/70‑1 (MQ=255)
                                                                aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTc    <  1:806378/71‑1 (MQ=255)
                                                                aaCGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCa         >  1:3293269/1‑66 (MQ=255)
                                                                aaCGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCa         >  1:1578734/1‑66 (MQ=255)
                                                                 aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                             >  1:3815623/1‑45 (MQ=255)
                                                                 aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                     <  1:2504854/53‑1 (MQ=255)
                                                                 aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTc    <  1:451451/70‑1 (MQ=255)
                                                                  cGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCt  >  1:3366901/1‑71 (MQ=255)
                                                                   gTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTc          <  1:1996702/62‑1 (MQ=255)
                                                                    ttCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCaa        >  1:317425/1‑63 (MQ=255)
                                                                     |                                                                   
GTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCT  >  W3110S.gb/1989284‑1989420

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: