Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,404,829 C→G 31.3% intergenic (+516/‑462) yhdH → / → accB predicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,404,8290CG31.3% 81.6 / 50.2 67intergenic (+516/‑462)yhdH/accBpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding/acetyl CoA carboxylase, BCCP subunit
Reads supporting (aligned to +/- strand):  ref base C (26/20);  new base G (10/11);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCGATGATGGCGCTAA  >  W3110S.gb/3404761‑3404895
                                                                    |                                                                  
cTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCtt                                                                  >  1:3284454/1‑71 (MQ=255)
   gTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCg                                                                >  1:779181/1‑70 (MQ=255)
    tCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGc                                                               <  1:2519527/70‑1 (MQ=255)
      gggATCCTGGCATCCCTACATTTTTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAgg                                                            <  1:284931/69‑1 (MQ=255)
       cgATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAgg                                                            <  1:541849/70‑1 (MQ=255)
         aTCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct                                                         <  1:960671/71‑1 (MQ=255)
         aTCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggct                                                         <  1:407808/71‑1 (MQ=255)
            cTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGggctgg                                                       <  1:3950550/70‑1 (MQ=255)
            cTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGggctgg                                                       <  1:3033417/70‑1 (MQ=255)
              ggCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgctt                                                    >  1:110712/1‑71 (MQ=255)
              ggCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgctt                                                    >  1:2342106/1‑71 (MQ=255)
               gCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttt                                                   <  1:2612832/71‑1 (MQ=255)
                    ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagc                                              <  1:2450972/71‑1 (MQ=255)
                    ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagc                                              >  1:2947892/1‑71 (MQ=255)
                    ccTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagc                                              >  1:2761726/1‑71 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagctttt                                          >  1:1134486/1‑71 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagctttt                                          >  1:1763465/1‑71 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagctt                                            <  1:2929794/69‑1 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAg                                                             >  1:630452/1‑52 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                                          >  1:1025237/1‑71 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                                          >  1:2814651/1‑71 (MQ=255)
                        cATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                                          >  1:533322/1‑71 (MQ=255)
                         attattTGTGGTCTGTGAAGAATATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagctttt                                          >  1:553761/1‑70 (MQ=255)
                          ttattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttttca                                        >  1:2700442/1‑71 (MQ=255)
                           tattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGggctg                                                        >  1:1449926/1‑54 (MQ=255)
                           tattTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagcttttca                                        >  1:406749/1‑70 (MQ=255)
                            attTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCAc                                       <  1:2461004/70‑1 (MQ=255)
                               tgtgGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCtt                                    >  1:326326/1‑70 (MQ=255)
                                   gTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGt                                <  1:3941537/70‑1 (MQ=255)
                                   gTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGt                                <  1:679963/70‑1 (MQ=255)
                                    tCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagcttttca                                        <  1:3128216/61‑1 (MQ=255)
                                    tCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGt                                <  1:3724326/69‑1 (MQ=255)
                                      tgtgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgcttttcagct                                             >  1:2441564/1‑54 (MQ=255)
                                      tgtgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGt                                >  1:432170/1‑67 (MQ=255)
                                      tgtgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat                             >  1:2938370/1‑70 (MQ=255)
                                        tgAAGATTATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat                             >  1:1186297/1‑68 (MQ=255)
                                        tgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtata                            <  1:1660773/69‑1 (MQ=255)
                                        tgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtata                            <  1:1533662/69‑1 (MQ=255)
                                        tgAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat                             <  1:702399/68‑1 (MQ=255)
                                           aGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctggct                                                     <  1:1328749/41‑1 (MQ=255)
                                              ttATCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagct                                             <  1:117550/46‑1 (MQ=255)
                                              ttATCTCATTGCAGCCCCTCATCTTCGCAGGGCAGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGtt                     >  1:2351032/1‑70 (MQ=255)
                                              ttATCTCATTGCAGCCCCTCATCTTCGCAGGGCAGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGt                      >  1:1273267/1‑69 (MQ=255)
                                                aTCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagcttt                                           >  1:1634419/1‑46 (MQ=255)
                                                aTCTCATTGCAGCCCCTCATGTTCGCAGGGCTGgcttttcagcttt                                           >  1:3894893/1‑46 (MQ=255)
                                                     aTTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTc                <  1:3310243/68‑1 (MQ=255)
                                                     aTTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTc                <  1:2370121/68‑1 (MQ=255)
                                                     aTTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTc                <  1:3047169/68‑1 (MQ=255)
                                                     aTTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgat             <  1:3255104/71‑1 (MQ=255)
                                                       tGCAGCCCCTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgatga           >  1:2836976/1‑71 (MQ=255)
                                                        gCAGCCCCTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgatgat          >  1:3595236/1‑71 (MQ=255)
                                                        gCAGCCCCTCATCTTCGAAGGGCTGGCTTTTCAGCTTTTCACATTACGTTATAAGAAGTTCCGTCgatga           >  1:3135594/1‑70 (MQ=255)
                                                         cAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgatgat          <  1:3052343/70‑1 (MQ=255)
                                                         cAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgatgat          <  1:1367316/70‑1 (MQ=255)
                                                             cccTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgat             <  1:3852411/63‑1 (MQ=255)
                                                             cccTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgat             <  1:2527212/63‑1 (MQ=255)
                                                             cccTCATGTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgat             <  1:1441325/63‑1 (MQ=255)
                                                                tCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgatgat          <  1:3962385/63‑1 (MQ=255)
                                                                tCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCgatgat          <  1:3229720/63‑1 (MQ=255)
                                                                tCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCGATTATGGCGCTaa  >  1:458343/1‑71 (MQ=255)
                                                                 cATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCGATGATgg        >  1:157287/1‑64 (MQ=255)
                                                                 cATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCGATGATgg        >  1:741996/1‑64 (MQ=255)
                                                                 cATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACATTACGTTATAAGAAGTTCCGTCGATGATgg        >  1:1907927/1‑64 (MQ=255)
                                                                  aTCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaagaag                       >  1:2348530/1‑48 (MQ=255)
                                                                  aTCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaagaag                       >  1:3739735/1‑48 (MQ=255)
                                                                   tCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaagaag                       >  1:2337159/1‑47 (MQ=255)
                                                                   tCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaagaag                       >  1:194821/1‑47 (MQ=255)
                                                                    |                                                                  
CTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAAGAAGTTCCGTCGATGATGGCGCTAA  >  W3110S.gb/3404761‑3404895

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