Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,453,575 T→C H23R (CAC→CGC)  paaZ ← oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,453,5750TC97.6% 136.9 / ‑4.8 43H23R (CAC→CGC) paaZoxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/16);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGTGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGAATCATAAAAAGCAAAG  >  NC_000913/1453445‑1453708
                                                                                                                                  |                                                                                                                                     
cgccgcACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCACTAATACCGTGGCGAATCAAAc                                                                                                                              >  1:268010/1‑140 (MQ=255)
 gccgcACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACg                                                                                                                             >  1:654736/1‑140 (MQ=255)
 gccgcACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACg                                                                                                                             >  1:446731/1‑140 (MQ=255)
     gcgtttCGATAAAGGTCATAGCGCGAAGGGGGGGGGCGCCTTTTACAATGGCAAACGGGGGGGGAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTa                                                                                                                         <  2:133633/136‑1 (MQ=255)
     cACGTTCGTTAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTa                                                                                                                         >  2:526979/1‑140 (MQ=255)
     cACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTa                                                                                                                         >  2:574576/1‑140 (MQ=255)
     cACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCAAAACACCTCGCCGCTAATAGCGCGGCGAATCAAACGGCTa                                                                                                                         >  2:236129/1‑140 (MQ=255)
          tCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGcc                                                                                                                    >  2:442212/1‑140 (MQ=255)
             aTAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCg                                                                                                                 >  1:27409/1‑140 (MQ=255)
                aaGGTCATAGCGCGAGGGGGGGGGGCCCTTTTTTCAATGGCAACCGGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGcc                                                                                                              <  2:31201/140‑1 (MQ=255)
                   gTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGa                                                                                                           >  1:662239/1‑140 (MQ=255)
                        ggcgcGAGGGGGGGGGGCTCTTTTTTAAATGGAAAACTGGGGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGcc                                                                                                      <  2:662239/139‑1 (MQ=255)
                             gAAGGGCGGGGGCACATTTTTAGATCGCACACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTa                                                                                                 <  1:312706/140‑1 (MQ=255)
                              aaGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAGTAGCGTGGGTAAGCAAACGGCTACGGCCCCTGCCAGACTGCCAGGTAc                                                                                                >  2:179875/1‑140 (MQ=255)
                                           ccTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTAAGGCCCCGGCCAAACGGCCAGGTACCGGATAAGAAACt                                                                                   >  1:328201/1‑140 (MQ=255)
                                               tttCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGACCCCGGCGCGACTGCCTGGTACCGGATAAGAAACTGGCt                                                                               >  2:247454/1‑140 (MQ=255)
                                                      ggCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAActgct                                                                        >  2:126621/1‑140 (MQ=255)
                                                           aCTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATc                                                                   >  2:66263/1‑140 (MQ=255)
                                                             tGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGc                                                                 <  1:372923/140‑1 (MQ=255)
                                                                     acccgccATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTcc                                                         <  1:574576/138‑1 (MQ=255)
                                                                     aGCCGCCATATCAAGACCTTCAGTGGTCACTTCCCATCACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTcc                                                         <  1:236129/140‑1 (MQ=255)
                                                                          ccATATCAAGACCTTCACTGGTAACCTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATg                                                    <  1:206343/140‑1 (MQ=255)
                                                                                  aGACCTTCACTGGTCCCTTCCCCTACCGCGTCGCCGCGAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACAt                                            <  1:205841/140‑1 (MQ=255)
                                                                                    aCCTTCACTGGTCACCTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACAttt                                          <  1:247454/140‑1 (MQ=255)
                                                                                       ttCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGTGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGACACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTg                                       <  2:429543/140‑1 (MQ=255)
                                                                                        tCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCGTCGCTACTCTCCAAATGTTTCACATTTCTGt                                      >  2:236420/1‑140 (MQ=255)
                                                                                        tCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGt                                      >  1:704/1‑140 (MQ=255)
                                                                                        tCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGt                                      >  2:115001/1‑140 (MQ=255)
                                                                                         cACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGtt                                     >  1:389285/1‑140 (MQ=255)
                                                                                            tGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCt                                  <  2:446731/140‑1 (MQ=255)
                                                                                              gTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGGTAACTGCTGCATCGCTACTCTCCAGGTGTTGCACATTTCTGTTGCTaa                                >  2:267683/1‑140 (MQ=255)
                                                                                              gTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTaa                                >  1:597301/1‑140 (MQ=255)
                                                                                              gTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTaa                                >  1:422609/1‑140 (MQ=255)
                                                                                              gTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTaa                                >  1:475082/1‑140 (MQ=255)
                                                                                                   ttCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCt                                  >  2:530980/1‑133 (MQ=255)
                                                                                                   ttCCCATAACGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCt                                  <  1:530980/133‑1 (MQ=255)
                                                                                                           aCGCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATcgcg                   <  1:442212/140‑1 (MQ=255)
                                                                                                             gCCTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGaa                 >  1:34475/1‑140 (MQ=255)
                                                                                                               cTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGaa                 >  2:441976/1‑138 (MQ=255)
                                                                                                               cTCGCCGCTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGaa                 <  1:441976/138‑1 (MQ=255)
                                                                                                                   ccgcTAATAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGAATCATaa           <  2:422609/140‑1 (MQ=255)
                                                                                                                         aTAGCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGAATCATAAAAAGCa     <  1:115001/140‑1 (MQ=255)
                                                                                                                            gCGTGGCGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGAATCATAAAAAGCAAAg  >  2:536487/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                     
CGCCGCACGTTCGATAAAGGTCATAGCGCGAAGGGCGGGGGCACCTTTTTCAATGGCAAACTGGCGGGCAGCCGCCATATCAAGACCTTCACTGGTCACTTCCCATAACGCCTCGCCGCTAATAGCGTGGTGAATCAAACGGCTACGGCCCCGGCCAGACTGCCAGGTACCGGATAAGAAACTGGCTAACTGCTGCATCGCTACTCTCCAGATGTTTCACATTTCTGTTGCTAATAGTTAAATCGCGAATCATAAAAAGCAAAG  >  NC_000913/1453445‑1453708

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