Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,702,966 T→A 11.3% R597R (CGT→CGA 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,702,9660TA11.3% 679.9 / 12.1 356R597R (CGT→CGArsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (121/170);  new base A (30/4);  total (171/185)
Fisher's exact test for biased strand distribution p-value = 1.89e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 388 total aligned reads displayed.

GCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGC  >  CP009273/1702901‑1703030
                                                                 |                                                                
 cTAATGCGGACCCAGAACACGTTGCCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:28437751/65‑1 (MQ=11)
 cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:12045959/65‑1 (MQ=35)
 cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:125508/65‑1 (MQ=35)
 cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:1477082/65‑1 (MQ=35)
 cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:32784364/65‑1 (MQ=35)
 cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGGCGAAGCCGCTATTGCCCGt                                                                  <  1:7989148/65‑1 (MQ=255)
 cTAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:17648494/65‑1 (MQ=25)
 cTAATCCGGAACCAGAACAACAGGTAGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                  <  1:7767942/65‑1 (MQ=11)
    ctGCCGAACCAGAACAACAGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                            <  1:31073306/67‑1 (MQ=11)
      gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTGGCAGCCGCTATTGCCCGAGCCAAAg                                                           >  1:31606245/1‑67 (MQ=11)
      gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAg                                                           >  1:2276520/1‑67 (MQ=25)
      gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGc                                                                >  1:20101519/1‑62 (MQ=25)
      gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCaacg                                                           >  1:8573521/1‑65 (MQ=11)
      gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGACGAAGCCGCTATTGCCCGTGCCAAAg                                                           >  1:6435626/1‑67 (MQ=255)
      gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCACGAGCCAAAg                                                           >  1:19111205/1‑67 (MQ=11)
        ggAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgca                                                         <  1:1588681/67‑2 (MQ=16)
           aCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCCGCTATTGCCCGCGCCAAGGCCCGCa                                                     >  1:16863855/1‑68 (MQ=11)
             cAGAACAAGAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTTCCCGTGCCAAAGCGCGCAAg                                                   <  1:26687019/68‑1 (MQ=255)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                   <  1:20558712/68‑1 (MQ=37)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                   <  1:26602988/68‑1 (MQ=37)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                   <  1:1026441/68‑1 (MQ=37)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                   <  1:7002767/68‑1 (MQ=37)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                   <  1:9361161/68‑1 (MQ=37)
             cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                   <  1:6201871/68‑1 (MQ=37)
                aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                >  1:27713282/1‑68 (MQ=40)
                aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGCAAGCTg                                                >  1:2210391/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCt                                                 <  1:24630748/66‑1 (MQ=39)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCt                                                 <  1:3059643/66‑1 (MQ=39)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTgg                                               >  1:32208780/1‑68 (MQ=255)
                 acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTgg                                               >  1:15528584/1‑68 (MQ=18)
                  caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                              <  1:4156493/68‑1 (MQ=37)
                  caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                              >  1:14353622/1‑68 (MQ=37)
                   aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                             <  1:21757204/68‑1 (MQ=35)
                   aacaGGGCGATCCGCGCAAAGCTGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                             <  1:2302630/68‑1 (MQ=255)
                    acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCGAAAGCGAGCAAGCTGGCAc                                            >  1:17820382/1‑68 (MQ=11)
                      aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa                                           >  1:29921488/1‑67 (MQ=35)
                      aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa                                           >  1:29808825/1‑67 (MQ=35)
                        gTCGATCCGCGCAAAGCCGCCGTGGAAGCCGCTATTGCCCGAGCGAAAGCGCGAAAGCTGGCACAGCa                                        >  1:8182138/1‑68 (MQ=11)
                         tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                        >  1:26953886/1‑67 (MQ=35)
                         tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCa                                        >  1:25022263/1‑67 (MQ=25)
                             tCCGGGCAAAGCCGCCGTCGAAGCCGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                   <  1:34422804/68‑1 (MQ=255)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                   <  1:11564426/68‑1 (MQ=35)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                   <  1:8862879/68‑1 (MQ=35)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                   <  1:24419288/68‑1 (MQ=35)
                             tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                   <  1:25587760/68‑1 (MQ=35)
                              ccGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAggg                                  >  1:14560224/1‑67 (MQ=11)
                              ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGc                                  <  1:12987935/68‑1 (MQ=35)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                               <  1:28561491/67‑1 (MQ=35)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                              >  1:11638847/1‑68 (MQ=35)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                              <  1:11789260/68‑1 (MQ=35)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                              <  1:33493678/68‑1 (MQ=35)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                              <  1:3923756/68‑1 (MQ=35)
                                  gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                              >  1:2583654/1‑68 (MQ=25)
                                   cAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                              >  1:28719198/1‑67 (MQ=255)
                                    aaaGCCGCCGTCGAAGCCGCTATGGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATg                             <  1:25943771/67‑1 (MQ=255)
                                     aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCGTAAGCTCGCAAGCTGTAACAGCAACAGGCTAATGCg                           >  1:19992438/1‑68 (MQ=11)
                                       tccgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                         <  1:23425920/67‑1 (MQ=255)
                                       gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCggt                         >  1:10945322/1‑67 (MQ=39)
                                       gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCggt                         >  1:28017387/1‑67 (MQ=39)
                                       gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCgg                          >  1:10051936/1‑67 (MQ=39)
                                       gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                         >  1:23183967/1‑68 (MQ=40)
                                            cggcGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:21031711/64‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCag                    >  1:18244808/1‑68 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:33005216/67‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:15393601/67‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:5028765/67‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:17089085/67‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:33654877/67‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:30716429/67‑1 (MQ=255)
                                            cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:2543544/67‑1 (MQ=255)
                                            cgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                     <  1:13268488/67‑1 (MQ=255)
                                              tCGAAGCCGCTATTGCCCGAGCCAAAGCGCCCAAGCTGGAACAGCAACAGGCTAAGGAGGAGCCagaa                  >  1:18663060/1‑68 (MQ=11)
                                                gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga                <  1:27933029/68‑1 (MQ=255)
                                                gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga                <  1:19572711/68‑1 (MQ=255)
                                                    ccGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAg            >  1:9779048/1‑68 (MQ=255)
                                                    ccGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            >  1:27474039/1‑68 (MQ=255)
                                                     ctcTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            <  1:5578264/65‑1 (MQ=255)
                                                     cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGGGGAACCAGAAGAACAg            <  1:15060491/67‑1 (MQ=255)
                                                     cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            <  1:27207551/67‑1 (MQ=255)
                                                     cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            <  1:30390458/67‑1 (MQ=255)
                                                     agCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            <  1:14194460/66‑1 (MQ=255)
                                                     agCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            <  1:3333709/66‑1 (MQ=255)
                                                     agCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg            <  1:7507047/66‑1 (MQ=16)
                                                        tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt          <  1:33987114/66‑1 (MQ=255)
                                                         aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa       <  1:29000936/68‑1 (MQ=255)
                                                          ttGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt      <  1:24099377/68‑1 (MQ=255)
                                                           tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACATGTCgct      >  1:15999463/1‑65 (MQ=255)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc    >  1:4607316/1‑67 (MQ=40)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc    >  1:24751094/1‑67 (MQ=40)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc    >  1:19736120/1‑67 (MQ=40)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc    >  1:29494337/1‑67 (MQ=40)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATcc    >  1:2621271/1‑67 (MQ=255)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATcc    >  1:11374611/1‑67 (MQ=255)
                                                             cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATcc    >  1:30970779/1‑67 (MQ=255)
                                                              ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGTATCCAGAAGAACAGGAGGCTCcgc  >  1:13618858/1‑68 (MQ=14)
                                                              ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAgg           >  1:9610823/1‑59 (MQ=255)
                                                              ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCcgc  >  1:25137548/1‑68 (MQ=255)
                                                              ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc  >  1:1623540/1‑68 (MQ=40)
                                                              ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc  >  1:23626875/1‑68 (MQ=40)
                                                              ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCcgc  >  1:9476810/1‑68 (MQ=255)
                                                                 |                                                                
GCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGC  >  CP009273/1702901‑1703030

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: