Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,703,158 A→C 14.2% R661R (CGA→CGC rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,703,1580AC14.2% 565.8 / 28.7 315R661R (CGA→CGCrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (153/108);  new base C (44/6);  total (197/118)
Fisher's exact test for biased strand distribution p-value = 4.45e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGC  >  CP009273/1703093‑1703222
                                                                 |                                                                
gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTGGACGCCGCTATTGCCCGAGc                                                                >  1:24562991/1‑68 (MQ=255)
gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAg                                                                 >  1:10036643/1‑67 (MQ=255)
gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAg                                                                 >  1:9397208/1‑67 (MQ=255)
 cTAATTCCGACGCAGAACAACTGGTGGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:24273206/65‑1 (MQ=11)
 cTAATGCCGAGCCAGAACAAGAGGGCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:7305694/65‑1 (MQ=14)
 cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:25493418/65‑1 (MQ=255)
 cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:7580751/65‑1 (MQ=255)
 cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGca                                                               >  1:23115352/1‑67 (MQ=255)
 cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCGGCTATTGCCCGCGcc                                                               >  1:9135234/1‑68 (MQ=14)
 cTAATCCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:17189701/65‑1 (MQ=255)
 cTAATCCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:5182405/65‑1 (MQ=255)
 cTAATCACGACCCAGAACACCATGGAGATCCCCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAg                                                                 <  1:26034585/66‑1 (MQ=2)
 aaaaTCCCGAGCCAGAACAACATGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                  <  1:6476075/63‑1 (MQ=14)
     tGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCaac                                                            >  1:29975070/1‑66 (MQ=255)
     tGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAg                                                           <  1:31178041/68‑1 (MQ=255)
      gCCGAGCCAGACCAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:27780097/1‑68 (MQ=14)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:22800803/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:20866169/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:32014942/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:436091/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:25324654/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                          >  1:13389443/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAg                                                           >  1:2211027/1‑67 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAg                                                           >  1:2509852/1‑67 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                          >  1:1383853/1‑68 (MQ=255)
      gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAg                                                           >  1:8423611/1‑67 (MQ=255)
          aGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgcgc                                                      >  1:1355267/1‑68 (MQ=39)
            ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                                    >  1:6115141/1‑68 (MQ=35)
            ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCGAGCaa                                                    >  1:33661591/1‑68 (MQ=255)
            ccAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                                    >  1:16036426/1‑68 (MQ=255)
             cAGAACACCATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                   <  1:32144420/68‑1 (MQ=255)
             cAGAACAACATGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                   <  1:18784433/68‑1 (MQ=255)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                   <  1:17720062/68‑1 (MQ=37)
             cAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                   <  1:29458284/68‑1 (MQ=255)
             cAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                   <  1:26904442/68‑1 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:21368468/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:4943478/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:32783368/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:22378526/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:21751342/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:18302820/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                               >  1:6305798/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACg                                                 <  1:4672809/66‑1 (MQ=40)
                 acaacaGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCggg                                               >  1:20228496/1‑68 (MQ=255)
                      aGGTCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAACCGGGAAcg                                           >  1:15097265/1‑66 (MQ=11)
                      aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                           >  1:20585902/1‑67 (MQ=255)
                      aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                           >  1:5857648/1‑67 (MQ=255)
                      aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                           >  1:16860827/1‑67 (MQ=255)
                        gtagATCCGCGAAAATCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                                        <  1:32283424/65‑1 (MQ=14)
                        gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                                        >  1:22561365/1‑68 (MQ=255)
                         tCGATCCGCGCAAAGCCGCAGTAGCCGCCGCTATTGCCCGCGCAAATGCGCGCAAACGCGAACAGCa                                        >  1:19981320/1‑67 (MQ=14)
                            aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCg                                    >  1:13700696/1‑68 (MQ=255)
                             tCCTCGTAAACCTGCCGTGGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                   <  1:14217953/68‑1 (MQ=14)
                             tCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                   <  1:2907692/68‑1 (MQ=255)
                             tCCTCGCAAACCCGCCGTCGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                   <  1:28023327/68‑1 (MQ=14)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                   <  1:30774740/68‑1 (MQ=255)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                   <  1:16147060/68‑1 (MQ=255)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                   <  1:18018631/68‑1 (MQ=255)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCg                                    <  1:799062/67‑1 (MQ=255)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCg                                    <  1:15334837/67‑1 (MQ=255)
                                  gcAAAGCCGCCGTCGACGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGACCAGCAACCGGCTAAt                              >  1:26152425/1‑68 (MQ=255)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTatt                              >  1:30335686/1‑66 (MQ=255)
                                   cAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGCACAGCAACCGGCTagt                              >  1:33046105/1‑65 (MQ=255)
                                       gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                         >  1:27180997/1‑68 (MQ=255)
                                       gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                         >  1:12523879/1‑68 (MQ=255)
                                       gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                         >  1:12150739/1‑68 (MQ=255)
                                        ccgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                         >  1:9687343/1‑67 (MQ=255)
                                         cgcagTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGc                       <  1:7046089/68‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCag                    >  1:5633009/1‑68 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCa                     <  1:16768016/67‑1 (MQ=255)
                                            cgTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCa                     <  1:10438551/67‑1 (MQ=255)
                                             gTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                   <  1:9838441/68‑1 (MQ=255)
                                              tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaa                  >  1:19741988/1‑68 (MQ=255)
                                               ggAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                 <  1:34406013/67‑1 (MQ=255)
                                               ggAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                 <  1:3996831/67‑1 (MQ=255)
                                               cGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCagaag                 <  1:28190303/68‑1 (MQ=255)
                                               cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaa                  <  1:14450492/67‑1 (MQ=255)
                                                gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGTAACAGCAACCGGCTAATGCGGAGCCagaaga                <  1:32510003/68‑1 (MQ=255)
                                                gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                <  1:28331582/68‑1 (MQ=255)
                                                gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                <  1:15849329/68‑1 (MQ=255)
                                                gAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                <  1:9711270/68‑1 (MQ=255)
                                                gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                <  1:28661533/68‑1 (MQ=14)
                                                  aGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc              >  1:14822982/1‑68 (MQ=255)
                                                  aGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc              >  1:5046683/1‑68 (MQ=255)
                                                    ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg            >  1:23821416/1‑68 (MQ=255)
                                                    ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg            >  1:16396063/1‑68 (MQ=255)
                                                    ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg            >  1:29206990/1‑68 (MQ=255)
                                                     ggCTATTGCCCGTGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg            <  1:8679546/66‑1 (MQ=14)
                                                     cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg            <  1:14315212/67‑1 (MQ=255)
                                                     agCTATTGCCCGCGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg            <  1:30951900/66‑1 (MQ=14)
                                                        tATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg        >  1:11773280/1‑68 (MQ=255)
                                                         aTTGCCCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGa       <  1:26400937/68‑1 (MQ=255)
                                                           tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGAt      >  1:10813649/1‑67 (MQ=255)
                                                             cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc    >  1:18238168/1‑67 (MQ=255)
                                                             cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc    >  1:27877903/1‑67 (MQ=255)
                                                             cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc    >  1:11597838/1‑67 (MQ=255)
                                                              ccGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc  <  1:15191797/68‑1 (MQ=255)
                                                              ccGAGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc  <  1:10108792/68‑1 (MQ=255)
                                                              ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc  >  1:361670/1‑68 (MQ=255)
                                                              ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc  >  1:31484618/1‑68 (MQ=255)
                                                                 |                                                                
GCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGC  >  CP009273/1703093‑1703222

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: