Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 502,322 T→G 28.0% intergenic (‑262/‑106) estP ← / → trpG Esterase EstP/Anthranilate synthase component 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440502,3220TG28.0% 66.4 / 4.2 58intergenic (‑262/‑106)estP/trpGEsterase EstP/Anthranilate synthase component 2
Reads supporting (aligned to +/- strand):  ref base T (14/26);  new base G (16/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGGGACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAACACCATGTTACTGATGATCGACAATTACGACT  >  KT2440/502195‑502455
                                                                                                                               |                                                                                                                                     
tgCCGGGACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCCCCCGGccccgccccggt                                                                                                                                >  4:171437/1‑133 (MQ=255)
 gCCGGGACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGttt                                                                                                                               <  2:241854/135‑1 (MQ=255)
   cGGGACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCCGTTGc                                                                                                                             >  2:385472/1‑135 (MQ=255)
     ggACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCGGTTGccg                                                                                                                           >  3:275171/1‑134 (MQ=255)
     ggACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCCGCCCGGCCCCGTTTGCCt                                                                                                                           >  3:124465/1‑135 (MQ=255)
      gACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCGGTTGCCGtc                                                                                                                         >  8:348364/1‑136 (MQ=255)
          tATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGGTTGGCGTCTCg                                                                                                                      >  4:253703/1‑135 (MQ=255)
             gAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGttt                                                                                                                               >  5:234434/1‑123 (MQ=255)
             gAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGttt                                                                                                                               <  6:234434/123‑1 (MQ=255)
                       tATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCCCCCGGCCCGGCCCCGGTGTGCGTCGCGGTTACCGCCCGtt                                                                                                         >  5:145496/1‑135 (MQ=255)
                       tATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCCCCCCGCCCCGCCCCGTTTGGCGTCGAGTTTTGCGCCCCTtg                                                                                                        >  7:437704/1‑136 (MQ=255)
                            tgtgTTGATCGTCAGCAAACTCTTGATGCGCGTCATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCCCCCGGGCCCTACCCGTTTTCGTTGTCGTTTAACGCCCGTAGTgg                                                                                                     >  1:55462/1‑134 (MQ=37)
                            tgtgTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCCCCCCGTGCCCCTTTGCAGTCGGGTTTAAAGCCCGATGTggg                                                                                                    >  3:441629/1‑134 (MQ=255)
                                 tGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTGCCCGGTTGCCGTCTCGGTTACCGCCCGTTGTGGCAAAGc                                                                                               >  2:10178/1‑135 (MQ=255)
                                 tGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCCGCCCCCGTTGGCGTGTCGGTTAACCCCCGGTTGGGCAAAGc                                                                                               >  5:333617/1‑135 (MQ=255)
                                 tGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCCGCCCGGCCCCCGTTGCCGTGGAGGTTTACGCCCGGTGTGGCGAAGc                                                                                               >  1:332460/1‑135 (MQ=255)
                                 tGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCCCCCCGCCCCTCCCCGTTTGGCATGACGTTTTACGCCCGGTGTGGGAAAGc                                                                                               >  6:341884/1‑135 (MQ=255)
                                   aTCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCGTTTTCCGTATCGTTTTCAGCCCGTTGGGGCAAAGccc                                                                                             >  3:334782/1‑134 (MQ=255)
                                      gTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTGCCCGGTTGGCTGCTCGTTTAACGCCCGTTGTGGCAAAGCCGCGcc                                                                                         >  8:233209/1‑136 (MQ=255)
                                      gTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGGTTGCCTTCTCGTTTACCGCCCCTTGTGGGAAAGCCgcgc                                                                                          >  5:180342/1‑135 (MQ=255)
                                      gTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGGTTGCCTTCTCGTTTAACGCCCCTTGTGGCAAAGCCCCGcc                                                                                         >  8:14437/1‑136 (MQ=255)
                                      gTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCCGCCCGGCCCCGGTTGGCTTCTCGTTTTCCGCCCGGTGTGGGCAAGCCGCGcc                                                                                         >  7:4657/1‑136 (MQ=255)
                                      gTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCCCCCGGCCCCGCCCCGGTTGGGGGGTCGGTTTCCGCCCGGTGTGGCGAAGcccccg                                                                                          >  5:31918/1‑131 (MQ=255)
                                           cAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCCGCCCCTCCCCCATTTCCTTCTCGGTTAACGCCCCGTGTGGGAAAGCCCCGCCGGTgg                                                                                    >  8:322203/1‑136 (MQ=255)
                                                        cgcgTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGaa                                                                        <  7:14437/135‑1 (MQ=255)
                                                         gcgTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGaa                                                                        <  2:259048/134‑1 (MQ=255)
                                                             aaTCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCCGCCCCGTTTGGCGTCTGGTTTACCGCCCGTTGTGGCAAATCAGCGCCGGTTGGTGGTGCCGGAAGGtt                                                                    >  1:50433/1‑134 (MQ=255)
                                                             aaTCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCCGCCCGGCCCCGTTTGGCTTCTCGTTTTCCGCCCGGTGTGGCAAAGCCGCGCCGGGGGGGGGGGCCGGGAgggtg                                                                   >  3:280860/1‑132 (MQ=255)
                                                               tCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGGCCGCCCGGCCCGGCCCCCGTTGGCGTCGCGGTTTCCGCCCGGTGTGGGAAACCCGCCCCGGTGGCGGGTGGCGGGGGGTTCggg                                                                >  7:389044/1‑134 (MQ=255)
                                                                 aaGGCCCCCGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAAt                                                               <  8:424202/135‑1 (MQ=255)
                                                                 aaGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAAt                                                               <  8:430811/135‑1 (MQ=255)
                                                                   ggCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCa                                                             <  4:229560/135‑1 (MQ=255)
                                                                      ccACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCAGACCCGCCCCGCCCCGTTTGGAATGTCGATTACCCCCAGGTGTTGAAAAGCCGCGCCTCTGGGTGGGGCCGGTAGGTTCGAATCAttga                                                         >  8:45660/1‑133 (MQ=255)
                                                                        aCGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCGCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGTCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGTAg                                                        <  3:149636/135‑1 (MQ=255)
                                                                           ccggttggcgttggcGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCa                                                     <  3:38893/132‑1 (MQ=255)
                                                                           cccGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCGTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGGAGGGTCGAATGATGAAGTCa                                                     >  6:417241/1‑135 (MQ=255)
                                                                            ccGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGGTTGGGTTGGCGTTTTCCGCCCGTTGTGGGCAAGCCGCGCCTGTGGGTGGGGGCGGGAGGTTGGATTCATGAAGTAAg                                                    >  4:147572/1‑135 (MQ=255)
                                                                                  gcgttggcgCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTc                                             <  8:295587/136‑1 (MQ=255)
                                                                                  gcgttggcgCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCCGTTGGCGTGTCGGTTACCGCCCGGTGTGGCCAAGCCGCGCCGGTGGGGGGGGCCGGCAGGGTCGAATCATGTAGTGAGTAGAtt                                              >  8:78874/1‑135 (MQ=255)
                                                                                   cgttggcgCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTc                                             <  5:221284/135‑1 (MQ=255)
                                                                                   cgttggcgCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTc                                             <  6:311623/135‑1 (MQ=255)
                                                                                   cgttggcgCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTc                                             <  4:5490/135‑1 (MQ=255)
                                                                                   cggtggcgCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTGCTCGTTTACAGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTc                                             <  2:398951/132‑1 (MQ=255)
                                                                                        gcgcgcCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAagag                                        <  7:322203/135‑1 (MQ=255)
                                                                                        gcgcgcCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAgg                                       <  7:348364/136‑1 (MQ=255)
                                                                                          gcgcCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGGTTTGCGAAGCTTTTTACTCCCGTGTTTGCAAATCCGCGCGATTGTCTGGTGGCGGGGAGTTGGGATAATGTAATCAGAAGATTGAAGAggg                                      >  2:350730/1‑134 (MQ=255)
                                                                                            gcCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTa                                   <  8:437704/136‑1 (MQ=255)
                                                                                             ccTTACTTTATCTCCTGCCGGCCGGCCGGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTa                                   <  1:189642/135‑1 (MQ=255)
                                                                                              cTGACTTTATCTCCTGCCGGCCCGCCGTGCCCGTCGCCGTTTGCCTTATCGTTTACCGCGCGTTGTGGCAAAGCCGCGCGTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTaa                                  <  3:171437/135‑1 (MQ=255)
                                                                                              cTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTaa                                  <  1:271116/135‑1 (MQ=255)
                                                                                              cTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTaa                                  <  2:332460/135‑1 (MQ=255)
                                                                                               tGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAac                                 <  2:396164/135‑1 (MQ=255)
                                                                                               tGACTTTATCTCCTGCCGGCCCGCCCGGCCCGGCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGCCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGGTTAac                                 <  1:296427/135‑1 (MQ=255)
                                                                                                               cGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAACACCATGTTACtgat                   <  4:389142/133‑1 (MQ=255)
                                                                                                               cGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAACACCATGTTACtgat                   <  4:334782/133‑1 (MQ=255)
                                                                                                                ggCCCGCCCGGCCGGTCCCCGTTTGCCTTCGCGTTTACCGCCTGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCTAATCATGAAGTCAGTAGATTCATGAGTTTTAACACCATGTTACTgatgat                <  4:124465/135‑1 (MQ=255)
                                                                                                                ggCCCGCCCGGCCCGTCCCCGTTTGGCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCGGTGGGTGGTGGCGGGAGGTTCGATTCATGGAGTCAGTAGATTGAAGAGGTTTAACACCAGGTTACtgattctc               >  8:348492/1‑132 (MQ=255)
                                                                                                                    cgcccgGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAACACCATGTTACTGATGATCGAc            <  4:275171/135‑1 (MQ=255)
                                                                                                                              gTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAACACCATGTTACTGATGATCGACAATTACGACt  <  3:253703/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                     
TGCCGGGACGTATGAGCGAGACCTATCCTGTGTTGATCGTCAGCAAACTCTTGATGCGCGTAATCAAGGCCCACGCCCGTTGGCGTTGGCGCGCCTGACTTTATCTCCTGCCGGCCCGCCCGGCCCGTCCCCGTTTGCCTTCTCGTTTACCGCCCGTTGTGGCAAAGCCGCGCCTGTGGCTGGTGCCGGAAGGTTCGAATCATGAAGTCAGTAGATTCAAGAGGTTTAACACCATGTTACTGATGATCGACAATTACGACT  >  KT2440/502195‑502455

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