Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,227,068 A→C 7.5% intergenic (+26/+19) fadH → / ← higA 2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked/antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,227,0680AC 7.5% 646.4 / 2.5 211intergenic (+26/+19)fadH/higA2,4‑dienoyl‑CoA reductase, NADH and FMN‑linked/antitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (124/66);  new base C (2/17);  total (127/84)
Fisher's exact test for biased strand distribution p-value = 3.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

CGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTT  >  CP009273/3227001‑3227135
                                                                    |                                                                   
cgGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCaa                                                                     <  1:29521844/68‑1 (MQ=255)
  ggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAAc                                                                    >  1:6182873/1‑67 (MQ=255)
  ggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAAc                                                                    >  1:16209215/1‑67 (MQ=255)
  ggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAAc                                                                    >  1:3287784/1‑67 (MQ=255)
    cAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 >  1:9793533/1‑68 (MQ=255)
    cAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 >  1:27543589/1‑68 (MQ=255)
    cAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 >  1:4487344/1‑68 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGACGCGCCGGGGGCAACGCAt                                                                <  1:5756019/68‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 <  1:32288363/67‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 <  1:11813147/67‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 <  1:4664635/67‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 <  1:23712511/67‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 <  1:15191351/67‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCa                                                                 <  1:767947/67‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                <  1:3921078/68‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                <  1:15867502/68‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                <  1:9657772/68‑1 (MQ=255)
     aaTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                <  1:30825346/68‑1 (MQ=255)
      aTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                >  1:32711587/1‑67 (MQ=255)
      aTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                >  1:8962962/1‑67 (MQ=255)
      aTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                >  1:21959486/1‑67 (MQ=255)
      aTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                >  1:20572302/1‑67 (MQ=255)
      aTTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCAt                                                                >  1:11899740/1‑67 (MQ=255)
        tGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCg                                                             >  1:30182576/1‑68 (MQ=255)
        tGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCg                                                             >  1:24896901/1‑68 (MQ=255)
        tGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCg                                                             >  1:19163082/1‑68 (MQ=255)
        tGCCCAGGGAACAAGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGAATCCg                                                             >  1:26169544/1‑68 (MQ=255)
         gCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCg                                                             >  1:22055759/1‑67 (MQ=255)
          cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGc                                                           >  1:11021009/1‑68 (MQ=255)
           ccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGc                                                           >  1:23175546/1‑67 (MQ=255)
           ccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGc                                                           >  1:21060721/1‑67 (MQ=255)
             aGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACg                                                        >  1:696390/1‑68 (MQ=255)
               ggAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:1852063/1‑68 (MQ=255)
               ggAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:18686017/1‑68 (MQ=255)
               ggAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:21492375/1‑68 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGTCGGATGCAACGCATCCGGCACGTc                                                      >  1:23779142/1‑67 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCACCGCATCCGGCACGTc                                                      >  1:5547670/1‑67 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:10382940/1‑67 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:7896930/1‑67 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:17211233/1‑67 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:14773270/1‑67 (MQ=255)
                gAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTc                                                      >  1:24478974/1‑67 (MQ=255)
                 aaCCCGGCGGGCGCTGGATATTTAAATCTTTATTTCACGCGCGGTATGCAACGCATCGGGCACGTcgc                                                    <  1:19678218/68‑1 (MQ=255)
                 aaCACGGGTGGCGCTGGATATTTAAATGGTTATTTGAAGCGCCTGATGAAACGCATCCGGCACGTcgc                                                    <  1:34852383/68‑1 (MQ=255)
                  acacGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTcgcg                                                   >  1:14950556/1‑68 (MQ=255)
                   cacGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGt                                                  >  1:22407365/1‑68 (MQ=255)
                     cGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGCACGTCGCGTTa                                                >  1:13330163/1‑68 (MQ=38)
                     cGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTa                                                >  1:26379777/1‑68 (MQ=255)
                     cGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTa                                                >  1:3616582/1‑68 (MQ=255)
                       gggggCGCTGTATTTTAAAATGTTTTTTTCACGCGCCGTATCCACCCCAGCGGGCACGTCGCGTTaa                                               <  1:7211226/64‑1 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:272926/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:13563082/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:28390557/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:18200077/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:33567598/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:15544696/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:5211240/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                               >  1:19668073/1‑67 (MQ=255)
                       gCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATACGGCACGTCGCGTTaa                                               >  1:10013573/1‑67 (MQ=255)
                         tGGCGGTGTATATTAAAATCTTTTTTTGACGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATc                                             <  1:13701438/67‑1 (MQ=25)
                         tGGCGCTGTATTTTAAAATCGTTATTTGACCCGCCGGATGCAACGAAGCGGGCACGTCGCGTTAATc                                             <  1:3163432/67‑1 (MQ=255)
                         tGGCGCGGGATATTAAAAGCGTTATTTGAAGCGCGGGAGGCACCGCACCCGGCACGTCGCGTTAATc                                             <  1:34576387/67‑1 (MQ=255)
                         gggCGCTGGATATTAAAATCGTTATTTCACGCGCGGGAGCCAACGCACCCGGCACGTCGCGTTAATCa                                            <  1:2017348/67‑1 (MQ=255)
                            cgcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATa                                         >  1:16493875/1‑68 (MQ=255)
                             gcTGTATATTTAAATGTTTTTAAGAACCGGCGCTTGAAAGGCAGCCGGCACGTCGCGTTAATCAATaa                                        <  1:25351349/68‑1 (MQ=255)
                             gcTGTATATTTAAATGGTTATTTGAAGCGCCGTATGCAACGCAGCGGGCACGTCGCGTTAATCAATaa                                        <  1:4847045/68‑1 (MQ=255)
                             gcTGGATATTTAAATCGTTATTTGAAGCGCGGGATGCAACGCATCCGGCACGTCGCGTTAATCAATaa                                        <  1:33287794/68‑1 (MQ=255)
                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATaa                                        >  1:9811006/1‑68 (MQ=255)
                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATaa                                        >  1:28737155/1‑68 (MQ=255)
                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATACGGCCCGTCGCGTTAATCAATaa                                        >  1:32440625/1‑68 (MQ=255)
                               tGGATATTTAAATCTTTATTTGAAGCGCCGGATGCAACGCATCGGGCACGTCGCGTTAATCAATaaa                                       <  1:16647124/67‑1 (MQ=255)
                                  ttaTTTAAATGGTTATTTTACGCCCGGTATGAACGGCTGCGGGCACGTCGCGTTAATCAATAAACAAg                                   <  1:25644649/66‑1 (MQ=255)
                                  ataTTTAAATGTTTTTTTCACGCGCGGTATCCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAg                                   <  1:22770271/66‑1 (MQ=255)
                                  ataTTTAAATCTTTTTTTGAAGCGCGGTATGAAACGAATCGGGCACGTCGCGTTAATCAATAAACAAg                                   <  1:27125652/66‑1 (MQ=255)
                                  agaTTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAg                                   >  1:16918958/1‑68 (MQ=255)
                                   taTTTAAATCTTTATTTGACGCGCGGTTTGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAgg                                  <  1:33507518/67‑1 (MQ=255)
                                    ttttAAATTGTTTTTTTCACGCGCGGTATGCAACGAATCCGGCACGTCGCGTTAATCAATAAACAAgg                                  <  1:24370490/67‑1 (MQ=255)
                                     tttAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCg                                >  1:12253847/1‑68 (MQ=255)
                                        aaaTCTTTTTTTTACGCGCGGTTTGCACCGCAGCGGGCACGTCGCGTTAATCAATAAACAAGGCGGga                             <  1:14827307/68‑1 (MQ=255)
                                        aaaTCTTTATTTGAAGCGCGGGTTGCAACGAATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGga                             <  1:11527411/68‑1 (MQ=255)
                                        aaaTCTTTATTTCACGCGCCGGATCCAACCCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGga                             <  1:26865771/68‑1 (MQ=255)
                                          aTCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgaga                           >  1:24033761/1‑68 (MQ=255)
                                          aTCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgaga                           >  1:18008547/1‑68 (MQ=255)
                                            cGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGaaa                         >  1:14413436/1‑68 (MQ=255)
                                            cGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGaaa                         >  1:10268817/1‑68 (MQ=255)
                                            cGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGaaa                         >  1:29966894/1‑68 (MQ=255)
                                             gTTATTTGAAGCGCGGTATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAAt                        <  1:31596418/68‑1 (MQ=255)
                                                     gAAGCGCCGGATGCAACGCATCAGGCCCGGCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCgca                  >  1:34066975/1‑65 (MQ=255)
                                                       aGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGGGGGAGAAATGCCGAATcg               >  1:28210357/1‑68 (MQ=255)
                                                          gcCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGt            >  1:3512713/1‑68 (MQ=255)
                                                          gcCGGATGCAACGCATCCGGCACGTCGCGGTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGt            >  1:25766118/1‑68 (MQ=255)
                                                              taTGCAACCCAGCGGGCACGCCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGc         <  1:23337334/66‑1 (MQ=255)
                                                               aTGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCa       >  1:24704766/1‑68 (MQ=255)
                                                               aTGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCa       >  1:25612324/1‑68 (MQ=255)
                                                                 gaaaCGCTTCGGGCACGTCGCGTTAATCAATAAACAAGGCTGGAGAAATGCCGAATCGCGTTGCCAAt     <  1:17393743/66‑1 (MQ=255)
                                                                 gCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAt     >  1:17041732/1‑68 (MQ=255)
                                                                 gCAACCCAGCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAt     <  1:20973608/68‑1 (MQ=255)
                                                                   aaCGAATCGGCCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttt   <  1:5835486/68‑1 (MQ=255)
                                                                    aCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAtttt  >  1:11980150/1‑68 (MQ=255)
                                                                    aCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAtttt  >  1:28898188/1‑68 (MQ=255)
                                                                    |                                                                   
CGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTT  >  CP009273/3227001‑3227135

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: