Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,298,030 G→C 10.4% E83Q (GAG→CAG) ‡ pvdQ → acyl‑homoserine lactone acylase PvdQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,298,0300GC10.4% 139.3 / 7.2 48E83Q (GAG→CAG) ‡pvdQacyl‑homoserine lactone acylase PvdQ
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (2/3);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

GACGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAG  >  NC_002947/3297883‑3298167
                                                                                                                                                   |                                                                                                                                         
gACGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCgtgagcggtga                                                                                                                                          <  3:231659/149‑1 (MQ=255)
  cGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCgtgagcggtgagc                                                                                                                                        <  4:312122/149‑1 (MQ=255)
          cccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCACg                                                                                                                                 <  2:207772/148‑1 (MQ=255)
           ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGc                                                                                                                                >  3:232364/1‑148 (MQ=255)
           ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCACggt                                                                                                                               <  8:199835/149‑3 (MQ=255)
               aTCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTAttt                                                                                                                           <  7:26718/149‑1 (MQ=255)
                  cGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCgg                                                                                                                        >  4:77948/1‑149 (MQ=255)
                         cccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAg                                                                                                                 >  6:43009/1‑149 (MQ=255)
                          ccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAgg                                                                                                                <  5:43009/149‑1 (MQ=255)
                          ccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAgg                                                                                                                >  5:155910/1‑149 (MQ=255)
                             ttACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGcc                                                                                                             <  6:155910/149‑1 (MQ=255)
                                 ggCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAc                                                                                                         >  6:113782/1‑149 (MQ=255)
                                 ggCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAc                                                                                                         >  5:324700/1‑149 (MQ=255)
                                   cGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAccc                                                                                                       <  7:209999/149‑1 (MQ=255)
                                    gTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCt                                                                                                      <  5:8013/149‑1 (MQ=255)
                                       ccGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTg                                                                                                     >  4:95515/1‑147 (MQ=255)
                                       ccGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                   <  5:113782/149‑1 (MQ=255)
                                         gCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAg                                                                                                                 <  7:172947/133‑1 (MQ=255)
                                           atatCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACATATGCCTGCTCGCGAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACACTGGagcag                                                                                               >  3:125954/1‑149 (MQ=255)
                                                gTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGa                                                                                          <  3:77948/149‑1 (MQ=255)
                                                    ccAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAAc                                                                                      >  1:225096/1‑149 (MQ=255)
                                                      aaGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCt                                                                                    >  8:323319/1‑149 (MQ=255)
                                                           tGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCa                                                                               <  4:232364/149‑1 (MQ=255)
                                                                 tGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGAcc                                                                         <  7:323319/149‑1 (MQ=255)
                                                                   gCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTg                                                                       <  6:324700/149‑1 (MQ=255)
                                                                      tGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttc                                                                    >  1:3092/1‑149 (MQ=255)
                                                                        ggCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttc                                                                    <  5:179856/147‑1 (MQ=255)
                                                                                    ggCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCaa                                                      <  2:111414/149‑1 (MQ=255)
                                                                                         tGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGag                                                                                                  >  3:130036/1‑100 (MQ=255)
                                                                                         tGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGag                                                                                                  <  4:130036/100‑1 (MQ=255)
                                                                                             tACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgc                                             >  3:287035/1‑149 (MQ=255)
                                                                                             tACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgc                                             >  7:292753/1‑149 (MQ=255)
                                                                                                      gACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGa                                    >  4:272562/1‑149 (MQ=255)
                                                                                                       aCAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGAt                                   <  2:225096/149‑1 (MQ=255)
                                                                                                        cAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTc                                      >  7:196263/1‑145 (MQ=255)
                                                                                                        cAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTc                                      <  8:196263/145‑1 (MQ=255)
                                                                                                           ccTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGAGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCt                               >  4:117715/1‑149 (MQ=255)
                                                                                                            cTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTg                                                                       <  8:322772/108‑1 (MQ=255)
                                                                                                            cTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTg                                                                       >  7:322772/1‑108 (MQ=255)
                                                                                                              aTGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTcc                            >  8:189098/1‑149 (MQ=255)
                                                                                                                 ccTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCGAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGc                         <  7:274507/149‑1 (MQ=255)
                                                                                                                         ccAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGTCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATTCCTTCCTGCAGGCGCaa                 <  4:125954/149‑1 (MQ=255)
                                                                                                                          cAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGAt                                   >  2:224920/1‑130 (MQ=255)
                                                                                                                          cAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGAt                                   <  1:224920/130‑1 (MQ=255)
                                                                                                                                  gaaatagcGTGAGCGCTCACCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCg        >  8:137733/8‑149 (MQ=255)
                                                                                                                                    cTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCgg       <  3:117715/148‑1 (MQ=255)
                                                                                                                                      gACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTg     >  7:95482/1‑148 (MQ=255)
                                                                                                                                        ccGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAg  <  4:287035/149‑1 (MQ=255)
                                                                                                                                         cGTGAGCGCTCACCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTc                             >  2:270084/1‑121 (MQ=255)
                                                                                                                                         cGTGAGCGCTCACCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTc                             <  1:270084/121‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GACGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAG  >  NC_002947/3297883‑3298167

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