Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 571,256 A→G 21.4% T194A (ACA→GCA)  ybcM → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb571,2560AG21.4% 141.7 / 35.7 70T194A (ACA→GCA) ybcMpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/24);  new base G (5/10);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAA  >  W3110S.gb/571187‑571326
                                                                     |                                                                      
ctgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCta                                                                        <  1:1430326/70‑1 (MQ=255)
ctgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCta                                                                        <  1:3571743/70‑1 (MQ=255)
   ctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtgcat                                                                     <  1:190752/70‑1 (MQ=255)
     aTCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacatta                                                                   >  1:958767/1‑70 (MQ=255)
      tCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacattac                                                                  <  1:2229949/70‑1 (MQ=255)
      tCCGAACACGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacattac                                                                  <  1:2732700/70‑1 (MQ=255)
        cGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATc                                                               >  1:1458243/1‑71 (MQ=255)
           aCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtgca                                                                      >  1:2341005/1‑61 (MQ=255)
           aCCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTGCATTACATCGcc                                                            >  1:1290360/1‑71 (MQ=255)
               ggCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCt                                                        <  1:2052440/71‑1 (MQ=255)
                gCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacat                                                                     <  1:2453276/57‑1 (MQ=255)
                 cacaAAATGGACAGCAAACTAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCg                                                              <  1:724799/63‑1 (MQ=255)
                  acaAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacattac                                                                  <  1:4280860/58‑1 (MQ=255)
                  acaAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacat                                                                     <  1:2892081/55‑1 (MQ=255)
                  acaAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAg                                                      >  1:343371/1‑70 (MQ=255)
                   caAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTGCATTACATc                                                               <  1:1177075/60‑1 (MQ=255)
                   caAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGcc                                                            >  1:3432677/1‑63 (MQ=255)
                     aaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAg                                                  >  1:2371494/1‑71 (MQ=255)
                     aaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTa                                                       >  1:317659/1‑66 (MQ=255)
                     aaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCaa                                                   >  1:3477542/1‑70 (MQ=255)
                     aaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAg                                                  >  1:1412977/1‑71 (MQ=255)
                      aaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGc                                                 <  1:4426142/71‑1 (MQ=255)
                      aaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGc                                                 <  1:844771/71‑1 (MQ=255)
                      aaTGGACAGCAAACAAGGTTGCTCGATATCTATACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGc                                                 >  1:571871/1‑71 (MQ=255)
                       aTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCg                                                <  1:4511547/71‑1 (MQ=255)
                       aTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCg                                                <  1:2815675/71‑1 (MQ=255)
                       aTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCg                                                <  1:1925089/71‑1 (MQ=255)
                          gACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacat                                                                     <  1:453884/47‑1 (MQ=255)
                          gACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAg                                                      <  1:1952635/62‑1 (MQ=255)
                          gACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                             <  1:2292861/71‑1 (MQ=255)
                           aCAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacatt                                                                    <  1:489821/47‑1 (MQ=255)
                              gCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGcc                                                            <  1:1384268/52‑1 (MQ=255)
                              gCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAggggg                                          >  1:1240268/1‑70 (MQ=255)
                               cAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtgcattaca                                                                 <  1:1300094/46‑1 (MQ=255)
                               cAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacatta                                                                   >  1:1596114/1‑44 (MQ=255)
                               cAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCtacatta                                                                   >  1:2922384/1‑44 (MQ=255)
                               cAAACAAGGGTGCTCGATATCTCTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGc                                                 <  1:784762/62‑1 (MQ=255)
                                 aaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATc                                                               >  1:254766/1‑46 (MQ=255)
                                 aaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                             >  1:3937654/1‑64 (MQ=255)
                                 aaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAg                                      >  1:169874/1‑71 (MQ=255)
                                  aCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGc                                                 <  1:4400593/59‑1 (MQ=255)
                                  aCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCg                                                <  1:1017532/60‑1 (MQ=255)
                                  aCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGt                                     <  1:1561389/71‑1 (MQ=255)
                                   cAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAg                                                  >  1:3643905/1‑57 (MQ=255)
                                   cAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAg                                                  >  1:660547/1‑57 (MQ=255)
                                   cAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAg                                                  >  1:665964/1‑57 (MQ=255)
                                   cAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGtt                                    <  1:4065197/71‑1 (MQ=255)
                                   cAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGtt                                    >  1:1455605/1‑71 (MQ=255)
                                    aaGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCg                                                >  1:1214813/1‑58 (MQ=255)
                                     aGGTTGCTCGATATCTCTACATTTCTGTTTCtgcattaca                                                                 >  1:574003/1‑40 (MQ=255)
                                     aGGTTGCTCGATATCTCTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTc                                  >  1:1003719/1‑71 (MQ=255)
                                        ttGCTCGATATCTCTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCaa                               <  1:1572608/71‑1 (MQ=255)
                                        ttGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGcc                                                            >  1:3443123/1‑42 (MQ=255)
                                         tGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCa                                >  1:4159031/1‑69 (MQ=255)
                                          gCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCaaa                              <  1:1834328/70‑1 (MQ=255)
                                           cTCGATATCTCTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAg                             <  1:2022369/70‑1 (MQ=255)
                                             cGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGtt                           >  1:1920804/1‑69 (MQ=255)
                                                tatCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAg                                              <  1:1559010/48‑1 (MQ=255)
                                                tatCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACt                       >  1:1872640/1‑71 (MQ=255)
                                                 atCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGc                                                     >  1:4120693/1‑40 (MQ=255)
                                                  tctcTACATTTCTGTTTCTGCATTACATCGCCGTCTAGc                                                     <  1:4257072/39‑1 (MQ=255)
                                                  tctcTACATTTCTGTTTCTGCATTACATCGCCGTCTAGc                                                     <  1:1930868/39‑1 (MQ=255)
                                                  tctcTACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTg                      <  1:4503973/70‑1 (MQ=255)
                                                  tctcTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCg                                                <  1:2373913/44‑1 (MQ=255)
                                                  tctcTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTg                      <  1:956742/70‑1 (MQ=255)
                                                   ctctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGa                    >  1:971645/1‑71 (MQ=255)
                                                    tctACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCGa                                               >  1:2864823/1‑43 (MQ=255)
                                                    tctACATTTCTGTTTCTACATTACATCGCCGTCTAGCa                                                    <  1:808391/38‑1 (MQ=255)
                                                    tctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGAc                   >  1:3700129/1‑71 (MQ=255)
                                                      tACATTTCTGTTTCTGCATTACATCGCCGTCTAGCAAGCGAg                                              <  1:4077149/42‑1 (MQ=255)
                                                      tACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAg                                                  >  1:3260086/1‑38 (MQ=255)
                                                      tACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCg                                                <  1:2377175/40‑1 (MQ=255)
                                                       aCATTTCTGTTTCTACATTACATCGCCGTCTAGAAAGCGAGGGGGTAAGTTTCCAAAGta                           >  1:1256083/1‑60 (MQ=255)
                                                        cATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAg                                              >  1:1714413/1‑40 (MQ=255)
                                                        cATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGa                 <  1:4370908/69‑1 (MQ=255)
                                                          tttctgtttctgCATTACATCGCCGTCTAGCAAGCg                                                <  1:1328201/36‑1 (MQ=255)
                                                          tttctgtttctgCATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGtt                                    >  1:3478737/1‑48 (MQ=255)
                                                          tttCTGTTTCTACATTACATCGCCGTCTAGCAAGc                                                 <  1:3861568/35‑1 (MQ=255)
                                                          tttCTGTTTCTACATTACATCGCCGTCTAGCAAGc                                                 <  1:4548070/35‑1 (MQ=255)
                                                          tttCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGAtgt              >  1:1466129/1‑70 (MQ=255)
                                                           ttCTGTTTCTGCATTACATCGCCGTCTAGCAAGCg                                                <  1:1413752/35‑1 (MQ=255)
                                                           ttCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGta                           <  1:3054718/56‑1 (MQ=255)
                                                            tCTGTTTCTACATTACATCGCCGTCTAGCAAGCGa                                               >  1:183141/1‑35 (MQ=255)
                                                            tCTGTTTCTACATTACATCGCCGTCTAGCAAGCGa                                               >  1:3187085/1‑35 (MQ=255)
                                                            tCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAg           >  1:3209842/1‑71 (MQ=255)
                                                             cTGTTTCTGCATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATAc                        >  1:3600105/1‑57 (MQ=255)
                                                             cTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAgg          >  1:1985428/1‑71 (MQ=255)
                                                                  tCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTa       <  1:3218828/69‑1 (MQ=255)
                                                                     acattacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaataa  >  1:1852721/1‑71 (MQ=255)
                                                                     |                                                                      
CTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAA  >  W3110S.gb/571187‑571326

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