Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,883,454 C→G 19.0% intergenic (+86/+34) quiP → / ← metQ_2 Acyl‑homoserine lactone acylase QuiP/putative D‑methionine‑binding lipoprotein MetQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,883,4540CG19.0% 63.6 / 2.9 35intergenic (+86/+34)quiP/metQ_2Acyl‑homoserine lactone acylase QuiP/putative D‑methionine‑binding lipoprotein MetQ
Reads supporting (aligned to +/- strand):  ref base C (18/10);  new base G (3/4);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCA  >  ACB122/5883320‑5883584
                                                                                                                                        |                                                                                                                                  
gggggTCATCCGGGCCCAATGTACCTTCGGCGGGGACCCCCCTAAGAAAGCGCGGGCTTGGGGGCCGCTTCGCACAAAAGGGCTGCTCCGGGGGTGGGGGGGCGGGCCTCAAAAAAACCGCCGACACAAGGTGTCg                                                                                                                                   <  4:56168/136‑1 (MQ=255)
           cGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGcac                                                                                                                         >  5:199141/1‑135 (MQ=255)
           cGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGGTGGcac                                                                                                                         >  1:248453/1‑135 (MQ=255)
               gcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACatat                                                                                                                     >  8:281543/1‑135 (MQ=255)
                                  gTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAATGTGTCGGCGTTGGGACATATCAGGCGTTGGTTGCGaaa                                                                                                   >  6:139252/1‑134 (MQ=255)
                                   gtgcccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAt                                                                                                  >  8:197043/4‑135 (MQ=255)
                                      cccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTCCGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGGGTGGGGGTTGGCCCATATGCGGCGGTGCGAGCGCAATccg                                                                                               >  5:317063/1‑133 (MQ=255)
                                      cccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCa                                                                                                >  2:248336/1‑134 (MQ=255)
                                                    ggCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAACGTGTGGGCGTTGGCCCATATAAGGCCTTGCGAGCGCAATCAGAACGCCGGCAGTACg                                                                                >  8:31148/1‑136 (MQ=255)
                                                      ctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCCCCACAAGGTGTCGGCGGTGGCCCCTATCAGGCGGTGCTAGCGAAATCCGCACGCCGGCAGGACGGc                                                                              >  8:70714/1‑136 (MQ=255)
                                                            tGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGCCACAAGGTGTCGGCGTTGGCCCATATCAGGCGTTGCTAGCGAAATAAGAACGCCGGCAGTACGGCAccgt                                                                         >  1:357845/1‑133 (MQ=255)
                                                                 ccGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCAGGTGTCGGCGGTGGCACCTATCAGGCGGTGCTAGCGAAATCAGAACGCCGGCAGTACGGCCCCGTgggtct                                                                   >  7:234440/1‑132 (MQ=255)
                                                                      tCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGGTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCCCCTTGGTACTTCTTg                                                              >  4:128582/1‑136 (MQ=255)
                                                                                 ggCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGCAATCAGAACGCCGGCAGTACGGCCCCTTGGTACTTCTTGGCaatgaatg                                                    >  2:205349/1‑135 (MQ=255)
                                                                                           gggcggtgCGGGCAGGGCTGAAGAAAAACCCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATTAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                         <  7:226874/131‑1 (MQ=255)
                                                                                           gggcggtgCGGCCTGGGGTCAACAAACACCCCGAACAATGGTGTGGGCGTTGGCACATATCAGGCGTTGCTAGCTAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                         <  7:370875/131‑1 (MQ=255)
                                                                                           gggCTGGGCGGGCAGGGCTGAACAAAACCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATTAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                         <  7:232658/134‑1 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCt                                                                          >  7:115467/1‑97 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCt                                                                          <  8:115467/97‑1 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCt                                    >  1:246844/1‑135 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTg                                   >  7:297789/1‑136 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAACGTGTCCGCGTTGGCACATATCACGGGTTGCTAGCGAAATCACAACGCCGGCAGTACGGCACCATGGTACTTCTTGGCAATGAATGCCATGACTCACGGGCt                                    >  6:68535/1‑135 (MQ=255)
                                                                                                    gggCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTTGCAATGAATGCCTTGACTTCCGGGCTgg                                  >  5:171564/1‑134 (MQ=255)
                                                                                                           gCTGAAGAAAAACGCCGACACAAGGGGTCGGCGGTGGCACATATCACGCGGTTCTAGCGAAATAACCAACCAGGCACTACGGCACCGTGGTAATTCGTGGCAATGAAAGCCTAGACTTCCGGGCTGGTCAGGCCCt                         >  4:377670/1‑136 (MQ=255)
                                                                                                            tgcataaaaaaCGCCCACACATGGTGGCGGCGTTGGCACATATCAGGCGTTGCTAGCTAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                        <  7:31148/130‑1 (MQ=255)
                                                                                                             ttaagaaAAACGCCGACACAGGGTGGCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                        <  6:171564/133‑1 (MQ=255)
                                                                                                             ttaacaaaaaCGCCGACACAAGGTGGCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                        <  6:199141/130‑1 (MQ=255)
                                                                                                              gaagaaAAACGCCCACACAATGTGTCGGCGTTTGCACATATCAGGCGTTTCTAGCGAAATCAGAACGCCGGCAGTACGGCCACTTGGTACTTCTTGGCAATGAATGCCTTGACGTCCGGGCTGGTCAGGGCCTTg                       >  2:95974/1‑135 (MQ=255)
                                                                                                              cacaaaaacCCCCCACACATGGTGTGGGGGTTGGAACAAATCAGGCGTTGATAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTg                       <  1:248336/131‑1 (MQ=255)
                                                                                                                       cgccgcCACAAGGTGTGGGCGGTGGGACATATCAGGCGGTGGGAGCGAAATCAGAACGCCGGCAGTACGGCCCCGTGGGTCTTCTTGGCAATGTATGGCGTGGCGTCCGGGCTGGTCCGGGGCTTGGGCAGgttg              >  6:10229/1‑134 (MQ=255)
                                                                                                                       cccccacacacGGTGGGGGGGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCt             <  7:197043/131‑1 (MQ=255)
                                                                                                                           gACACCGGGTGCGGGCGTTGGCACATATCAGGCGTTGCTAGCTAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTGTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGa          <  5:68535/135‑1 (MQ=255)
                                                                                                                              acaGGGGGTCGGCGTTGGCTCATATCAGGGGTTGATAGCGAAATCAGATCGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc      <  7:45546/136‑1 (MQ=255)
                                                                                                                                aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCt    <  8:297789/136‑1 (MQ=255)
                                                                                                                                   ggTGCCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCa  <  7:281543/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCA  >  ACB122/5883320‑5883584

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: