Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,348,914 C→G 21.9% G142G (GGC→GGG)  U712_17210 → Para‑nitrobenzyl esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,348,9140CG21.9% 78.8 / 6.1 43G142G (GGC→GGG) U712_17210Para‑nitrobenzyl esterase
Reads supporting (aligned to +/- strand):  ref base C (10/23);  new base G (6/3);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 6.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTCATGGTGTGGATTCACGGAGGCGCTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTAAC  >  CP006881/3348780‑3349039
                                                                                                                                      |                                                                                                                                
ccTGTCATGGTGTGGATTCACGGAGGCGCTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAATTCATTGTCGTTTCTTTGAACTATCGGCTGGGGCCGTTTTGCgttt                                                                                                                              >  7:12931/1‑139 (MQ=255)
          tgtgGATTCACGGAGGCGCTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑c                                                                                                                   <  8:362497/139‑1 (MQ=255)
          tgtgGATTCACGGAGGCGCTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑c                                                                                                                   <  3:354745/139‑1 (MQ=255)
                    cGGAGGCGCTTTTTATCTCGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAAt                                                                                                         <  8:185198/139‑1 (MQ=255)
                    cGGAGGCGCTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAAt                                                                                                         <  1:237783/139‑1 (MQ=255)
                            cTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTa                                                                                                 <  1:155186/139‑1 (MQ=255)
                               tttATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCTGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTAtt                                                                                               <  7:50061/138‑1 (MQ=255)
                                   tctcGGAGCGGGGAGTGAGGCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGa                                                                                           <  1:342092/138‑1 (MQ=255)
                                   tCTTGGAGCGGTCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGTGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGa                                                                                           <  5:213510/138‑1 (MQ=255)
                                   tCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGa                                                                                           <  5:21602/138‑1 (MQ=255)
                                   tCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGa                                                                                           <  5:133125/138‑1 (MQ=255)
                                       ggAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTGTGGCTTTTTTGCCTTTT‑GGTCAATTTATGAGGCGGTTTATTAT‑AAc                                                                                      >  4:221471/1‑138 (MQ=255)
                                          gCGGGCAGTGAGCCATTTTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTg                                                                                 <  4:237821/139‑1 (MQ=255)
                                          gCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTg                                                                                 <  7:235121/139‑1 (MQ=255)
                                                aGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTGGGCttt                                                                           <  8:165060/139‑1 (MQ=255)
                                                              ttaCGTATAAAAAGTTGGGGTACAGTGTGATGTCTTTGTCGTTACATTTAACTATGGGCTGGGGCCGTTTGGGTTTTTGCACTTGT‑CTTCATTAATTGTGGCGTATTGTGAT‑AACC‑TTGGGCTTTTAGACCAAgccgc                                                              <  4:153481/136‑1 (MQ=255)
                                                              tGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCt                                                                                                                                 >  4:227029/1‑74 (MQ=255)
                                                              tGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCt                                                                                                                                 <  3:227029/74‑1 (MQ=255)
                                                              tGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTTGACGTGT‑CTTAATTTTATGAGGCGTTTTATGAT‑AAAC‑TTGGGCTTTTAGACCAAgccgcc                                                             >  6:273018/1‑139 (MQ=255)
                                                                aCGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTGGGCTTTTAGACCAAGCcgccgc                                                           <  5:155239/139‑1 (MQ=255)
                                                                aCGGATAAAAACTTGCGGCACAGGTAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCGGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATTAGGCGTATTCTGAT‑AACC‑TTGGGCTTTTAGACCAAGCcgccgc                                                           <  5:61098/139‑1 (MQ=255)
                                                                    aTCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGGTTTTTGCACTTTT‑GTTCATTTAATGAGGCGGATTCTGAT‑AAAC‑ATGGGGTTTTAGACCAAGCCGCCccgcg                                                        >  1:54961/1‑137 (MQ=255)
                                                                    aTCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTTGACTTGT‑CTTACATTAATTAGGCGTATTGAGAT‑AACA‑CTGGGGTTTTAGACCCACCCCCCGCGCt                                                        >  3:153481/1‑138 (MQ=255)
                                                                    aTCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑GTTCATTTAATGAGGCGTTTTCTGAT‑AACC‑TTGGGGTTTTAGACCAAGCCGCCGCGCt                                                        >  7:246849/1‑138 (MQ=255)
                                                                    aTCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTGTGTGTTTTTTGACTTGTGTTTCATTTATTAGGGGGTTTTTTATT‑AAAA‑TT‑GGCTTTTATAACAAGCCGCCccgcgg                                                       >  5:296073/1‑137 (MQ=255)
                                                                    aTCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTGTGGGGTTTTTGACGTTT‑CTTAATTTAATTAGGCGGATTATGAT‑AACC‑TTGGGGTTTTTGACCAACCCCCCccgcgg                                                       >  3:303910/1‑137 (MQ=255)
                                                                              gCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTTT‑CTTAAATTAATGAGGCGTAATCAGAT‑AACC‑ATGGGCTTTTATACCAAGCCGCCGCGCTTAAATGGGAGc                                             >  6:129069/1‑139 (MQ=255)
                                                                              gCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTGTGGGTTTTTTGACGTGT‑GTTGATTTTATTTGGCGGTTTGTGGT‑TACCGTTGGGGTTTT‑TACCAAGCCGCCGCGCGGCAATtggggg                                             >  3:227202/1‑136 (MQ=255)
                                                                                    ccgggATATGTCATTGTGGTAAATTTGAAGTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑ATTAATTTAATGAGGCTTATTCTGAT‑AACT‑TTGGGATTTTCGACAAAGCCGCCGCGCTGAAATGGGTGCgagaga                                       <  8:156902/137‑1 (MQ=255)
                                                                                                              gAACTATCGGCTGGGGCCGTTTGGCGTTTTTGACTTGT‑ATTAATTTAATTAGGGGGATTTTTATAAACA‑TTGGGATTTT‑TACCAAACCCCCCCCCGGAAAAGTGTTCGAGAGAAAATTTTAAGGTTTTGGGGTgtttcc             >  4:314188/1‑134 (MQ=255)
                                                                                                               aaCTATCGGCTGGGGCCGTGTGGGGTTTTTGGCTTGT‑GTTGATTTTATTAGGCGGATTATTAT‑AAAC‑ATTGGCGTTTAGACCCAACCCCCGCCCGGAAATGGGTGGGATAGAAAAATTAAGCGGTTTGCGGGGATaccc            >  4:61747/1‑138 (MQ=255)
                                                                                                                aCTATCGGCTGGGGCGGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGCGGCGTATTCTGAT‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGa           <  7:45242/139‑1 (MQ=255)
                                                                                                                    tCGGCTGGGGCCGTTTGGCTTTTTGGGCGTTT‑GTTGATTTAATAGGGCGTGATATTAT‑AACC‑AGTTGCGTTTTGACCAAGCCCCCGCCCGGAAATAGGTTGCAGAGAAAATTTCAGCGTTTGGCGGTGATAGCGATa         >  8:314094/1‑137 (MQ=255)
                                                                                                                    tCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGCTCCCGATAAc       >  1:23466/1‑139 (MQ=255)
                                                                                                                      ggggggggCCTTTTGGCGTTTTTCACTTGC‑TTTCTTTAAATGCGGCTTATTCTAAT‑ACCC‑TTGGGCTTTTAGCCCACGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGt     <  6:12314/135‑1 (MQ=255)
                                                                                                                      ggCGGGGGCCGTTTGGTTTTTGGCACTGGT‑CTTCATTTAATGCGGCTTATTCTGAA‑ACCC‑TTGGGCTTTTCCACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGt     <  4:265212/139‑1 (MQ=255)
                                                                                                                       gCTGGGGCCGTTTGGCTTTTTGCAGTTGT‑CTTCATTTAATGAGGCTTTTTCTTAT‑AACC‑TTGGGCTTTTAGACAAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTa    <  2:54961/139‑1 (MQ=255)
                                                                                                                         tGGGGCCGTTTTGGTTTTTTCACTTGT‑GTTGATTTAATTAGGGGTATTGTGAT‑AACC‑TTGGGGTTTTATACCAAGCCCCCGCGCGGAAATTGGTGCGAGAGAAAATTTAAGCGTTTGGCGGTGATCCCGATAACCTaa   >  6:99516/1‑138 (MQ=255)
                                                                                                                         tGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCTTTTAATGAGGCGTATTCTTTT‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTAAc  <  5:129069/139‑1 (MQ=255)
                                                                                                                         tGGGGCCGTTTGGCTTTTGGCACTTGT‑ATACTTTTACTGCGCCGTATTCTGAT‑CACG‑TTGGGCTTTTCGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTAAc  <  3:221471/139‑1 (MQ=255)
                                                                                                                         tGGGGCCGGTTGGGTTTTTGCAATTGG‑TTTAATTTAATGAGGCGTATTCTTAA‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGGTTGGCGGTGATCCCGATAACGTAAc  <  3:61747/139‑1 (MQ=255)
                                                                                                                         gggggCGGTTTGGGTTTTTGCAGTTGT‑CTTCATTTAATGAGGCGTATGCTTAT‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTAAc  <  7:296737/138‑1 (MQ=255)
                                                                                                                         gggggCCGTTGGGCTTTTGGCCCTTTT‑CTTAATTTAATGAGGCGTATTCTGAT‑AACC‑TTGGGATTTTACACCACGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTAAc  <  1:164086/138‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
CCTGTCATGGTGTGGATTCACGGAGGCGCTTTTTATCTAGGAGCGGGCAGTGAGCCATTGTATGACGGATCAAAACTTGCGGCACAGGGAGAAGTCATTGTCGTTACATTGAACTATCGGCTGGGGCCGTTTGGCTTTTTGCACTTGT‑CTTCATTTAATGAGGCGTATTCTGAT‑AACC‑TTGGGCTTTTAGACCAAGCCGCCGCGCTGAAATGGGTGCGAGAGAATATTTCAGCGTTTGGCGGTGATCCCGATAACGTAAC  >  CP006881/3348780‑3349039

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: ‑