Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,703,036 G→C 5.0% A621P (GCC→CCC)  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,0360GC 5.0% 1271.5 / 8.0 515A621P (GCC→CCC) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base G (285/191);  new base C (2/27);  total (288/227)
Fisher's exact test for biased strand distribution p-value = 8.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 589 total aligned reads displayed.

AAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGA  >  CP009273/1702970‑1703103
                                                                  |                                                                   
aaaGCGCGCAAGCTTGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgc                                                                    <  1:37045468/68‑1 (MQ=255)
aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgc                                                                    <  1:20380650/68‑1 (MQ=35)
  aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGATCCGCGCAAAgcg                                                                   >  1:10069752/1‑66 (MQ=255)
  aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGCTCCGCGCAAAgcc                                                                   >  1:52281247/1‑67 (MQ=25)
  aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcg                                                                   >  1:23936509/1‑66 (MQ=25)
  aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                   >  1:48511092/1‑67 (MQ=35)
  aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                   >  1:4682057/1‑67 (MQ=35)
  aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgcc                                                                   >  1:42136099/1‑67 (MQ=25)
   gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                  >  1:27064063/1‑67 (MQ=35)
   gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgccgc                                                                 >  1:34727456/1‑68 (MQ=25)
       gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGt                                                              <  1:35114599/67‑1 (MQ=35)
       gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGta                                                             >  1:41851490/1‑67 (MQ=11)
             tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                       <  1:4583007/68‑1 (MQ=35)
             tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                       <  1:29676890/68‑1 (MQ=35)
                aaCAGCAACAGGCTAATGCGGAACCAGAAGAACATGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:25887897/67‑1 (MQ=255)
                aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:53461048/67‑1 (MQ=35)
                aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:38494004/67‑1 (MQ=35)
                 aCAGCAACAGGCTAATGCGTATCGATAAGAACCGTTAGCTCCGCGCATAGCCGACGGCGAAGCCGCTa                                                   >  1:47019740/1‑68 (MQ=2)
                 aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTAGCCGCCGCTa                                                   >  1:31164625/1‑68 (MQ=11)
                 aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTa                                                   >  1:34092388/1‑68 (MQ=255)
                 aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                   >  1:23194843/1‑68 (MQ=35)
                 aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGACGCCGCTa                                                   >  1:58108093/1‑68 (MQ=255)
                 aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTa                                                   >  1:16362331/1‑68 (MQ=255)
                  cAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                  >  1:12564377/1‑68 (MQ=35)
                   aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTAtt                                                 >  1:52742302/1‑68 (MQ=255)
                      accAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                               <  1:54361361/65‑1 (MQ=25)
                         aGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCCGCTATTGCCCgc                                           >  1:28877020/1‑67 (MQ=11)
                         aGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                           >  1:12962734/1‑68 (MQ=35)
                           gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCTTg                                          >  1:33605525/1‑67 (MQ=25)
                           gCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTCGCCGCCGCTATTGCCCgc                                           >  1:5550550/1‑65 (MQ=11)
                           gCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCgag                                          >  1:727701/1‑65 (MQ=11)
                            cTAATGCGGAACCAGAAGACCAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGcc                                        >  1:46447869/1‑68 (MQ=11)
                            cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTAGAAGCCGCTATTGCCCGAGc                                         >  1:19028181/1‑67 (MQ=11)
                            cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                           <  1:6281376/65‑1 (MQ=35)
                            cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                          <  1:18859892/66‑1 (MQ=35)
                            cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                          <  1:14526560/66‑1 (MQ=35)
                            cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                          <  1:57188776/66‑1 (MQ=35)
                            cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGcc                                        >  1:43437968/1‑68 (MQ=25)
                              aaTGCGGAACCAGAAGAACAGCTCGATCCGCGCAAAGCCGCCGTGGACGCCGCTATTGCCCGAGCCaa                                      >  1:49990176/1‑68 (MQ=11)
                                tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                     >  1:47738475/1‑67 (MQ=25)
                                tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTCGACGCCGCTATTGCCCGGGCCaaa                                     >  1:32120841/1‑67 (MQ=11)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTAGCAGCCGCTATTGCCCGCGCaaagg                                    >  1:28275858/1‑65 (MQ=11)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTAGAAGCCGCTATTGCCCGAGCGAAAg                                    >  1:33069831/1‑67 (MQ=11)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCTTGCGAACGc                                   >  1:42859506/1‑68 (MQ=11)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                     >  1:30691092/1‑66 (MQ=25)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGc                                   >  1:37654407/1‑68 (MQ=25)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCCCAAg                                    >  1:34352218/1‑67 (MQ=11)
                                 gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGc                                   >  1:20988902/1‑68 (MQ=11)
                                      aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCa                              >  1:9147142/1‑68 (MQ=255)
                                      aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCCCGCa                              >  1:44461507/1‑68 (MQ=11)
                                      aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCACGCa                              >  1:11944952/1‑68 (MQ=255)
                                       ccAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTAGAAGCCGCTATTGCCCGGGCGAAGGCCCGCaa                             >  1:36539178/1‑68 (MQ=11)
                                       ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCGAAGGCCCGCaa                             >  1:10813377/1‑68 (MQ=11)
                                       ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCaa                             >  1:21713843/1‑68 (MQ=25)
                                        cAGAATAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:35892867/68‑1 (MQ=255)
                                        cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:31580661/68‑1 (MQ=40)
                                        cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:5623503/68‑1 (MQ=40)
                                           aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAACGCACGCAAGCTg                         >  1:33359767/1‑68 (MQ=255)
                                           aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTg                         >  1:11358741/1‑68 (MQ=255)
                                           aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTg                         >  1:40060844/1‑68 (MQ=255)
                                           aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCACGCAAGCTg                         >  1:51769290/1‑68 (MQ=255)
                                           aagaaCAGGTCGATCCGCGCAAAGCCGCCGTAGCAGCCGCTATTGCCCGCGCCAAAGCCAGCAAGcgg                         >  1:24942280/1‑66 (MQ=11)
                                           aagaaCAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTg                         >  1:38992917/1‑68 (MQ=14)
                                             gaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGAGCCAAAGCACGCAAGCTGGa                       >  1:25350019/1‑68 (MQ=255)
                                             caaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGa                       <  1:49317489/67‑1 (MQ=18)
                                              aaCAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCGGCTATTGCCCGGGCCAAGGCACGCAAGCGGGa                       >  1:50562413/1‑67 (MQ=11)
                                              aaCAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGGGCCAAAGCACGCAACCTCGa                       >  1:18193302/1‑67 (MQ=11)
                                               tcAGTTCGATCCGCGCAAAGCCGACGTCGAAGCCGCTATCGCCCGTGCCAAAGCACGCAAGCTGGa                       <  1:15645175/65‑1 (MQ=11)
                                               ccAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAc                     <  1:6872612/67‑1 (MQ=255)
                                               aCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAAc                     >  1:32568774/1‑68 (MQ=25)
                                                cAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                    <  1:24242078/68‑1 (MQ=25)
                                                    tCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                 >  1:8431067/1‑67 (MQ=25)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                    <  1:43538238/64‑1 (MQ=35)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGcc                 >  1:15232233/1‑66 (MQ=25)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGc                  <  1:797600/66‑1 (MQ=35)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGc                  <  1:18651931/66‑1 (MQ=35)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                 >  1:36241955/1‑67 (MQ=35)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                 >  1:41211136/1‑67 (MQ=35)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGTCCAAAGAACGCAAGCGGGACCCGCa                 >  1:51485195/1‑67 (MQ=11)
                                                    tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCa                 >  1:23582150/1‑67 (MQ=25)
                                                    gcGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGc                  <  1:33902139/65‑1 (MQ=25)
                                                        tCCTCGTAAAGCAGCCGTCGAAGCCGCTATGGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:13263234/68‑1 (MQ=11)
                                                        tCCGCGTAAACCCGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:35488620/68‑1 (MQ=11)
                                                        tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:50916571/68‑1 (MQ=35)
                                                        tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:384058/68‑1 (MQ=35)
                                                        tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:22132035/68‑1 (MQ=35)
                                                        tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:44000555/68‑1 (MQ=35)
                                                        tCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:26314279/68‑1 (MQ=255)
                                                        tCCGCGCAAACCCGCCGTGGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAg             <  1:30154826/67‑1 (MQ=255)
                                                        tCCCCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:17198950/68‑1 (MQ=25)
                                                         ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAggg           >  1:53746601/1‑67 (MQ=25)
                                                         ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAg             >  1:13726060/1‑66 (MQ=35)
                                                          cgcgCAAAGCCGCCGTCGCAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt          >  1:26740593/1‑68 (MQ=255)
                                                             gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAAt       >  1:45249499/1‑68 (MQ=35)
                                                               aaaCCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATg      <  1:19601291/67‑1 (MQ=25)
                                                               aaaCCCGCCGTCGAAGCCGCTATCGCCCGGGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATg      <  1:20108005/67‑1 (MQ=11)
                                                                aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc    >  1:48922153/1‑68 (MQ=37)
                                                                  gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:22699094/1‑68 (MQ=40)
                                                                  gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:56900096/1‑68 (MQ=40)
                                                                  gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:48142197/1‑68 (MQ=255)
                                                                  |                                                                   
AAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGA  >  CP009273/1702970‑1703103

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: