Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,885,984 C→T 75.0% intergenic (‑109/+250) ygcB ← / ← cysH conserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,885,9840CT75.0% 121.1 / 42.5 68intergenic (‑109/+250)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base T (20/31);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCC  >  W3110S.gb/2885918‑2886051
                                                                  |                                                                   
gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAaag                                                                 <  1:557246/71‑1 (MQ=255)
    tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAaagaa                                                               >  1:1792885/1‑69 (MQ=255)
    tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAaagaa                                                               >  1:3038587/1‑69 (MQ=255)
    tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAaagaa                                                               >  1:3030359/1‑69 (MQ=255)
    tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAaag                                                                 <  1:2623572/67‑1 (MQ=255)
    tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAaag                                                                 <  1:1399358/67‑1 (MQ=255)
    tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAaag                                                                 <  1:2682088/67‑1 (MQ=255)
      tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                             <  1:90202/69‑1 (MQ=255)
      tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                             <  1:178427/69‑1 (MQ=255)
      tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                             <  1:1675540/69‑1 (MQ=255)
      tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                             <  1:1528495/69‑1 (MQ=255)
      tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                             <  1:3082496/69‑1 (MQ=255)
       ttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGAc                                                          >  1:3500831/1‑71 (MQ=255)
        tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAgaga                                                           >  1:288236/1‑69 (MQ=255)
        tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAgaga                                                           >  1:3421998/1‑69 (MQ=255)
        tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAgaga                                                           >  1:1143549/1‑69 (MQ=255)
        tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAgaga                                                           >  1:2263479/1‑69 (MQ=255)
        tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGAc                                                          >  1:282455/1‑70 (MQ=255)
         ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         >  1:3278749/1‑70 (MQ=255)
         ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         >  1:3274445/1‑70 (MQ=255)
         ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         >  1:1435604/1‑70 (MQ=255)
         ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         >  1:520186/1‑70 (MQ=255)
         ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         >  1:2523487/1‑70 (MQ=255)
         ttCGATGATGCTCATAAATATCGAATTACAGTGATATACACATGTATTGCCTCCCTGTAAAGAAGAGACt                                                         <  1:791748/70‑1 (MQ=255)
               gctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTc                                                  >  1:963649/1‑71 (MQ=255)
               gctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTc                                                  >  1:369568/1‑71 (MQ=255)
                      aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGt                                           <  1:267302/71‑1 (MQ=255)
                      aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTg                                            <  1:135191/70‑1 (MQ=255)
                      aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTg                                            <  1:3652317/70‑1 (MQ=255)
                      aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTg                                            <  1:1562826/70‑1 (MQ=255)
                      aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTg                                            <  1:2958557/70‑1 (MQ=255)
                      aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTg                                            <  1:2331466/70‑1 (MQ=255)
                           tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                                      <  1:1008363/71‑1 (MQ=255)
                             tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                                      <  1:2484186/69‑1 (MQ=255)
                             tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                                      <  1:382434/69‑1 (MQ=255)
                             tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                                      <  1:421623/69‑1 (MQ=255)
                                  ttACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         <  1:186835/45‑1 (MQ=255)
                                  ttACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         <  1:2449449/45‑1 (MQ=255)
                                  ttACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACt                                                         <  1:3454372/45‑1 (MQ=255)
                                    aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTaa                                                       >  1:1356269/1‑45 (MQ=255)
                                    aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTaa                                                       >  1:1417332/1‑45 (MQ=255)
                                    aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTaa                                                       >  1:3127122/1‑45 (MQ=255)
                                    aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg                              <  1:3639024/70‑1 (MQ=255)
                                    aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg                              <  1:931901/70‑1 (MQ=255)
                                    aCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                             >  1:2475800/1‑71 (MQ=255)
                                         gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTg                                            <  1:2886130/51‑1 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCg                                  <  1:3382354/57‑1 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCg                                  <  1:1192199/57‑1 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCg                                  <  1:2222446/57‑1 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCg                                  <  1:2314480/57‑1 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCg                                  <  1:1712088/57‑1 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                    >  1:1986743/1‑71 (MQ=255)
                                             taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                    >  1:2200779/1‑71 (MQ=255)
                                               cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg                                            >  1:551332/1‑45 (MQ=255)
                                                   tGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGa                                 <  1:2048215/52‑1 (MQ=255)
                                                   tGTCCTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGa                                 <  1:3576247/52‑1 (MQ=255)
                                                     tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGaa                                               <  1:2557891/36‑1 (MQ=255)
                                                     tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGaa                                               <  1:2263004/36‑1 (MQ=255)
                                                     tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGaa                                               <  1:1941030/36‑1 (MQ=255)
                                                     tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa               >  1:3264576/1‑68 (MQ=255)
                                                     tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa               >  1:1926121/1‑68 (MQ=255)
                                                       ttGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa           >  1:426919/1‑70 (MQ=255)
                                                         gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa         >  1:3065051/1‑70 (MQ=255)
                                                         gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa           <  1:288769/68‑1 (MQ=255)
                                                         gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa           <  1:1715416/68‑1 (MQ=255)
                                                         gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa           <  1:59024/68‑1 (MQ=255)
                                                         gCCTCCCTACAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa         <  1:628792/70‑1 (MQ=255)
                                                               cTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc  <  1:2256509/71‑1 (MQ=255)
                                                                  |                                                                   
GGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCC  >  W3110S.gb/2885918‑2886051

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: