Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,984,502 T→G 23.2% V270G (GTG→GGG)  U712_09805 → Putative 4‑hydroxyphenylacetate 3‑monooxygenase yoaI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,984,5020TG23.2% 65.3 / 4.1 37V270G (GTG→GGG) U712_09805Putative 4‑hydroxyphenylacetate 3‑monooxygenase yoaI
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base G (7/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAGGCGTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGA  >  CP006881/1984367‑1984640
                                                                                                                                       |                                                                                                                                           
cAAAAGGCTTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGTATTCAATAGTGGTTTTTGATCATGTGTAAGTGCCGTGGAATCgtgtgtg                                                                                                                                          >  2:107882/1‑138 (MQ=255)
  aaaGGCGTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGtttt                                                                                                                                        >  2:132584/1‑139 (MQ=255)
   aaGGCGTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGttttt                                                                                                                                       >  8:366367/1‑139 (MQ=255)
     ggCGTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGttttttt                                                                                                                                     >  7:367119/1‑139 (MQ=255)
      gCGTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGttttttt                                                                                                                                     >  7:336728/1‑138 (MQ=255)
        gTAAAATCAATCACAACCTAATCATTTGTATAAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACg                                                                                                                                  <  8:186963/139‑1 (MQ=255)
                                        aaGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGAATCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGGGTTTTTTTAAGATAATGTTGAAGCAGCGAAAGAtttttttac                                                                                                   >  1:139893/1‑133 (MQ=255)
                                         aGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGGTTTTTTAGGGTAATTGTGAAGCAGCTAAAGATTTTTTGACa                                                                                                  >  7:5450/1‑138 (MQ=255)
                                         aGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGGGGTTTTTTTAGAATATGTGTGAGCAGCGAAAGAgttttttata                                                                                                  >  5:233898/1‑132 (MQ=255)
                                           cGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACaaa                                                                                                >  8:288504/1‑138 (MQ=255)
                                             agacagccTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAg                                                                                             <  3:151083/133‑1 (MQ=255)
                                              ctctTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGGGGTTTTTTTTGATATAGTGGTAGCACCGTAAGAGTTTATGGCAAAAAGc                                                                                            >  6:370461/1‑139 (MQ=255)
                                              ctctTCCTTTAATCATCCCTTAAGCTCAAGGAATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGGGTTTTTTTTAGATTATGTTTAAGAAACTTAATATTTTTTTTCaaaaata                                                                                            >  7:14995/1‑137 (MQ=255)
                                                 ttCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGct                                                                                           <  7:11140/137‑1 (MQ=255)
                                                      ttAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCa                                                                                    <  3:53653/139‑1 (MQ=255)
                                                       tAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCa                                                                                    <  3:191945/138‑1 (MQ=255)
                                                       tAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCa                                                                                    <  5:49337/138‑1 (MQ=255)
                                                               ccTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTAc                                                                          <  2:169657/139‑1 (MQ=255)
                                                                      cTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCAc                                                                   <  8:336728/139‑1 (MQ=255)
                                                                      cTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTAAGATAATGTGGAAGCAGCATAATATTTTATTACAAAAAAGAGTTTTTAAGCCTTTACAtttgaa                                                                   >  1:78392/1‑136 (MQ=255)
                                                                        cAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGATTTTTTAAGTTAATTTGGATGGCCGGAAACATTTTATTTCAAAAATCTCATTTTATGCCTTTACCTTTAAccc                                                                 >  4:90310/1‑138 (MQ=255)
                                                                          aGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAgg                                                               <  3:213138/139‑1 (MQ=255)
                                                                                     aaTGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGGGTTTTTTTACGATAATGTGGCTGCAGCGAAAGATTTTATTTCAAAAAGCTCGTTTTATTCCGTTTCCTTTCACCACGGTGTCGTTTgg                                                    >  5:6547/1‑139 (MQ=255)
                                                                                      aTGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGGGTTTTTTTTAGTTAATGTGGTAGGCGCGAAAAAATTTTTTAAAAAAATCACTTTTTATTGCATTTACAGTTACCAGGATGTGATTAgtg                                                   >  5:230322/1‑137 (MQ=255)
                                                                                              aaTTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGAt                                            <  1:132584/138‑1 (MQ=255)
                                                                                                        tttGATCATGTGTTAGTGCCGTGGAATCGTGGGTTTTTTTAAGATAAAGTTGGTGCAGCAAAAAATTTTTTTTAAAAAAACGCGTTTAATGGCGTTAAAGTTGCACACGTTTTAATTAGGGAAAAGATaaaaaattaaa                                 >  4:277242/1‑133 (MQ=255)
                                                                                                        acagATCATGTGTTAGTGCCGTGTAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAAt                                 <  1:183904/136‑1 (MQ=255)
                                                                                                         ttGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAAt                                 <  6:233898/138‑1 (MQ=255)
                                                                                                         ttGATCATGGGGTATTGCCGGGGAATCGGGTGTTTTTTTACGATTATGTTGAATCATCTAAAGATTTTATTACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAAt                                 <  2:78392/138‑1 (MQ=255)
                                                                                                              cATGTGGTAGTGCCGTGGACTCGAGGTTCTTTTTACGATAAAGTTGAAGCAGCTAAAGATGTTATGACACAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGCTAAAACTTGAATTTCTa                            <  4:218851/138‑1 (MQ=255)
                                                                                                               aTGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTAtt                          <  6:230322/139‑1 (MQ=255)
                                                                                                               aTGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTAtt                          <  7:288504/139‑1 (MQ=255)
                                                                                                                  tgtTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAg                        <  6:245560/138‑1 (MQ=255)
                                                                                                                          gcGTGGAATAGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGca                <  1:107882/137‑1 (MQ=255)
                                                                                                                          ccGTGGAATCGTGGGGTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTATCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGca                <  8:80420/138‑1 (MQ=255)
                                                                                                                          cagTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTCAATTTCTATTAGGTGTTGcac               <  8:367119/137‑1 (MQ=255)
                                                                                                                            gTGGAATCGTGTGTTTTTTTTAGATAATTTGGAAGCAGCGAAAGATTTTTTTACAAAAAGCTCTTTTGATTGCTTTTCCTTTGACCAGGTTGTGATTAGGCAAATTATAAAAATTTAATTTATATTTGGTTTTTGACAg             >  8:99482/1‑139 (MQ=255)
                                                                                                                             tGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGc            <  1:192081/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAATCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCTCCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTTCACAGCTTCTTGTCGa  <  1:212761/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCATCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  3:267074/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  7:99482/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  1:20522/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  5:370461/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  7:205485/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  3:277242/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  7:138265/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  3:277831/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  5:204956/139‑1 (MQ=255)
                                                                                                                                       tgtTTTTTTACGATAATGGGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGa  <  7:366367/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CAAAAGGCGTAAAATTCATCACAAGAGAATCATTTGTATTAAGCGACTCTTCCTTTAATCATCCCTTAAGCTCAAGGTATGAAGAAATGGATTCAATTGTGGTTTTTGATCATGTGTTAGTGCCGTGGAATCGTGTGTTTTTTTACGATAATGTGGAAGCAGCTAAAGATTTTATGACAAAAAGCTCTTTTCATGCCTTTACCTTTCACCAGGTTGTCATTAGGCAAATGATAAAAATTGAATTTCTATTAGGTGTTGCACAGCTTCTTGTCGA  >  CP006881/1984367‑1984640

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