Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 531,446 C→T 56.1% intergenic (‑1/‑229) ybbS ← / → allA transcriptional activator of the allD operon/ureidoglycolate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb531,4460CT56.1% 9.0 / 52.7 41intergenic (‑1/‑229)ybbS/allAtranscriptional activator of the allD operon/ureidoglycolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (6/12);  new base T (5/18);  total (11/30)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTA  >  W3110S.gb/531380‑531514
                                                                  |                                                                    
cGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCa                                                                  <  1:1427187/71‑1 (MQ=255)
  cTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAgg                                                                <  1:647063/71‑1 (MQ=255)
  cTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAgg                                                                <  1:1878352/71‑1 (MQ=255)
  cTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAgg                                                                <  1:3543929/71‑1 (MQ=255)
  cTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAgg                                                                <  1:2317764/71‑1 (MQ=255)
     ttGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGt                                                               <  1:1160998/69‑1 (MQ=255)
     ttGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGt                                                               <  1:2088048/69‑1 (MQ=255)
                ccTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTa                                                    <  1:468384/69‑1 (MQ=255)
                     ttCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCg                                              >  1:765590/1‑70 (MQ=255)
                     ttCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCg                                              >  1:1749263/1‑70 (MQ=255)
                       cAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCa                                            <  1:2900811/70‑1 (MQ=255)
                        aGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCaa                                           <  1:1239005/70‑1 (MQ=255)
                        aGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCaa                                           <  1:1232007/70‑1 (MQ=255)
                        aGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCaa                                           <  1:2832313/70‑1 (MQ=255)
                             ccGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCa                                                                  <  1:1685002/42‑1 (MQ=255)
                             ccGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCa                                                                  <  1:2965180/42‑1 (MQ=255)
                              cGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACAtttt                                                                    >  1:388959/1‑39 (MQ=255)
                              cGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACAtttt                                                                    >  1:111716/1‑39 (MQ=255)
                                cAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGa                                  <  1:2977052/71‑1 (MQ=255)
                                 aaTGAAAGTCCGCAAGGTTTCTGGATCGAACAtttt                                                                    <  1:2532336/36‑1 (MQ=255)
                                 aaTGAAAGTCCGCAAGGTTTCTGGATCGAACAtttt                                                                    <  1:1700478/36‑1 (MQ=255)
                                 aaTAAAAGTCCGCAAGGTTTCTGGATCGAACAtttt                                                                    <  1:1363139/36‑1 (MQ=37)
                                   tGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATcc                                                           <  1:2604631/43‑1 (MQ=255)
                                   tGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATcc                                                           <  1:115119/43‑1 (MQ=255)
                                        gTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:1583186/71‑1 (MQ=255)
                                         tCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTAt                         <  1:2736954/71‑1 (MQ=255)
                                          ccGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATc                        >  1:2394048/1‑71 (MQ=255)
                                           cGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATc                        <  1:124420/70‑1 (MQ=255)
                                               aGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCctct                    >  1:2873133/1‑70 (MQ=255)
                                                  tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:2818720/61‑1 (MQ=255)
                                                  tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:2780995/61‑1 (MQ=255)
                                                  tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:2502292/61‑1 (MQ=255)
                                                  tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:395393/61‑1 (MQ=255)
                                                  tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAtt                <  1:166331/71‑1 (MQ=255)
                                                  tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAATTGCGTGATTTTCTTATCCTCTGAtt                <  1:1432415/71‑1 (MQ=255)
                                                         aTCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGt         <  1:1742124/71‑1 (MQ=255)
                                                            gAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtt      >  1:2785666/1‑71 (MQ=255)
                                                            gAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtt      >  1:225090/1‑71 (MQ=255)
                                                               cATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:3123540/48‑1 (MQ=255)
                                                               cATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:331030/48‑1 (MQ=255)
                                                               cATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtttt    >  1:1847084/1‑70 (MQ=255)
                                                               cATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtttt    >  1:2977800/1‑70 (MQ=255)
                                                               cATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtttt    >  1:2030415/1‑70 (MQ=255)
                                                               cATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtttt    >  1:2484459/1‑70 (MQ=255)
                                                                aTCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTa  <  1:3644304/71‑1 (MQ=255)
                                                                 tCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTa  >  1:1986909/1‑70 (MQ=255)
                                                                  |                                                                    
CGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTA  >  W3110S.gb/531380‑531514

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