Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,473,932 G→T 11.9% pseudogene (6827/8622 nt)
pseudogene (1529/3324 nt)
ydbA →
ydbA →
ECK1398:JW5802+JW1402:b4492; predicted outer membrane protein
ECK1398:JW1402:b1405; predicted outer membrane protein, C‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,473,9320GT11.9% 235.2 / 10.7 75pseudogene (6827/8622 nt)
pseudogene (1529/3324 nt)
ydbA
ydbA
ECK1398:JW5802+JW1402:b4492; predicted outer membrane protein
ECK1398:JW1402:b1405; predicted outer membrane protein, C‑ter fragment
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base T (5/4);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCGGCGATCAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCGGCACCGG  >  W3110S.gb/1473863‑1474001
                                                                     |                                                                     
ccGGCGATCAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAg                                                                       <  1:3912644/70‑1 (MQ=255)
 cGGCGATCAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTg                                                                     <  1:3703059/71‑1 (MQ=255)
       tCAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGa                                                               <  1:1874549/71‑1 (MQ=255)
        cAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGa                                                               >  1:1174596/1‑70 (MQ=255)
        cAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGa                                                               >  1:2263089/1‑70 (MQ=255)
        cAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGa                                                               >  1:1429402/1‑70 (MQ=255)
         aGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGaaa                                                             <  1:609317/71‑1 (MQ=255)
           aaTGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGaa                                                              >  1:1310041/1‑68 (MQ=255)
                 ttcttcTGGAACCATGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCa                                                         >  1:4298348/1‑67 (MQ=255)
                 ttcttcTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCa                                                         >  1:2451214/1‑67 (MQ=255)
                  tcttctGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCa                                                         <  1:2672359/66‑1 (MQ=255)
                   cttctGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAg                                                                       >  1:4432631/1‑51 (MQ=255)
                   cttctGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAg                                                                       >  1:4010210/1‑51 (MQ=255)
                   cttctGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCca                                                     <  1:4309581/69‑1 (MQ=255)
                       tGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGaa                                                              >  1:4008917/1‑56 (MQ=255)
                          aaCCATGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACtt                                            <  1:3624/71‑1 (MQ=255)
                          aaCCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACtt                                            <  1:4169913/71‑1 (MQ=255)
                            ccAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTa                                          <  1:4288761/71‑1 (MQ=255)
                              aGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCATTGCCg                                                                  >  1:3090656/1‑45 (MQ=255)
                              aGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCCACAGCAACTTTAtt                                        >  1:985746/1‑71 (MQ=255)
                              aGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTAtt                                        >  1:3361355/1‑71 (MQ=255)
                                gATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACtt                                            <  1:1559293/65‑1 (MQ=255)
                                gATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGc                                      <  1:34745/71‑1 (MQ=255)
                                 aTAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGcc                                     <  1:200935/71‑1 (MQ=255)
                                    aCGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACtt                                            <  1:761168/61‑1 (MQ=255)
                                      ggggTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCAcc                                                       >  1:848157/1‑48 (MQ=255)
                                       gggTCATCAACTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAAt                               >  1:3466501/1‑71 (MQ=255)
                                       gggTCATCAACTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAAt                               >  1:319067/1‑71 (MQ=255)
                                       gggTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAAt                               <  1:2701054/71‑1 (MQ=255)
                                          tcatcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGc                                      <  1:1959334/61‑1 (MQ=255)
                                          tcatcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCa                                   <  1:917331/64‑1 (MQ=255)
                                          tcatcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCa                                   <  1:1720668/64‑1 (MQ=255)
                                          tcatcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAg                                  <  1:2528368/65‑1 (MQ=255)
                                          tcatcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGAc                            >  1:142595/1‑71 (MQ=255)
                                            atcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACgg                          >  1:2283822/1‑71 (MQ=255)
                                             tcaACTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCcac                                                    <  1:3322763/44‑1 (MQ=255)
                                             tcaACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCa                                   <  1:1407306/61‑1 (MQ=255)
                                              caACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTa                                          >  1:1748812/1‑53 (MQ=255)
                                              caACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGa                             >  1:1123529/1‑66 (MQ=255)
                                               aaCTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCaa                                               >  1:3227372/1‑47 (MQ=255)
                                               aaCTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACtt                                            <  1:2279918/50‑1 (MQ=255)
                                               aaCTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACtt                                            <  1:3100475/50‑1 (MQ=255)
                                               aaCTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAg                                  <  1:835402/60‑1 (MQ=255)
                                               aaCTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGaa                                >  1:1655971/1‑62 (MQ=255)
                                               aaCTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGcat                       >  1:1238990/1‑71 (MQ=255)
                                                 cTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAg                                  <  1:1026306/58‑1 (MQ=255)
                                                 cTTCGATCACGACAGCGCCATTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAg                                  <  1:3295499/58‑1 (MQ=255)
                                                 cTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAAc                                              <  1:3019225/46‑1 (MQ=255)
                                                 cTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAAc                                              <  1:4086919/46‑1 (MQ=255)
                                                  ttCGATCACGACAGCGCCAGTGCCGTGAAAGTCAcc                                                       <  1:228414/36‑1 (MQ=255)
                                                  ttCGATCACGACAGCGCCAGTGCCGTGAAAGTCAcc                                                       <  1:3346143/36‑1 (MQ=255)
                                                   tCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGcatcat                    >  1:4561494/1‑70 (MQ=255)
                                                   tCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGcatcat                    >  1:3647128/1‑70 (MQ=255)
                                                   tCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGcatcat                    >  1:4407615/1‑70 (MQ=255)
                                                    cGATCACGACAGCGCCATTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATg                   >  1:2301074/1‑70 (MQ=255)
                                                    cGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTTTTGCCCAGAATGACGGCATCATg                   <  1:4032640/70‑1 (MQ=255)
                                                    cGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATg                   <  1:2002606/70‑1 (MQ=255)
                                                    cGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATg                   <  1:3796949/70‑1 (MQ=255)
                                                    cGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATg                   <  1:1230268/70‑1 (MQ=255)
                                                    cGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATg                   <  1:168964/70‑1 (MQ=255)
                                                      aTCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTAtt                                        <  1:3718815/47‑1 (MQ=255)
                                                       tCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAAt                               <  1:3778356/55‑1 (MQ=255)
                                                       tCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACa               <  1:1191780/71‑1 (MQ=255)
                                                        cACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAAc                >  1:1776277/1‑69 (MQ=255)
                                                        cACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACAt              >  1:2330555/1‑71 (MQ=255)
                                                        cACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACAt              >  1:4393471/1‑71 (MQ=255)
                                                        cACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACAt              >  1:1263196/1‑71 (MQ=255)
                                                         aCGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACttt                                           <  1:3103325/41‑1 (MQ=255)
                                                         aCGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATg                   <  1:294207/65‑1 (MQ=255)
                                                         aCGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATc             <  1:1042954/71‑1 (MQ=255)
                                                             cAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCg         >  1:4083851/1‑71 (MQ=255)
                                                              aGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCgg        <  1:1149150/71‑1 (MQ=255)
                                                              aGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCgg        >  1:928157/1‑71 (MQ=255)
                                                               gcgcCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCGGc       <  1:3031984/71‑1 (MQ=255)
                                                               gcgcCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCGGc       <  1:265471/71‑1 (MQ=255)
                                                                 gcCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAg           >  1:4073844/1‑65 (MQ=255)
                                                                   cAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACACCAGCGGCACCg   >  1:4133456/1‑71 (MQ=255)
                                                                    aGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCGGCACCgg  <  1:2346090/71‑1 (MQ=255)
                                                                     gTGCCGTGAAAGTCACCCCCAGCAACTTTATTGCCCAGaa                                >  1:1266769/1‑40 (MQ=255)
                                                                     gTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGaa                                >  1:88506/1‑40 (MQ=255)
                                                                     gTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGaa                                >  1:901270/1‑40 (MQ=255)
                                                                     gTGCCGTGAAAGTCACCCACAGCAACTTTATCGCCCAg                                  <  1:8556/38‑1 (MQ=255)
                                                                     |                                                                     
CCGGCGATCAGAATGCCTTCTTCTGGAACCAGGATAACGGGGTCATCAACTTCGATCACGACAGCGCCAGTGCCGTGAAAGTCACCCACAGCAACTTTATTGCCCAGAATGACGGCATCATGAACATCAGCGGCACCGG  >  W3110S.gb/1473863‑1474001

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: