Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,665,749 T→G 5.1% L115R (CTG→CGG)  hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,665,7490TG 5.1% 1326.3 / 2.3 414L115R (CTG→CGG) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base T (83/301);  new base G (25/2);  total (111/303)
Fisher's exact test for biased strand distribution p-value = 7.04e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 426 total aligned reads displayed.

TAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTG  >  CP009273/2665683‑2665816
                                                                  |                                                                   
tAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCt                                                                     <  1:3515438/67‑1 (MQ=255)
 aaCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGc                                                                   <  1:22459949/68‑1 (MQ=255)
  tccATCGGCGAATGCTGTCACTGGGTTCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGcc                                                                  <  1:32176790/67‑1 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGCCCCCTGCCg                                                                 >  1:40518016/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGGCCGCTGCCg                                                                 >  1:33441963/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCTGGCg                                                                 >  1:3539029/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCTGCCg                                                                 >  1:39778220/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCGGCCg                                                                 >  1:19524872/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGc                                                                   >  1:17802724/1‑66 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCg                                                                 >  1:13398478/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCg                                                                 >  1:30226518/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCg                                                                 >  1:29218016/1‑68 (MQ=255)
   ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCGGCCg                                                                 >  1:26677115/1‑68 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCTCTGCCGTTGc                                                             <  1:33036601/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:23060507/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:24061136/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:29877980/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:1647122/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:686409/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:13653805/68‑1 (MQ=255)
       cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGc                                                             <  1:4114323/68‑1 (MQ=255)
        ggCGAAAGCTGGCACTGGGAGCAGCTTTTTATAGCAACCGGCGCGGCACCTAGACCCCAGCCGATgga                                                            >  1:20687439/1‑66 (MQ=255)
          cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCtt                                                           <  1:2843185/67‑1 (MQ=255)
          cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCtt                                                           <  1:15362990/67‑1 (MQ=255)
             aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATg                                                       <  1:15625560/68‑1 (MQ=255)
             aaGCTGGCACTGGGAGCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATg                                                       <  1:35973966/68‑1 (MQ=255)
              aGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGc                                                      <  1:18010857/68‑1 (MQ=255)
              aGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCACCTCGACCCCGGCCCTTGGGTGATGc                                                      >  1:33138387/1‑68 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCa                                                     <  1:17344136/66‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCa                                                     <  1:14214178/66‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCa                                                     <  1:38435214/66‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACt                                                   <  1:3654060/68‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACt                                                   <  1:36861827/68‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACt                                                   <  1:18157601/68‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACt                                                   <  1:10781187/68‑1 (MQ=255)
                 tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACt                                                   <  1:26876632/68‑1 (MQ=255)
                 tGGCACTGGGATCAGCGTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACt                                                   <  1:9608519/68‑1 (MQ=255)
                 tGGCACGGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGCTTGATGCa                                                     <  1:12878342/66‑1 (MQ=255)
                  ggCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTg                                                  <  1:19022451/68‑1 (MQ=255)
                     aCTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTggg                                                <  1:22707107/67‑1 (MQ=255)
                     aCTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGga                                               >  1:3239602/1‑68 (MQ=255)
                     aCTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGATGCACTGGga                                               >  1:14386697/1‑68 (MQ=255)
                     aCTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTGCTTGATGCACTGGga                                               >  1:8457575/1‑68 (MQ=255)
                     aCTGGGATCAGCTTCTTATAGCAACCGGCGCGGCCGCTCGCCCGCGTCCCTTGCTTGATTCACAgggc                                               >  1:34027160/1‑67 (MQ=255)
                         ggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAAc                                           >  1:931577/1‑68 (MQ=255)
                          gATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCTGCCGTTGCTTGATGCCCTGGGAGAAc                                           >  1:1657675/1‑67 (MQ=255)
                          gATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAAc                                           >  1:31044321/1‑67 (MQ=255)
                            tCAGCTTTTTTTAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACgct                                        <  1:29551034/68‑1 (MQ=255)
                            tCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACgct                                        <  1:33637071/68‑1 (MQ=255)
                            tCAGCTTTTTATACCAACCTGCGGGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGTGAACgct                                        <  1:36301006/68‑1 (MQ=255)
                             cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACgctg                                       <  1:37586263/68‑1 (MQ=255)
                               gCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCTGCCGTTGCTTGATGCACTGGGAGCACgcggcg                                     >  1:1576711/1‑67 (MQ=255)
                               gCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCTGCCGGTGCTTGATTCCCTGGGGGAACgctgct                                     >  1:38174370/1‑68 (MQ=255)
                                 tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGGCGGTGCTTGCTGCACTGGGACAACGCTGCttt                                   >  1:24581394/1‑68 (MQ=255)
                                 tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACgctgc                                      <  1:10910026/65‑1 (MQ=255)
                                 tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACgctgc                                      <  1:40854126/65‑1 (MQ=255)
                                 tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACgctgc                                      <  1:6510159/65‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCAGTGGGAGAACGCTGCTTTAc                                 <  1:41223812/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:41845291/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:4000170/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:41657200/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:5129676/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:14972149/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:13823310/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:42979980/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:4912493/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:20113487/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:20249976/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:3467536/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:2540463/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:32847201/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:12195539/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:27946713/68‑1 (MQ=255)
                                   tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAc                                 <  1:3137554/68‑1 (MQ=255)
                                    ttATAGCAACCGGCGCGGCAGCTCGACCGCTTCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAcc                                <  1:18686566/68‑1 (MQ=255)
                                    ttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAcc                                <  1:43623309/68‑1 (MQ=255)
                                    ttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTAcc                                <  1:889849/68‑1 (MQ=255)
                                        aGCAACCGGCGCGGCCGCGCGCCCGCTGCCGTTGCTTGATGCACGGGGTGAACGCTGGTTTACCCTa                             >  1:5169289/1‑67 (MQ=255)
                                        aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTAc                            <  1:43211042/68‑1 (MQ=255)
                                        aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTAc                            <  1:34400531/68‑1 (MQ=255)
                                         gCAACCGGCGCGGCAGCTCGGCCGCTGCCGTTGCTTGATGCACGGGGAGAACGCTGCGTTACCCTAc                            >  1:20350720/1‑67 (MQ=255)
                                         gCAACCGGCGCGGCAGCTCGCCCGCTGCCGTTGGTTGATGCACTGGGGGAACGCGGCTTTTCCCTAc                            >  1:37148613/1‑67 (MQ=255)
                                           aaCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACg                           <  1:40387684/66‑1 (MQ=255)
                                           aaCCGGCGCGGCAGCTCGACCGCTGCCGGTGGGTGATGCACTGGGAGAAAGCTGGATTACCCTACGcc                         >  1:25706595/1‑68 (MQ=255)
                                           aaCCGGCGCGGCAGCGCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGc                          <  1:42246902/67‑1 (MQ=255)
                                             ccGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCAt                       <  1:28850174/68‑1 (MQ=255)
                                               ggCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGc                     <  1:21582918/68‑1 (MQ=255)
                                               ggCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGc                     <  1:43537847/68‑1 (MQ=255)
                                                      cAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGa               <  1:39467346/67‑1 (MQ=255)
                                                      cAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGAt              >  1:39181948/1‑68 (MQ=255)
                                                         cgcGACCGCTGCCGTTGCTTGTTGCACTGGGACAACGCTGCTTTACCCTACGCCAAGCCGGTGATgcc           <  1:21259766/66‑1 (MQ=255)
                                                         cTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgcc           >  1:10288472/1‑68 (MQ=255)
                                                          tCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgcc           <  1:38711292/67‑1 (MQ=255)
                                                          tCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgcc           <  1:6164913/67‑1 (MQ=255)
                                                          tCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgcc           <  1:25049809/67‑1 (MQ=255)
                                                          tCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgcc           <  1:11813013/67‑1 (MQ=255)
                                                          tCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATgcc           <  1:10029873/67‑1 (MQ=255)
                                                               cGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGa     <  1:7796739/68‑1 (MQ=255)
                                                                  tGCCGTTGCTTGTTGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTg  <  1:27600750/68‑1 (MQ=255)
                                                                  tGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTg  <  1:30640489/68‑1 (MQ=255)
                                                                  |                                                                   
TAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTG  >  CP009273/2665683‑2665816

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