Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,045,238 A→G 25.0% intergenic (+203/+367) yodB → / ← serU predicted cytochrome/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,045,2380AG25.0% 53.7 / 17.1 32intergenic (+203/+367)yodB/serUpredicted cytochrome/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (15/9);  new base G (2/6);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAAAAATACAGGTTAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCTTTAATCTGA  >  W3110S.gb/2045170‑2045300
                                                                    |                                                              
aTCAAAAATACAGGTTAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCaaa                                                               >  1:2155773/1‑70 (MQ=255)
 tCAAAAATACAGGTTAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCaaaa                                                              >  1:841730/1‑70 (MQ=255)
   aaaaaTACAGGTTAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAGAAAAg                                                           >  1:625856/1‑71 (MQ=255)
        tACAGGTTAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCg                                                      <  1:258110/71‑1 (MQ=255)
               tAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTg                                                        >  1:895879/1‑62 (MQ=255)
               tAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTg                                                        >  1:519179/1‑62 (MQ=255)
               tAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTc                                                >  1:3704008/1‑70 (MQ=255)
                       gTATTTAAACCTTTTTTACTTGTTAGTATGATATAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAActc                                        <  1:2731130/70‑1 (MQ=255)
                       gTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAActc                                        <  1:2556849/70‑1 (MQ=255)
                       gTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAAct                                         <  1:2912965/69‑1 (MQ=255)
                        tATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAActc                                        <  1:4084543/69‑1 (MQ=255)
                          tttAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAGAAAAGCTGCGAAATTcaaca                                            <  1:3930516/63‑1 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAGAAAAGCTGCGAAATTCAACAACTCTTGGa                                   <  1:1883953/70‑1 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTgg                                    >  1:2793777/1‑69 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTgg                                    >  1:3405190/1‑69 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   >  1:427230/1‑70 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   <  1:1353745/70‑1 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   <  1:3687299/70‑1 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   >  1:3454354/1‑70 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   >  1:3276409/1‑70 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   >  1:691262/1‑70 (MQ=255)
                            tAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGa                                   <  1:2348580/70‑1 (MQ=255)
                                ccTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAGAAAAGCTGCGAAATTCAACAACTCTTGGAAATTa                              <  1:1668078/71‑1 (MQ=255)
                                   tttttACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATg                            >  1:3635476/1‑70 (MQ=255)
                                        aCTTGTTAGTATGATACAAAGGCTTTCAGAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACt                        <  1:2134052/69‑1 (MQ=255)
                                                   tGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTa                              >  1:144856/1‑52 (MQ=255)
                                                   tGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTa                              >  1:2283630/1‑52 (MQ=255)
                                                   tGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCt           >  1:297966/1‑71 (MQ=255)
                                                      tACAAAGGCTTTCAGAAAAGCTGCTAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCttt         <  1:507144/70‑1 (MQ=255)
                                                      tACAAAGGCTTTCAGAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCttt         <  1:1338095/70‑1 (MQ=255)
                                                      tACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCTTTa        <  1:245345/71‑1 (MQ=255)
                                                          aaGGCTTTCAGAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCcg                >  1:3103783/1‑59 (MQ=255)
                                                          aaGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCTTTAATCt    >  1:1073292/1‑71 (MQ=255)
                                                            ggCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCTTTAATCTGa  >  1:1401416/1‑71 (MQ=255)
                                                                    |                                                              
ATCAAAAATACAGGTTAATTTTGGTATTTAAACCTTTTTTACTTGTTAGTATGATACAAAGGCTTTCAAAAAAGCTGCGAAATTCAACAACTCTTGGAAATTATGGACTTTGTCCCGCGGCTTTAATCTGA  >  W3110S.gb/2045170‑2045300

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