Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,881,004 G→A 9.6% M1I (ATG→ATA) † quiP → Acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,881,0040GA 9.6% 121.1 / 4.1 42M1I (ATG→ATA) †quiPAcyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base G (22/15);  new base A (2/2);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TAGAGCTGGCCGGTGAGCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCTGCAAGCGCCGGTCAGCGT  >  ACB122/5880872‑5881138
                                                                                                                                    |                                                                                                                                      
tAGAGCTTGCCTGTGAGCATTTTGTGTTTTTCGCTGAGGTACATGGAAACTACTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCtgaatgaa                                                                                                                                      <  6:54244/135‑1 (MQ=255)
tAGAGCTGGCCGGTGAGCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCtgaatgaa                                                                                                                                      <  6:54248/135‑1 (MQ=255)
      tGGCCGGTGAGCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCCCATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAg                                                                                                                                >  7:290549/1‑135 (MQ=255)
         ccGGTGAGCATTTTGTGTTGTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGcc                                                                                                                              <  5:267615/134‑1 (MQ=255)
           ggTGAGCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGTATGAAGCGCAGCCTGa                                                                                                                           >  4:118823/1‑135 (MQ=255)
           ggTGAGCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGa                                                                                                                           >  5:140442/1‑135 (MQ=255)
           ggTGAGCATTTAGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTATTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGa                                                                                                                           >  7:219918/1‑135 (MQ=255)
                gCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCct                                                                                                                       <  1:273479/134‑1 (MQ=255)
                      tggggttttCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGc                                                                                                                 <  8:286883/129‑1 (MQ=255)
                              tCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATTAAGCGCAGCCTGACACTGCTGGCGGTGGTAg                                                                                                         >  4:236264/1‑134 (MQ=255)
                                          aTGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGc                                                                                            <  1:293509/135‑1 (MQ=255)
                                            ggAAACTCCTTCCTGGGGCTATGCTCTGGTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTg                                                                                          <  2:252447/135‑1 (MQ=255)
                                                    cttcctGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGc                                                                                            >  7:163307/1‑125 (MQ=255)
                                                    cttcctGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGc                                                                                            <  8:163307/125‑1 (MQ=255)
                                                       cctGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCgg                                                                                >  5:86601/1‑134 (MQ=255)
                                                           ggggTATGGTCTGTTCGAGATTCCGCGGATGATCGGACCTGATTTTTCGTTTGAATGCACTAGGAGTGTGTATGAAGCGCAGCCTGACCCTGCTTGCGGTGGTAGTCAGCGTGGCGGCTTGCGCCTCCGGGTACt                                                                           <  8:219918/135‑1 (MQ=255)
                                                           gggCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACt                                                                           <  2:25251/135‑1 (MQ=255)
                                                            ggCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATTCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTg                                                                          <  5:34481/135‑1 (MQ=255)
                                                            ggCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACt                                                                           <  1:161456/134‑1 (MQ=255)
                                                             gCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGGAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTgg                                                                         >  3:244856/1‑135 (MQ=255)
                                                               tATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTa                                                                       >  7:111500/1‑135 (MQ=255)
                                                                aTGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTAc                                                                      >  2:251916/1‑135 (MQ=255)
                                                                   ctctGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGt                                                                    >  4:203088/1‑134 (MQ=255)
                                                                    tctGTTCGCGAATCCGCCGATGATCGGACCACATATTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTCTCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGCTACGt                                                                    >  6:192361/1‑133 (MQ=255)
                                                                       gTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCa                                                                 <  5:95772/133‑1 (MQ=255)
                                                                       gTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCa                                                                 <  8:303573/133‑1 (MQ=255)
                                                                         tCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGATGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAggg                                                              >  1:200424/1‑134 (MQ=255)
                                                                         tCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCggtg                                                                                                             <  7:34264/87‑1 (MQ=255)
                                                                         tCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCggtg                                                                                                             >  8:34264/1‑87 (MQ=255)
                                                                         tCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAggg                                                              >  6:13463/1‑134 (MQ=255)
                                                                         tCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATCCTCAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAgg                                                               <  8:6986/133‑1 (MQ=255)
                                                                          cgcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGGACGTTCAggg                                                              >  5:165349/1‑133 (MQ=255)
                                                                          cgcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGCTACTGGGACGTTCAggg                                                              >  5:168860/1‑133 (MQ=255)
                                                                           gcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCCGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCaa                                                           >  5:221742/1‑135 (MQ=255)
                                                                           gcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGCCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCaa                                                           >  2:239249/1‑135 (MQ=255)
                                                                           gcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATAAAGCGCAGCCTGACCctgc                                                                                                                     >  8:95651/1‑77 (MQ=255)
                                                                           gcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATAAAGCGCAGCCTGACCctgc                                                                                                                     >  8:100120/1‑77 (MQ=255)
                                                                           gcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATAAAGCGCAGCCTGACCctgc                                                                                                                     <  7:95651/77‑1 (MQ=255)
                                                                           gcgAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATAAAGCGCAGCCTGACCctgc                                                                                                                     <  7:100120/77‑1 (MQ=255)
                                                                              aaTCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCt                                                        >  6:181532/1‑135 (MQ=255)
                                                                                     cgATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAg                                                 >  7:232779/1‑135 (MQ=255)
                                                                                              gACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCgcgagggc                                        >  1:112437/1‑135 (MQ=255)
                                                                                                                 aTGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGGCGCCGTTGCCGGcc                     <  1:268436/135‑1 (MQ=255)
                                                                                                                   gCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCt                    >  7:141942/1‑134 (MQ=255)
                                                                                                                               tgactgaaGCGCAGCCGGACCCTGCTGGCGGTGGTAGGCACCGTGGCCGCTGGCGCCGCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCTGCAAGCGCCGGTc       <  7:260978/135‑1 (MQ=255)
                                                                                                                                  atgaaGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCTGCAAGCGCCGGTCAGc    >  7:223595/1‑135 (MQ=255)
                                                                                                                                    gaaGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCTGCAAGCGCCGGTCAGCGt  >  7:126226/1‑135 (MQ=255)
                                                                                                                                    gaaGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCTGCAAGCGCCGGTCAGCGt  >  3:52530/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TAGAGCTGGCCGGTGAGCATTTTGTGTTTTTCGCTGAGGGACATGGAAACTCCTTCCTGGGGCTATGCTCTGTTCGCGAATCCGCCGATGATCGGACCACATTTTTCGTTTGAATGCACAAGGAGTCTGAATGAAGCGCAGCCTGACCCTGCTGGCGGTGGTAGTCACCGTGGCCGCTGGCGCCTCCGGGTACTGGTACGTTCAGGGCAAGCTGCCGCAGCGCGAGGGCGAGGTCGCCGTTGCCGGCCTGCAAGCGCCGGTCAGCGT  >  ACB122/5880872‑5881138

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