Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,902,330 A→T 36.4% intergenic (+300/+1073) ygcE → / ← ygcF predicted kinase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,902,3300AT36.4% 51.1 / 57.2 55intergenic (+300/+1073)ygcE/ygcFpredicted kinase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base T (10/10);  total (28/27)
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.95e-01

AATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATT  >  W3110S.gb/2902267‑2902390
                                                               |                                                            
aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCgg                                                       >  1:1530304/1‑71 (MQ=255)
  ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGt                                                      >  1:1352776/1‑70 (MQ=255)
  ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGt                                                      >  1:3088489/1‑70 (MQ=255)
    ttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCgg                                                       >  1:4432188/1‑67 (MQ=255)
    ttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTAt                                                   >  1:1402591/1‑71 (MQ=255)
      cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATaga                                                            >  1:1960486/1‑60 (MQ=255)
        cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGt                                                      <  1:2669737/64‑1 (MQ=255)
        cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTaa                                               <  1:3578172/71‑1 (MQ=255)
        cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTaa                                               <  1:2859153/71‑1 (MQ=255)
           aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAAc                                            >  1:1786813/1‑71 (MQ=255)
           aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTaa                                               <  1:3727513/68‑1 (MQ=255)
           aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAc                                            >  1:895482/1‑71 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGtata                                                              >  1:501137/1‑47 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGtata                                                              >  1:698911/1‑47 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGtata                                                              <  1:269141/47‑1 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTc                                       <  1:2234410/70‑1 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCa                                      <  1:21216/71‑1 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCa                                      <  1:1302075/71‑1 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCa                                      <  1:2762660/71‑1 (MQ=255)
                 tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACACCGtata                                                              <  1:601445/47‑1 (MQ=255)
                     cTGGATATCTTTATTAATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTaa                                               <  1:2146992/58‑1 (MQ=255)
                        gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATaga                                                            >  1:3344365/1‑42 (MQ=255)
                        gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGc                                 <  1:1482786/69‑1 (MQ=255)
                        gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCt                               >  1:2886008/1‑71 (MQ=255)
                          tatCTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTaa                                               <  1:1387081/53‑1 (MQ=255)
                           atCTTTATTCATAAGAACTATTCATCAACAGCGTATtga                                                            >  1:3718111/1‑39 (MQ=255)
                            tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                            <  1:93296/70‑1 (MQ=255)
                            tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                            <  1:954002/70‑1 (MQ=255)
                            tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                            >  1:32129/1‑70 (MQ=255)
                             cTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCaaaa                          >  1:4518697/1‑71 (MQ=255)
                             cTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCaaaa                          >  1:2532206/1‑71 (MQ=255)
                             cTTTATTCATAAGAACTATTCATCAACAGCGTATTGAGGAGGTTATGTAAAACCACTCATTAGCCTCaaaa                          >  1:3032692/1‑71 (MQ=255)
                             cTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaaaa                          >  1:2764734/1‑71 (MQ=255)
                               ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTa                                                    >  1:3523175/1‑43 (MQ=255)
                                tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCAt                                     >  1:638494/1‑57 (MQ=255)
                                tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                            <  1:3146199/66‑1 (MQ=255)
                                  ttCATAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCAtt                                    >  1:1575210/1‑56 (MQ=255)
                                    cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTa                                                >  1:1253695/1‑42 (MQ=255)
                                     aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACa                  <  1:4314271/71‑1 (MQ=255)
                                      tAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTaaa                                              >  1:1732440/1‑42 (MQ=255)
                                      tAAGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACa                  <  1:2991657/70‑1 (MQ=255)
                                        aGAACTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAAc               <  1:2508392/71‑1 (MQ=255)
                                        aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaaa                <  1:631564/70‑1 (MQ=255)
                                           aCTATTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTc            <  1:2875807/71‑1 (MQ=255)
                                              aTTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaaa                >  1:3480796/1‑64 (MQ=255)
                                              aTTCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTc            >  1:4328188/1‑68 (MQ=255)
                                               ttCATCAACAGCGTATTGAGGCGGTTATGTAAAACCACTCa                                      <  1:4219699/41‑1 (MQ=255)
                                               ttCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTc            >  1:155281/1‑67 (MQ=255)
                                                tcatcaACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACaaa                      <  1:1622163/56‑1 (MQ=255)
                                                tcatcaACAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTcc           >  1:4302472/1‑67 (MQ=255)
                                                tcatcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGAt       >  1:2085051/1‑71 (MQ=255)
                                                 catcaACAGCGTATATAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGAtt      >  1:3998782/1‑71 (MQ=255)
                                                  atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGAttt     <  1:1174674/71‑1 (MQ=255)
                                                   tcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                            <  1:1952979/47‑1 (MQ=255)
                                                   tcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGa                    <  1:3057458/55‑1 (MQ=255)
                                                   tcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaa                 <  1:537629/58‑1 (MQ=255)
                                                     aaCAGCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATttat   >  1:3857647/1‑70 (MQ=255)
                                                     aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATttatt  <  1:4232928/71‑1 (MQ=255)
                                                     aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATttat   >  1:4534203/1‑70 (MQ=255)
                                                     aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCAAGAt       >  1:3852384/1‑66 (MQ=255)
                                                     aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCAAGATttat   >  1:4396570/1‑70 (MQ=255)
                                                         gCGTATTGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATttat   <  1:1090097/66‑1 (MQ=255)
                                                               |                                                            
AATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATT  >  W3110S.gb/2902267‑2902390

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: