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

GGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTA‑CCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCA  >  ACB122/5883320‑5883584
                                                                                                                                        |                                                                                                                                  
gggggTCATCCGGGCCCAATGTACCTT‑CGGCGGGGACCCCCCTAAGAAAGCGCGGGCTTGGGGGCCGCTTCGCACAAAAGGGCTGCTCCGGGGGTGGGGGGGCGGGCCTCAAAAAAACCGCCGACACAA‑GGTGTCg                                                                                                                                   <  4:56168/136‑1 (MQ=255)
           cGGCGCAAAGTACCTTACGGCTGGTA‑CCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGcac                                                                                                                         >  5:199141/1‑135 (MQ=255)
           cGGCGCAAAGTACCTTACGGCTGGTA‑CCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGGTGGcac                                                                                                                         >  1:248453/1‑135 (MQ=255)
               gcAAAGTACCTTACGGCTGGTA‑CCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACatat                                                                                                                     >  8:281543/1‑135 (MQ=255)
                                  gTA‑CCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑TGTGTCGGCGTTGGGACATATCAGGCGTTGGTTGCGaaa                                                                                                   >  6:139252/1‑134 (MQ=255)
                                   gtgcccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAAt                                                                                                  >  8:197043/4‑135 (MQ=255)
                                      cccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTCCGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGGGTGGGGGTTGGCCCATATGCGGCGGTGCGAGCGCAATccg                                                                                               >  5:317063/1‑133 (MQ=255)
                                      cccccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCa                                                                                                >  2:248336/1‑134 (MQ=255)
                                                    ggCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑CGTGTGGGCGTTGGCCCATATAAGGCCTTGCGAGCGCAATCAGAACGCCGGCAGTACg                                                                                >  8:31148/1‑136 (MQ=255)
                                                      ctgctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCCCCACAA‑GGTGTCGGCGGTGGCCCCTATCAGGCGGTGCTAGCGAAATCCGCACGCCGGCAGGACGGc                                                                              >  8:70714/1‑136 (MQ=255)
                                                            tGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGCCACAA‑GGTGTCGGCGTTGGCCCATATCAGGCGTTGCTAGCGAAATAAGAACGCCGGCAGTACGGCAccgt                                                                         >  1:357845/1‑133 (MQ=255)
                                                                 ccGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCA‑GGTGTCGGCGGTGGCACCTATCAGGCGGTGCTAGCGAAATCAGAACGCCGGCAGTACGGCCCCGTgggtct                                                                   >  7:234440/1‑132 (MQ=255)
                                                                      tCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGGTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCCCCTTGGTACTTCTTg                                                              >  4:128582/1‑136 (MQ=255)
                                                                                 ggCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGCAATCAGAACGCCGGCAGTACGGCCCCTTGGTACTTCTTGGCaatgaatg                                                    >  2:205349/1‑135 (MQ=255)
                                                                                           gggcggtgCGGGCAGGGCTGAAGAAAAACCCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATTAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                         <  7:226874/131‑1 (MQ=255)
                                                                                           gggcggtgCGGCCTGGGGTCAACAAACACCCCGA‑ACAATGGTGTGGGCGTTGGCACATATCAGGCGTTGCTAGCTAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                         <  7:370875/131‑1 (MQ=255)
                                                                                           gggCTGGGCGGGCAGGGCTGAACAAAACCGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATTAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTcc                                         <  7:232658/134‑1 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCt                                                                          >  7:115467/1‑97 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCt                                                                          <  8:115467/97‑1 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCt                                    >  1:246844/1‑135 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTg                                   >  7:297789/1‑136 (MQ=255)
                                                                                                 tgCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑CGTGTCCGCGTTGGCACATATCACGGGTTGCTAGCGAAATCACAACGCCGGCAGTACGGCACCATGGTACTTCTTGGCAATGAATGCCATGACTCACGGGCt                                    >  6:68535/1‑135 (MQ=255)
                                                                                                    gggCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTTGCAATGAATGCCTTGACTTCCGGGCTgg                                  >  5:171564/1‑134 (MQ=255)
                                                                                                           gCTGAAGAAAAACGCCGACACAA‑GGGGTCGGCGGTGGCACATATCACGCGGTTCTAGCGAAATAACCAACCAGGCACTACGGCACCGTGGTAATTCGTGGCAATGAAAGCCTAGACTTCCGGGCTGGTCAGGCCCt                         >  4:377670/1‑136 (MQ=255)
                                                                                                            tgcataaaaaaCGCCCACACAT‑GGTGGCGGCGTTGGCACATATCAGGCGTTGCTAGCTAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                        <  7:31148/130‑1 (MQ=255)
                                                                                                             ttaagaaAAACGCCGACACAG‑GGTGGCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                        <  6:171564/133‑1 (MQ=255)
                                                                                                             ttaacaaaaaCGCCGACACAA‑GGTGGCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCtt                        <  6:199141/130‑1 (MQ=255)
                                                                                                              gaagaaAAACGCCCACACAA‑TGTGTCGGCGTTTGCACATATCAGGCGTTTCTAGCGAAATCAGAACGCCGGCAGTACGGCCACTTGGTACTTCTTGGCAATGAATGCCTTGACGTCCGGGCTGGTCAGGGCCTTg                       >  2:95974/1‑135 (MQ=255)
                                                                                                              cacaaaaacCCCCCACACAT‑GGTGTGGGGGTTGGAACAAATCAGGCGTTGATAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTg                       <  1:248336/131‑1 (MQ=255)
                                                                                                                       cgccgcCACAA‑GGTGTGGGCGGTGGGACATATCAGGCGGTGGGAGCGAAATCAGAACGCCGGCAGTACGGCCCCGTGGGTCTTCTTGGCAATGTATGGCGTGGCGTCCGGGCTGGTCCGGGGCTTGGGCAGgttg              >  6:10229/1‑134 (MQ=255)
                                                                                                                       cccccacacac‑GGTGGGGGGGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCt             <  7:197043/131‑1 (MQ=255)
                                                                                                                           gACACCG‑GGTGCGGGCGTTGGCACATATCAGGCGTTGCTAGCTAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTGTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGa          <  5:68535/135‑1 (MQ=255)
                                                                                                                              acaG‑GGGGTCGGCGTTGGCTCATATCAGGGGTTGATAGCGAAATCAGATCGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc      <  7:45546/136‑1 (MQ=255)
                                                                                                                                aa‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCt    <  8:297789/136‑1 (MQ=255)
                                                                                                                                   ggTGCCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCa  <  7:281543/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTA‑CCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAA‑GGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCA  >  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: ‑