Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,703,158 A→C 15.6% 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,1580AC15.6% 808.2 / 45.2 463R661R (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 (219/159);  new base C (56/23);  total (278/185)
Fisher's exact test for biased strand distribution p-value = 4.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 493 total aligned reads displayed.

AGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGC  >  CP009273/1703091‑1703224
                                                                   |                                                                  
  gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAg                                                                   >  1:34832409/1‑67 (MQ=255)
   cTAATTCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                    <  1:35047930/65‑1 (MQ=255)
   cTAATGCCGAGCCAGAACACCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                    <  1:20388812/65‑1 (MQ=255)
   cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCgcgca                                                                 >  1:7490425/1‑67 (MQ=14)
   cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGCCGCGGCTATTGCCCgcgca                                                                 >  1:38674443/1‑67 (MQ=11)
   cTAATGACGAGCCAGAACAACATGGCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGa                                                                    <  1:29723316/65‑1 (MQ=14)
   atAATTCCGACCCAGAACAACAGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAg                                                                   <  1:23330639/65‑1 (MQ=14)
   atAATGCCGACGCAGAACAAGAGGCGGATCCTCGCAAAGCAGCCGTCGAAGCAGCTATTGCCCGa                                                                    <  1:4381921/64‑1 (MQ=255)
   atAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                    <  1:10661835/64‑1 (MQ=255)
   atAATCCCGACCCAGAACACCAGGTGGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                    <  1:15753180/64‑1 (MQ=11)
        gacGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAACgc                                                            >  1:2947396/3‑68 (MQ=14)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                            >  1:22634000/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                            >  1:7989619/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                            >  1:7271202/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAg                                                             >  1:9244288/1‑67 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAACgc                                                            >  1:32367017/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                            >  1:26823261/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                            >  1:26667800/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                            >  1:2111811/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAg                                                             >  1:5281911/1‑67 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGACCGAGCCCAAg                                                             >  1:8522549/1‑67 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCAAAgc                                                            >  1:38451024/1‑68 (MQ=255)
        gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCGGCTATTGCCCGCGCCAAAg                                                             >  1:31324719/1‑67 (MQ=14)
         cggAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcc                                                           >  1:2447657/3‑67 (MQ=255)
          cGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgccc                                                          >  1:21607911/1‑66 (MQ=255)
             accAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                                       >  1:22914336/2‑68 (MQ=18)
              ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGaaa                                                      >  1:26406137/1‑68 (MQ=25)
               cAGAACATGAGGTCGATCCTCGGATAGGAGCCGTCTAACCGCCAATTGCCCGAGCCAAAGCGCGCaaa                                                     <  1:20001408/68‑1 (MQ=2)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                     <  1:4887820/68‑1 (MQ=37)
                   acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                                 >  1:11998015/1‑68 (MQ=255)
                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGa                                                >  1:17931756/1‑67 (MQ=25)
                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGa                                                >  1:9461026/1‑67 (MQ=25)
                     aacaGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAACGCCCGCAAGCgggc                                                >  1:28667200/1‑66 (MQ=11)
                      acaGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCGAAAGCCAGCAAGCGGGacc                                              >  1:28039530/1‑66 (MQ=11)
                       caGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                             <  1:1545544/68‑1 (MQ=255)
                        aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                             >  1:26216989/1‑67 (MQ=255)
                        aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                             >  1:24858260/1‑67 (MQ=255)
                        aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                             >  1:33108459/1‑67 (MQ=255)
                         ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                           >  1:2926819/1‑68 (MQ=255)
                          gtagATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCACCAAACGCGAACAGc                                           >  1:37348688/4‑67 (MQ=14)
                          gctgATCCTCGTAAACCAGCCGTTGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                                          <  1:25855190/65‑1 (MQ=1)
                          gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGACCAGCa                                          >  1:2997711/1‑68 (MQ=255)
                          gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                                          >  1:4049749/1‑68 (MQ=255)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCa                                          >  1:282041/1‑67 (MQ=25)
                           tCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCAAAAGCCCGCAAGCGGGCACAGCa                                          >  1:6155140/1‑67 (MQ=11)
                           tCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCACGCGCCAAAGCGCGCAAACGGGAACAGCa                                          >  1:28978226/1‑67 (MQ=14)
                             gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAAcc                                       >  1:5811331/1‑68 (MQ=255)
                               tCCTCGTAAACCTGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                     <  1:32131005/68‑1 (MQ=14)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAAcag                                      >  1:10006014/1‑65 (MQ=9)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                     <  1:14751597/68‑1 (MQ=255)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                                     <  1:24280665/68‑1 (MQ=255)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCg                                      >  1:35833675/1‑67 (MQ=255)
                               tCCACGAAAACCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAcccg                                     <  1:8752239/68‑3 (MQ=17)
                                ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAACCg                                      >  1:27383193/1‑66 (MQ=255)
                                ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                    >  1:34301893/1‑68 (MQ=255)
                                 cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCt                                   <  1:13785840/68‑1 (MQ=255)
                                 cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCt                                   >  1:33835980/1‑68 (MQ=255)
                                 cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCt                                   <  1:1520162/68‑1 (MQ=255)
                                 cgcgCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCCACGCGCGCAAACGGGAACAGCACCCGGCt                                   >  1:5519714/1‑68 (MQ=17)
                                   cgtaaaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTaa                                 <  1:14044209/65‑1 (MQ=255)
                                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                                >  1:30489860/1‑68 (MQ=255)
                                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                                >  1:23811943/1‑68 (MQ=255)
                                         gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCgg                            >  1:20775668/1‑67 (MQ=255)
                                         gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                           >  1:14489604/1‑68 (MQ=255)
                                          ctgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg                          <  1:35500632/66‑1 (MQ=255)
                                           tgccgTCGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGc                         <  1:30855468/67‑1 (MQ=14)
                                           cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg                          <  1:11690396/67‑1 (MQ=255)
                                           cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGc                         <  1:11474210/68‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCa                       <  1:22349038/67‑1 (MQ=255)
                                               gTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                     <  1:637597/68‑1 (MQ=255)
                                                 cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                   >  1:25379218/1‑68 (MQ=255)
                                                 agAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCTGGAACAGCAACCGGCTAATGCGGAGCCagaag                   <  1:18868801/67‑1 (MQ=14)
                                                  gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                  <  1:16103807/68‑1 (MQ=255)
                                                  gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                  <  1:181604/68‑1 (MQ=255)
                                                  gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                  <  1:23538432/68‑1 (MQ=255)
                                                  gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCCAATGCGGAGCCagaaga                  <  1:17449966/68‑1 (MQ=255)
                                                      ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              >  1:17147017/1‑68 (MQ=255)
                                                      ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              >  1:833077/1‑68 (MQ=255)
                                                      ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              >  1:9688224/1‑68 (MQ=255)
                                                       cGCTATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:19753363/67‑1 (MQ=255)
                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:33256460/67‑1 (MQ=255)
                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:6908898/67‑1 (MQ=255)
                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:9826923/67‑1 (MQ=255)
                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:3142731/67‑1 (MQ=255)
                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:34637009/67‑1 (MQ=255)
                                                       cGCTATAGCCCGCGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:37778281/67‑1 (MQ=14)
                                                       agCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:38184083/66‑1 (MQ=255)
                                                       agCTATCGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg              <  1:3181246/66‑1 (MQ=255)
                                                          tATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGt            <  1:16560855/66‑1 (MQ=255)
                                                               cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc      >  1:15910590/1‑67 (MQ=255)
                                                               cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc      >  1:22159914/1‑67 (MQ=255)
                                                               cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc      >  1:12718141/1‑67 (MQ=255)
                                                               cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc      >  1:1342817/1‑67 (MQ=255)
                                                                ccGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc    <  1:10718094/68‑1 (MQ=255)
                                                                ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc    >  1:32610936/1‑68 (MQ=255)
                                                                ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc    >  1:34116145/1‑68 (MQ=255)
                                                                ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc    >  1:19415423/1‑68 (MQ=255)
                                                                ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc    >  1:15017285/1‑68 (MQ=255)
                                                                ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcgc    >  1:26139145/1‑68 (MQ=255)
                                                                  gAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCgcgc  >  1:36575212/1‑68 (MQ=255)
                                                                   |                                                                  
AGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGC  >  CP009273/1703091‑1703224

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: