Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 306,593 T→C 44.7% 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,5930TC44.7% 6.6 / 49.2 38N655S (AAC→AGC) yagXpredicted aromatic compound dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (13/8);  new base C (4/13);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 2.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATT  >  W3110S.gb/306523‑306663
                                                                      |                                                                      
aCCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGt                                                                        >  1:1508628/1‑71 (MQ=255)
      gtcCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTAc                                                                  <  1:2846106/71‑1 (MQ=255)
      gtcCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTAc                                                                  <  1:2978666/71‑1 (MQ=255)
      gtcCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTAc                                                                  <  1:1316627/71‑1 (MQ=255)
      gtcCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTAc                                                                  <  1:3661824/71‑1 (MQ=255)
           ccAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTa                                                                   >  1:3566625/1‑65 (MQ=255)
                   gtgtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGAc                                                     <  1:4069697/71‑1 (MQ=255)
                     gtTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGc                                                   >  1:1892798/1‑71 (MQ=255)
                      ttGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     <  1:1036363/52‑1 (MQ=255)
                      ttGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtt                                                                       <  1:3836604/50‑1 (MQ=255)
                        gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGc                                                >  1:2426114/1‑71 (MQ=255)
                        gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGc                                                >  1:3973315/1‑71 (MQ=255)
                        gAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGc                                                >  1:104612/1‑71 (MQ=255)
                         aaTATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGcc                                               <  1:1436318/71‑1 (MQ=255)
                           tatCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     <  1:2479302/47‑1 (MQ=255)
                                  ccAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGt                                      >  1:687961/1‑71 (MQ=255)
                                     gCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCa                                              <  1:3234384/60‑1 (MQ=255)
                                     gCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCa                                              <  1:2939235/60‑1 (MQ=255)
                                      cGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:1822193/1‑36 (MQ=255)
                                      cGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:544627/1‑36 (MQ=255)
                                      cGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:2079081/1‑36 (MQ=255)
                                      cGTTGCCATCAGTCCGCCCGCTGGCAGCAATGtttt                                                                     >  1:240409/1‑36 (MQ=255)
                                        ttGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGtt                                >  1:3689678/1‑71 (MQ=255)
                                          gCCATCAGTCCGCCCGCTGGCAGCAATGTTCTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              >  1:2657923/1‑71 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              <  1:3874174/70‑1 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGt                             <  1:2529328/71‑1 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGt                             <  1:1168841/71‑1 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGt                             <  1:3365970/71‑1 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              <  1:935587/70‑1 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              <  1:912041/70‑1 (MQ=255)
                                           ccATCAGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTgg                              <  1:481100/70‑1 (MQ=255)
                                              tCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTAtt                          >  1:2077784/1‑71 (MQ=255)
                                               cAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTAtt                          >  1:3142095/1‑70 (MQ=255)
                                                aGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        >  1:1977493/1‑71 (MQ=255)
                                                aGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        >  1:1460433/1‑71 (MQ=255)
                                                aGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        >  1:1484242/1‑71 (MQ=255)
                                                aGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        >  1:1575462/1‑71 (MQ=255)
                                                aGTCCGCCCGCTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        <  1:26479/71‑1 (MQ=255)
                                                    cgcccgcTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACg                      >  1:1835795/1‑69 (MQ=255)
                                                    cgcccgcTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACg                      >  1:3621446/1‑69 (MQ=255)
                                                          cTGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTc              <  1:3684632/71‑1 (MQ=255)
                                                           tGGCAGCAATGCTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACg                      <  1:1475376/62‑1 (MQ=255)
                                                             gcagcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTccgc           >  1:3337509/1‑71 (MQ=255)
                                                                gcaATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGa                        <  1:3606134/55‑1 (MQ=255)
                                                                  aaTGTTTTTACCCTGCCAGCCGACGCTGACATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCgcgc      <  1:2301690/71‑1 (MQ=255)
                                                                    tGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTc                                     <  1:201068/38‑1 (MQ=255)
                                                                      tttttACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTAtt                          <  1:3032673/47‑1 (MQ=255)
                                                                      tttttACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTAtt  >  1:401811/1‑71 (MQ=255)
                                                                      tttttACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTAtt  >  1:1147200/1‑71 (MQ=255)
                                                                      |                                                                      
ACCGTCGTCCTCCAGCCCGGTGTTGAATATCACCCCAGCGTTGCCATCAGTCCGCCCGCTGGCAGCAATGTTTTTACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATT  >  W3110S.gb/306523‑306663

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: