Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,653,495 G→A 11.8% intergenic (+61/+37) PP_3218 → / ← PP_3219 NtaA/SnaA/SoxA family monooxygenase/alkansulfonate monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,653,4950GA11.8% 105.6 / 4.3 42intergenic (+61/+37)PP_3218/PP_3219NtaA/SnaA/SoxA family monooxygenase/alkansulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (3/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AGGCACGTGGGGTGTATCGCACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGAAGCCACG  >  NC_002947/3653349‑3653642
                                                                                                                                                  |                                                                                                                                                   
aGGCACGTGGGGTGTATCGCACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTg                                                                                                                                                     >  5:61251/1‑147 (MQ=255)
   catgttgggTGTATCTCTCGGCTTATGAGGGCACGACCCTTCTTGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCTCCTGGGCCGGCCTCTTCGCGTCTAACAACCGTTACTACAGGAATCTTCTATCACCTGTAGGa                                                                                                                                                <  6:328004/143‑1 (MQ=255)
           gtgtATCGCACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCtt                                                                                                                                        >  2:65912/1‑149 (MQ=255)
               aTCGCACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTtcgtc                                                                                                                                    >  8:35930/1‑149 (MQ=255)
                  gCACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTcgc                                                                                                                                  >  6:328520/1‑148 (MQ=255)
                   cACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTATAGGAACGTCCTTCGTcgcg                                                                                                                                 >  3:147852/1‑148 (MQ=255)
                    aCGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTCCAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGaa                                                                                                                               <  1:65912/149‑1 (MQ=255)
                      ggCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAaga                                                                                                                             <  2:107527/149‑1 (MQ=255)
                              aGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTATAGGAACGTCCTTCGTCGCGAag                                                                                                                              <  6:720/140‑1 (MQ=255)
                              aGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTATAGGAACGTCCTTCGTCGCGAag                                                                                                                              >  5:720/1‑140 (MQ=255)
                                ggCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCa                                                                                                                   <  6:67814/149‑1 (MQ=255)
                                     gACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCa                                                                                                              <  7:217778/149‑1 (MQ=255)
                                       cccTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTATAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAgg                                                                                                            <  4:147852/149‑1 (MQ=255)
                                        ccTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGgc                                                                                                           <  1:295808/149‑1 (MQ=255)
                                          ttCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGgcgc                                                                                                         >  7:256816/1‑149 (MQ=255)
                                             gggAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGc                                                                                                      <  8:85260/149‑1 (MQ=255)
                                                 aaGCTTGGGGGTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGt                                                                                                  <  2:58508/149‑1 (MQ=255)
                                                             ggCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTGGCGGCTAACAACCGTTCCGACAGGAATCTTCAATCACCTGTAGGAACGGCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGc                                                                                      >  7:276087/1‑149 (MQ=255)
                                                              gCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCt                                                                                     <  3:280956/149‑1 (MQ=255)
                                                                gATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTc                                                                                   >  7:268338/1‑149 (MQ=255)
                                                                 aTACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTcc                                                                                  <  6:61251/149‑1 (MQ=255)
                                                                 aTACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTcc                                                                                  <  2:5269/149‑1 (MQ=255)
                                                                    ccGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTgg                                                                               <  7:35930/149‑1 (MQ=255)
                                                                     cgcgGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGc                                                                                                      >  8:77610/1‑125 (MQ=255)
                                                                     cgcgGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGc                                                                                                      <  7:77610/125‑1 (MQ=255)
                                                                          tCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCt                                                                         <  2:237535/149‑1 (MQ=255)
                                                                              gCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCa                                                                     <  8:247456/149‑1 (MQ=255)
                                                                                          ccTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGt                                                                                                                      <  6:85232/88‑1 (MQ=255)
                                                                                              ggccggccTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCgg                                                     <  8:256816/149‑1 (MQ=255)
                                                                                                     ctctTGGCGGCTAACAACCGTTCCTACACGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAgc                                              >  1:191901/1‑149 (MQ=255)
                                                                                                     ctctTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAgc                                              >  4:231888/1‑149 (MQ=255)
                                                                                                               cTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCa                                     >  2:60623/1‑148 (MQ=255)
                                                                                                                  acaacCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATTAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAgcgc                                                                 <  6:66171/117‑1 (MQ=255)
                                                                                                                    aacCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTc                                >  2:40606/1‑148 (MQ=255)
                                                                                                                      ccGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCAc                             <  3:231888/149‑1 (MQ=255)
                                                                                                                             tACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACg                            <  7:129744/143‑1 (MQ=255)
                                                                                                                             tACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCAGGTGCAGATCAGCCGCTGCGCGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCACGAATCCGCGGAATCAGAGCCTGACCATACTCCACg                            >  8:129744/1‑143 (MQ=255)
                                                                                                                                     tCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACt               >  4:235983/1‑148 (MQ=255)
                                                                                                                                       ttCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGAt            >  2:37237/1‑149 (MQ=255)
                                                                                                                                        tCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATc           >  6:336440/1‑149 (MQ=255)
                                                                                                                                            tCACCTATAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGAAg       >  6:191302/1‑149 (MQ=255)
                                                                                                                                              aCCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGAAGc      >  6:306123/1‑148 (MQ=255)
                                                                                                                                                 tGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGAAGCCACg  >  2:298940/1‑149 (MQ=255)
                                                                                                                                                  gTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGaa        >  2:255233/1‑142 (MQ=255)
                                                                                                                                                  gTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGaa        <  1:255233/142‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGGCACGTGGGGTGTATCGCACGGCGTATGAGGGCACGACCCTTCGGGAAAGCTTGGGGTTGGCGATACCGCGGTCTGGCCGCTAGATCGCCTGGGCCGGCCTCTTCGCGGCTAACAACCGTTCCTACAGGAATCTTCAATCACCTGTAGGAACGTCCTTCGTCGCGAAGAGGCCGGTGCAGATCAGGCGCTGCGGGTGCTGCGACCGGCTTCCTGGTCTGCTACCAGCGCCCGAATCCGCGGAATCAGCGCCTGCCCATACTCCACGGCATCACGCAACTGATCGAAGCCACG  >  NC_002947/3653349‑3653642

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