Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,702,966 T→A 8.6% 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,9660TA 8.6% 840.3 / 9.2 422R597R (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 (122/231);  new base A (28/2);  total (171/251)
Fisher's exact test for biased strand distribution p-value = 1.45e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 451 total aligned reads displayed.

AGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGC  >  CP009273/1702899‑1703032
                                                                   |                                                                  
aGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGt                                                                    <  1:11054961/68‑1 (MQ=255)
   cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                    <  1:21719088/65‑1 (MQ=35)
   cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                    <  1:20246450/65‑1 (MQ=35)
   cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                    <  1:3686781/65‑1 (MQ=35)
   cTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                    <  1:559810/65‑1 (MQ=35)
   cTAATGCGGAACCAGAACAACAGGTCGAACCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                    <  1:21064048/65‑1 (MQ=25)
   cTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTAGAAGCGGCTATTGCCCGCGcc                                                                 >  1:6784221/1‑68 (MQ=11)
   cTAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                    <  1:38198984/65‑1 (MQ=25)
     aaTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCgcgcac                                                                >  1:34936387/1‑65 (MQ=11)
        gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTGGCAGCCGCTATTGCCCGCGCCaacg                                                             >  1:2614606/1‑65 (MQ=11)
        gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCGGCTATTGCCCGCGCCaacg                                                             >  1:10310750/1‑65 (MQ=11)
        gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGAGCCaatg                                                             >  1:34197626/1‑65 (MQ=11)
        gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgc                                                            >  1:36280322/1‑68 (MQ=35)
        gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCGGCTATTGCCCGCGCCAAAg                                                             >  1:28312593/1‑67 (MQ=11)
        gCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTAATGCCCGAGCTaacg                                                             >  1:20170015/1‑65 (MQ=11)
        gCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCAAAAg                                                             >  1:13584969/1‑67 (MQ=11)
          ggAACCAGAACCAGAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAgcac                                                          <  1:34710920/68‑3 (MQ=11)
          ggAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                          >  1:37381700/1‑68 (MQ=14)
             aCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                       >  1:33263682/1‑68 (MQ=37)
             aCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                       <  1:21367374/68‑1 (MQ=37)
              ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCaa                                                      >  1:12388627/1‑68 (MQ=35)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                     <  1:39128946/68‑1 (MQ=37)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                     <  1:15860328/68‑1 (MQ=37)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                     <  1:2377576/68‑1 (MQ=37)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                     <  1:22716450/68‑1 (MQ=37)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                     <  1:15026270/68‑1 (MQ=37)
               cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                     <  1:20664091/68‑1 (MQ=37)
                   acaacaGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGCGCCAAAGCCCACAAGcgcg                                                 >  1:10098504/1‑65 (MQ=11)
                   acaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTgg                                                 >  1:29556588/1‑68 (MQ=18)
                    caacaTGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                <  1:5129745/68‑1 (MQ=255)
                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                >  1:11864374/1‑67 (MQ=35)
                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                >  1:16333529/1‑67 (MQ=35)
                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGa                                                >  1:26417407/1‑67 (MQ=25)
                     aacaGGTCGATCCGCGCAAAGCCGCAGTCGCCGCCGCTATTGCCCGCGCCCACGCGAGCAAGGTggc                                                >  1:4872572/1‑66 (MQ=2)
                      acaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAAc                                              >  1:17945712/1‑68 (MQ=25)
                         ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGc                                           >  1:31190714/1‑68 (MQ=25)
                          gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAcctc                                           >  1:7924991/1‑64 (MQ=25)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa                                             <  1:18362110/64‑1 (MQ=35)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                           <  1:22497747/66‑1 (MQ=35)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                           <  1:14317614/66‑1 (MQ=35)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                           <  1:24478138/66‑1 (MQ=35)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                          <  1:9316871/67‑1 (MQ=35)
                           tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                          <  1:10752409/67‑1 (MQ=35)
                           tCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCa                                          <  1:28914440/67‑1 (MQ=11)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:9685774/68‑1 (MQ=35)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:34361629/68‑1 (MQ=35)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:38704364/68‑1 (MQ=35)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:33880623/68‑1 (MQ=35)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:36938009/68‑1 (MQ=35)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:14013789/68‑1 (MQ=35)
                               tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     <  1:24835126/68‑1 (MQ=35)
                                ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg                                      >  1:29903140/1‑66 (MQ=35)
                                ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAACAg                                      >  1:7653741/1‑66 (MQ=25)
                                 cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCt                                   >  1:18759817/1‑68 (MQ=35)
                                   cgaaaaTCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                 <  1:17450589/65‑1 (MQ=255)
                                    gcAAAGCCGCCGTGGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                                <  1:32299165/68‑1 (MQ=255)
                                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAAt                                >  1:19742588/1‑68 (MQ=25)
                                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGAAAGCTGGACCAGCAACAGGCTAAt                                >  1:31986457/1‑68 (MQ=11)
                                    gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                                >  1:4393931/1‑68 (MQ=25)
                                    gcAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAt                                <  1:2195729/68‑1 (MQ=255)
                                    gcAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGGGGGCACAGCAGCAGGCTAAt                                >  1:24380177/1‑68 (MQ=11)
                                          ccgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGaa                          <  1:24756299/68‑1 (MQ=255)
                                           cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGaa                          <  1:27218465/67‑1 (MQ=255)
                                           cgccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                         <  1:32997428/68‑1 (MQ=255)
                                              cgTTGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:17742673/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:12991416/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:9210395/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:7310295/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:19947641/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:860895/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:4133944/67‑1 (MQ=255)
                                              cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:25626303/67‑1 (MQ=255)
                                              cgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:6279439/67‑1 (MQ=255)
                                              cgTCGAAGCAGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                       <  1:23076499/67‑1 (MQ=255)
                                               gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga                     <  1:35076087/68‑1 (MQ=255)
                                                 cGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                   <  1:25161177/68‑1 (MQ=255)
                                                  gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga                  <  1:10625041/68‑1 (MQ=255)
                                                  gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaaga                  <  1:15437456/68‑1 (MQ=255)
                                                      ccGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAAGAACAg              >  1:36033812/1‑68 (MQ=255)
                                                       cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg              <  1:28441619/67‑1 (MQ=255)
                                                       cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg              <  1:11357290/67‑1 (MQ=255)
                                                        gCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt            <  1:23938212/68‑1 (MQ=255)
                                                          tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGtag          >  1:19324499/1‑66 (MQ=255)
                                                           aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGTCTAATGCGGAACCAGAAGAACAGGTCGa         <  1:17045748/68‑1 (MQ=255)
                                                           aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa         <  1:30675317/68‑1 (MQ=255)
                                                           aTTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGa         <  1:27680435/68‑1 (MQ=25)
                                                            ttGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGa         >  1:28023135/1‑67 (MQ=255)
                                                             tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGAt        >  1:32717799/1‑67 (MQ=255)
                                                             tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGAt        >  1:33618370/1‑67 (MQ=255)
                                                               cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                     >  1:29262520/1‑36 (MQ=40)
                                                               cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGaaac      >  1:38089857/1‑64 (MQ=255)
                                                               cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGATcc      >  1:12053227/1‑67 (MQ=255)
                                                               cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc      >  1:6643006/1‑67 (MQ=40)
                                                               cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATcc      >  1:16641651/1‑67 (MQ=255)
                                                                ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATcctc    >  1:14701518/1‑66 (MQ=255)
                                                                ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATcctc    >  1:31751141/1‑66 (MQ=255)
                                                                ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc    >  1:37223377/1‑68 (MQ=40)
                                                                ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc    <  1:19562696/68‑1 (MQ=40)
                                                                ccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCcgc    >  1:13325237/1‑68 (MQ=255)
                                                                cagTGCCATAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc    <  1:35584718/66‑1 (MQ=255)
                                                                   |                                                                  
AGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGC  >  CP009273/1702899‑1703032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: