Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 306,593 T→C 26.2% N655S (AAC→AGC)  yagX ← predicted aromatic compound dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb306,5930TC26.2% 66.0 / 26.1 42N655S (AAC→AGC) yagXpredicted aromatic compound dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (19/12);  new base C (7/4);  total (26/16)
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

CGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATT  >  W3110S.gb/306528‑306663
                                                                 |                                                                      
cgtcCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTa                                                                   >  1:3843670/1‑71 (MQ=255)
              gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:2895877/1‑55 (MQ=255)
              gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTAccc                                                                <  1:1312427/60‑1 (MQ=255)
                gtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTg                                                              >  1:2221455/1‑60 (MQ=255)
                gtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGc                                                   >  1:2592565/1‑71 (MQ=255)
                 ttGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCt                                                  >  1:1102436/1‑71 (MQ=255)
                 ttGAATATCACACCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCt                                                  >  1:1320154/1‑71 (MQ=255)
                   gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:2838998/1‑50 (MQ=255)
                   gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:435342/1‑50 (MQ=255)
                   gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGa                                                      >  1:4256144/1‑65 (MQ=255)
                   gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGc                                                >  1:852700/1‑71 (MQ=255)
                   gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGc                                                >  1:1587267/1‑71 (MQ=255)
                    aaTATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGc                                                >  1:3297202/1‑70 (MQ=255)
                    aaTATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGc                                                <  1:2163701/70‑1 (MQ=255)
                     atatCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGc                                                <  1:155707/69‑1 (MQ=255)
                     atatCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCa                                              >  1:2259845/1‑71 (MQ=255)
                     atatCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCa                                              >  1:1318739/1‑71 (MQ=255)
                      tatCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTAcc                                                                 >  1:3895851/1‑51 (MQ=255)
                       atCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGcc                                                        <  1:839421/59‑1 (MQ=255)
                        tCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCAtt                                            >  1:1952830/1‑70 (MQ=255)
                        tCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTg                                           >  1:2507876/1‑71 (MQ=255)
                          aCCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGc                                         >  1:4270568/1‑71 (MQ=255)
                          aCCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGc                                         >  1:3398662/1‑71 (MQ=255)
                           ccccAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     <  1:2993906/42‑1 (MQ=255)
                                gCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtt                                                                       <  1:1197444/35‑1 (MQ=255)
                                gCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtt                                                                       <  1:3444350/35‑1 (MQ=255)
                                 cGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAg                                  >  1:3767010/1‑71 (MQ=255)
                                      ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGt                             <  1:2984227/71‑1 (MQ=255)
                                      ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGt                             <  1:3087651/71‑1 (MQ=255)
                                         tCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTAtt                          >  1:1140110/1‑71 (MQ=255)
                                         tCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGCCGGTCAAGTTGGTAtt                          <  1:2509039/71‑1 (MQ=255)
                                             tCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              <  1:4528609/63‑1 (MQ=255)
                                              ccgcccgcTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACg                      >  1:2667903/1‑70 (MQ=255)
                                               cgcccgcTGGCAGCAATGTTTTTACCCTGCCAGCCGATGCTGCCATTGGCGGTCAAGTTGGTATTGACGTa                    <  1:2125504/71‑1 (MQ=255)
                                               cgcccgcTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGc                                                <  1:3411065/43‑1 (MQ=255)
                                               cgcccgcTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCAtt                                            >  1:3905872/1‑47 (MQ=255)
                                                gcccgcTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAg                   <  1:3692781/71‑1 (MQ=255)
                                                  ccgcTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTg                               >  1:2710775/1‑57 (MQ=255)
                                                   cgcTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTAtt                          <  1:926423/61‑1 (MQ=255)
                                                      tGGCAGCAATCTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTAtt                          <  1:452690/58‑1 (MQ=255)
                                                       ggCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTa                    >  1:304410/1‑63 (MQ=255)
                                                         cagcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGtt                                >  1:2978977/1‑49 (MQ=255)
                                                         cagcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGtt                                >  1:2792217/1‑49 (MQ=255)
                                                            caATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        >  1:1309295/1‑54 (MQ=255)
                                                             aaTGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTc              <  1:3143230/63‑1 (MQ=255)
                                                              aTGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTa                            >  1:439407/1‑48 (MQ=255)
                                                              aTGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCt     <  1:1332658/71‑1 (MQ=255)
                                                               tGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              <  1:1759576/45‑1 (MQ=255)
                                                                 tttttACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTa    >  1:310117/1‑69 (MQ=255)
                                                                 tttttACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTAtt  <  1:4119076/71‑1 (MQ=255)
                                                                 |                                                                      
CGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATT  >  W3110S.gb/306528‑306663

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: