Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,885,984 C→T 100% 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,9840CT85.1% 167.6 / 21.7 67intergenic (‑109/+250)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (10/0);  new base T (32/25);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCT  >  W3110S.gb/2885924‑2886052
                                                            |                                                                    
tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                              >  1:1168625/1‑69 (MQ=255)
tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAagaaga                                                              >  1:3587222/1‑69 (MQ=255)
  tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTaaa                                                                   >  1:2874149/1‑62 (MQ=255)
  tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTaa                                                                    >  1:1681853/1‑61 (MQ=255)
  tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAagag                                                             >  1:1112659/1‑68 (MQ=255)
            gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCg                                                  <  1:3581956/69‑1 (MQ=255)
               cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACCTTCGaa                                                <  1:480809/68‑1 (MQ=255)
               cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGaa                                                <  1:815402/68‑1 (MQ=255)
               cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGa                                                 <  1:2979922/67‑1 (MQ=255)
               cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGa                                                 <  1:1870292/67‑1 (MQ=255)
               cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAAt                                               >  1:941593/1‑69 (MQ=255)
                   aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAaagaa                                                                >  1:811051/1‑48 (MQ=255)
                   aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAaagaa                                                                >  1:2139017/1‑48 (MQ=255)
                   aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTa                                                         >  1:1130231/1‑55 (MQ=255)
                   aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTa                                                         >  1:2007342/1‑55 (MQ=255)
                   aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGtt                                           >  1:1690541/1‑69 (MQ=255)
                     tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTg                                         <  1:3191690/69‑1 (MQ=255)
                      atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTg                                         <  1:1941771/68‑1 (MQ=255)
                       tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTg                                         <  1:1245909/67‑1 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCg                                                  <  1:2752518/56‑1 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          >  1:3147197/1‑64 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          >  1:3132458/1‑64 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          <  1:1447047/64‑1 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          <  1:2770111/64‑1 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          >  1:2600523/1‑64 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          <  1:236173/64‑1 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          >  1:1587138/1‑64 (MQ=255)
                         gAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGttt                                          >  1:1603447/1‑64 (MQ=255)
                            ttACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGa                                  >  1:3490432/1‑69 (MQ=255)
                             tACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGaa                                 >  1:3720176/1‑69 (MQ=255)
                             tACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGaa                                 >  1:796997/1‑69 (MQ=255)
                                aGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg                               <  1:1395813/68‑1 (MQ=255)
                                aGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg                               <  1:1426293/68‑1 (MQ=255)
                                aGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg                               <  1:2918179/68‑1 (MQ=255)
                                 gTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg                               <  1:1869044/67‑1 (MQ=255)
                                    atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGc                          >  1:716313/1‑69 (MQ=255)
                                    atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGc                          >  1:210720/1‑69 (MQ=255)
                                      ataCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGc                          >  1:3351379/1‑67 (MQ=255)
                                      ataCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGc                          >  1:3202422/1‑67 (MQ=255)
                                       taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGATTTGTTTGGGTTCGAACGCTGGCCTc                       >  1:1944121/1‑69 (MQ=255)
                                       taCACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTc                       >  1:2556216/1‑69 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg                                             >  1:1265921/1‑45 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     >  1:1136463/1‑69 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     <  1:3668820/69‑1 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     >  1:1564156/1‑69 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     >  1:2191455/1‑69 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     <  1:2701896/69‑1 (MQ=255)
                                         cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     >  1:273226/1‑69 (MQ=255)
                                          acaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAg                     <  1:1699100/68‑1 (MQ=255)
                                          acaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGACTCAg                     <  1:1616138/68‑1 (MQ=255)
                                             tGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGt                                            >  1:3656294/1‑42 (MQ=255)
                                             tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGtt                  >  1:2581193/1‑68 (MQ=255)
                                             tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGtt                  >  1:3168123/1‑68 (MQ=255)
                                               tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCAAACGCTGGCCTCAGGtt                  >  1:657709/1‑66 (MQ=255)
                                               tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCAAACGCTGGCCTCAGGtt                  >  1:1890462/1‑66 (MQ=255)
                                               tCTTGCATCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGtt                  >  1:323590/1‑66 (MQ=255)
                                                cTTGCCTCCCTGTAAAGAAGAGACTAAGATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt               >  1:692951/1‑68 (MQ=255)
                                                 ttGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATa              >  1:3317088/1‑68 (MQ=255)
                                                    ccTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                <  1:2554831/63‑1 (MQ=255)
                                                    ccTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                <  1:3445001/63‑1 (MQ=255)
                                                    ccTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTGCGAACGCTGGCCTCAGGTTGa                <  1:3385129/63‑1 (MQ=255)
                                                     cTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa            <  1:12442/66‑1 (MQ=255)
                                                     cTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa            <  1:2727646/66‑1 (MQ=255)
                                                          tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGc    >  1:845981/1‑69 (MQ=255)
                                                           gCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa            >  1:2872421/1‑60 (MQ=255)
                                                           gCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa            >  1:3189592/1‑60 (MQ=255)
                                                           gCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa            >  1:2544159/1‑60 (MQ=255)
                                                            cAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCt  >  1:589635/1‑69 (MQ=255)
                                                            cAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCt  >  1:854416/1‑69 (MQ=255)
                                                            |                                                                    
TTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCT  >  W3110S.gb/2885924‑2886052

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: