Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,703,146 C→A 6.1% A657A (GCC→GCA 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,1460CA 6.1% 851.3 / 9.4 380A657A (GCC→GCArsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (219/120);  new base A (3/24);  total (235/145)
Fisher's exact test for biased strand distribution p-value = 4.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 423 total aligned reads displayed.

TGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGG  >  CP009273/1703079‑1703213
                                                                   |                                                                   
tGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                     <  1:13060347/68‑1 (MQ=25)
tGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                     <  1:11253572/68‑1 (MQ=11)
   aaCAGCAACAGGCTAATGCCGAGCCAGAACACCAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                   <  1:1156269/67‑1 (MQ=255)
   aaCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                   <  1:6831342/67‑1 (MQ=40)
   aaCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                   <  1:17413441/67‑1 (MQ=40)
   aaCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                   <  1:13790664/67‑1 (MQ=40)
   aaCAGCAACAGGCAAATGCCGACCCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                   <  1:9415969/67‑1 (MQ=11)
    aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTa                                                                 >  1:35659142/1‑68 (MQ=255)
    aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                 >  1:28398464/1‑68 (MQ=40)
    aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                 >  1:26755965/1‑68 (MQ=40)
    aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTa                                                                 >  1:34180663/1‑68 (MQ=255)
    aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTAGAAGCCGCTa                                                                 >  1:5113744/1‑68 (MQ=14)
    aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTAGAAGCCGCTa                                                                 >  1:31282458/1‑68 (MQ=14)
    aCAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                  <  1:32095298/67‑1 (MQ=25)
        cAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                              >  1:29787257/1‑67 (MQ=40)
         aaCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGCCGCGGCTATTGc                                                             >  1:18172974/1‑67 (MQ=14)
         aaCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGACGAAGCGGCTATTGc                                                             >  1:12561449/1‑67 (MQ=14)
          aCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCAGTCGAAGCCGCTATTGccc                                                           >  1:10358749/1‑68 (MQ=255)
             ggCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGa                                                         >  1:26127561/1‑67 (MQ=255)
              gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCgcg                                                        >  1:28139395/1‑65 (MQ=255)
               cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                         <  1:22136227/65‑1 (MQ=255)
               cTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCgcgc                                                       >  1:15040159/1‑67 (MQ=14)
                 aaTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCa                                                     >  1:6284276/1‑67 (MQ=255)
                    gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                 >  1:29147877/1‑68 (MQ=255)
                    gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgc                                                 >  1:20361392/1‑68 (MQ=255)
                    gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                 >  1:28847704/1‑68 (MQ=255)
                    gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                 >  1:12379502/1‑68 (MQ=255)
                    gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                 >  1:22634205/1‑68 (MQ=255)
                    gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCAAAAg                                                  >  1:33558283/1‑67 (MQ=14)
                     ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                 >  1:18893853/1‑67 (MQ=255)
                      cGATCCAGAACACCAGGGCGCTCCGGGAAAAGCGTACTTCTTAGAGGCTATTGCCCGAGCCAAAgcg                                                <  1:18591959/67‑1 (MQ=255)
                      cGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAgcgc                                               >  1:6503449/1‑68 (MQ=255)
                         accAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCa                                            >  1:7391130/2‑68 (MQ=18)
                          ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                           >  1:17959961/1‑68 (MQ=35)
                          ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaa                                           >  1:2510286/1‑68 (MQ=35)
                           cAGCACAACAGGTCGATCCCCGGAAAGCTGCCGTCGAAGCCGCTATGGCCCGAGCCAAAGCGCGCaaa                                          <  1:21210090/68‑1 (MQ=11)
                           cAGCACAACAGGGCGATCCGCGCAAACCCGCCGGCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                          <  1:16570545/68‑1 (MQ=11)
                           cAGAACCACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                          <  1:27858233/68‑1 (MQ=255)
                           cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                          <  1:15448243/68‑1 (MQ=37)
                               acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                      >  1:5839750/1‑68 (MQ=255)
                               acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACggg                                      >  1:30349641/1‑68 (MQ=255)
                               acaacaGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACg                                        <  1:11979207/66‑1 (MQ=255)
                                caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                     >  1:8233923/1‑68 (MQ=255)
                                 aacaGGTCGATCCGCGCAAAGCCGCAGTCGCAGCGGCTATTGCCCGAGCGAAAGCGCGAAAGCgggc                                     >  1:16426584/1‑66 (MQ=11)
                                  acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGAGCAAACGGGAAc                                   >  1:23837428/1‑68 (MQ=255)
                                   caGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                  >  1:10693916/1‑68 (MQ=255)
                                    aGGTCGATCTGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                  >  1:34945376/1‑67 (MQ=255)
                                    aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                  >  1:17101200/1‑67 (MQ=255)
                                    aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                  >  1:10477856/1‑67 (MQ=255)
                                    aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                  >  1:2002864/1‑67 (MQ=255)
                                      gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCa                               >  1:27051613/1‑68 (MQ=255)
                                       tCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGATCAGCa                               >  1:32430780/1‑67 (MQ=255)
                                       tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                <  1:5075039/66‑1 (MQ=255)
                                       tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                               >  1:2593931/1‑67 (MQ=255)
                                       tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCa                               <  1:13896499/67‑1 (MQ=255)
                                         gATCCGCGCAAAGCCGCCGTCGCCGCCGCTATTGCCCGGGCGAAGGCGCGCAAGCGGGACCAGCagca                            >  1:24582672/1‑65 (MQ=1)
                                           tCCTCGTAAACCTGCCGTGGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                          <  1:2479353/68‑1 (MQ=14)
                                           tCCTCGCAAACCTGCCGTCGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                          <  1:31850769/68‑1 (MQ=14)
                                           tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCg                           >  1:27566515/1‑67 (MQ=255)
                                           tCCGCGCAAAGCCGCCGTCGAAGCCGCTATAGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                          <  1:22834021/68‑1 (MQ=255)
                                           tCCGCGCAAAGCCGCCGTCGAAGCCGCTATAGCACCAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                          <  1:24107958/68‑1 (MQ=17)
                                           tCCGCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                          <  1:4785616/68‑1 (MQ=255)
                                           tCCCCGTAAACCCGCCGTAGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCg                           <  1:2432120/67‑1 (MQ=14)
                                           tCCACGCAAATCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCgg                          <  1:20599902/68‑1 (MQ=17)
                                            ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                         <  1:30836829/68‑1 (MQ=255)
                                             cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCGGGACCCGCAACCGGCt                        >  1:20243337/1‑68 (MQ=14)
                                               cgcAAAGCCGCCGGCGAAGCCGCTATTGCCCTAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTaa                      <  1:1625018/68‑1 (MQ=14)
                                                gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                     >  1:34337396/1‑68 (MQ=255)
                                                gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                     >  1:4295210/1‑68 (MQ=255)
                                                gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                     >  1:16927098/1‑68 (MQ=255)
                                                gaaaaCCCGCCGTAGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                     <  1:25658963/66‑1 (MQ=17)
                                                 cAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATg                    <  1:15765008/68‑1 (MQ=255)
                                                  aaaCCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATg                    <  1:33507703/67‑1 (MQ=255)
                                                  aaaCCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATg                    <  1:35272396/67‑1 (MQ=255)
                                                     gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                >  1:8833736/1‑68 (MQ=255)
                                                     gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                >  1:5277547/1‑68 (MQ=255)
                                                     gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                >  1:21488705/1‑68 (MQ=255)
                                                     gccgccGTCGAAGCCGCTATAGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                >  1:36500520/1‑68 (MQ=255)
                                                      ctgccgTTGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg               <  1:26493526/66‑1 (MQ=11)
                                                      ctgccgTTGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAg               <  1:5704506/66‑1 (MQ=11)
                                                       cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                <  1:14474112/66‑1 (MQ=255)
                                                          cgTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGc              <  1:9835011/65‑1 (MQ=255)
                                                          cgTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCa            <  1:21109163/67‑1 (MQ=255)
                                                           gtggAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga          <  1:29992173/65‑1 (MQ=255)
                                                            tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaa         >  1:20847980/1‑68 (MQ=255)
                                                             cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaa         <  1:37067588/67‑1 (MQ=255)
                                                             cGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag        <  1:9147261/68‑1 (MQ=255)
                                                              gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga       <  1:19553429/68‑1 (MQ=255)
                                                              gAAGCCGCTATCGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCagaaga       <  1:3787043/68‑1 (MQ=255)
                                                                aGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc     <  1:10060301/68‑1 (MQ=255)
                                                                  ccGCTATTGCCCTTGCCAAAGCGCGCAAGCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGCACAg   >  1:32657240/1‑68 (MQ=11)
                                                                  ccGCTATTGCCCGAGCCAAAGCGCGCACACGGGAACAGCACCCGGCTAATGCGGCGCAAGAAGAACAg   >  1:8346779/1‑68 (MQ=18)
                                                                   cGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg   <  1:25443299/67‑1 (MQ=255)
                                                                   cGCTATTGCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg   <  1:14670276/67‑1 (MQ=255)
                                                                   cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAgg  <  1:16064376/68‑1 (MQ=255)
                                                                   cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAgg  >  1:7645902/1‑68 (MQ=255)
                                                                   cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg   <  1:23298956/67‑1 (MQ=255)
                                                                   cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg   <  1:31023542/67‑1 (MQ=255)
                                                                   cGCTATTGCACGAGCCAAAGCGCGAAATCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg   <  1:15204078/67‑1 (MQ=14)
                                                                   cGCTATAGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg   <  1:17630664/67‑1 (MQ=255)
                                                                   |                                                                   
TGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGG  >  CP009273/1703079‑1703213

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: