Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 531,446 C→T 62.7% 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,4460CT62.7% 37.9 / 54.3 51intergenic (‑1/‑229)ybbS/allAtranscriptional activator of the allD operon/ureidoglycolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (10/9);  new base T (21/11);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTA  >  W3110S.gb/531384‑531514
                                                              |                                                                    
tttGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAgg                                                                <  1:1452909/69‑1 (MQ=255)
                        aCCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGtt                                        >  1:2774060/1‑69 (MQ=255)
                        aCCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGtt                                        >  1:1308277/1‑69 (MQ=255)
                           gCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTc                                                                   >  1:3395906/1‑39 (MQ=255)
                           gCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTc                                                                   >  1:2640501/1‑39 (MQ=255)
                           gCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTc                                                                   >  1:2757911/1‑39 (MQ=255)
                           gCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATTTTc                                                                   >  1:3542192/1‑39 (MQ=255)
                             aaTGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTg                                   <  1:586106/69‑1 (MQ=255)
                                gAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGAtt                                >  1:2841748/1‑69 (MQ=255)
                                gAAAGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGAtt                                >  1:2200148/1‑69 (MQ=255)
                                gAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCtttt                                                     >  1:3225897/1‑48 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTa                                                              <  1:403367/38‑1 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTa                                                              <  1:1968583/38‑1 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTa                                                              <  1:575408/38‑1 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTa                                                              <  1:598080/38‑1 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGAttt                               >  1:1452267/1‑69 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGAttt                               >  1:193745/1‑69 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGt                                    >  1:1113591/1‑64 (MQ=255)
                                 aaaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGt                                    >  1:3387789/1‑64 (MQ=255)
                                   aGTCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGAtttt                              >  1:2150342/1‑68 (MQ=255)
                                   aGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGAtttt                              <  1:3170467/68‑1 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           >  1:66270/1‑69 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           >  1:1841714/1‑69 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           >  1:342578/1‑69 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           >  1:2016518/1‑69 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCt                            >  1:130660/1‑68 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATACGCAAGTTGCGTGATTTTCtt                           >  1:2506928/1‑69 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           <  1:2657077/69‑1 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           <  1:748779/69‑1 (MQ=255)
                                     tCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           <  1:1570314/69‑1 (MQ=255)
                                      ccGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           <  1:2442170/68‑1 (MQ=255)
                                      ccGCAAGGTTTCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                          <  1:3321938/69‑1 (MQ=255)
                                              tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCt                            >  1:2692193/1‑59 (MQ=255)
                                              tttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:374347/1‑69 (MQ=255)
                                              tttCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:2755203/1‑69 (MQ=255)
                                               ttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:3678457/1‑68 (MQ=255)
                                               ttCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:262497/1‑68 (MQ=255)
                                               ttCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:964025/1‑68 (MQ=255)
                                               ttCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:1241859/1‑68 (MQ=255)
                                               ttCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:2961711/1‑68 (MQ=255)
                                               ttCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                  >  1:2914718/1‑68 (MQ=255)
                                                tCTGGATCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAtt                >  1:1160647/1‑69 (MQ=255)
                                                  tGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGa                                  <  1:3221539/49‑1 (MQ=255)
                                                      tCGAACATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAg          >  1:3295529/1‑69 (MQ=255)
                                                      tCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAg          <  1:3344440/69‑1 (MQ=255)
                                                      tCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAg          <  1:3134796/69‑1 (MQ=255)
                                                       cGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGt         <  1:407332/69‑1 (MQ=255)
                                                         aaCATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCtt                           <  1:1422179/49‑1 (MQ=255)
                                                          aCATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAt       <  1:2485568/68‑1 (MQ=255)
                                                          aCATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAt       <  1:2489940/68‑1 (MQ=255)
                                                           cATTTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAt       <  1:2006224/67‑1 (MQ=255)
                                                              cTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTa  >  1:1036144/1‑69 (MQ=255)
                                                              |                                                                    
TTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTA  >  W3110S.gb/531384‑531514

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: