Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,703,036 G→C 5.6% 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.6% 887.1 / 5.7 357A621P (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 (194/138);  new base C (2/21);  total (197/160)
Fisher's exact test for biased strand distribution p-value = 2.64e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 398 total aligned reads displayed.

AGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGA  >  CP009273/1702972‑1703103
                                                                |                                                                   
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                   >  1:1891439/1‑67 (MQ=25)
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTGGATCCGCGCAAAgcg                                                                   >  1:31390065/1‑66 (MQ=255)
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcg                                                                   >  1:14209966/1‑66 (MQ=25)
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGAGCAAAgcc                                                                   >  1:40608332/1‑67 (MQ=25)
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgcg                                                                   >  1:28224733/1‑66 (MQ=255)
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgcc                                                                   >  1:25475946/1‑67 (MQ=25)
aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgcc                                                                   >  1:11611348/1‑67 (MQ=25)
 gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcgg                                                                  >  1:33843560/1‑65 (MQ=25)
 gcgcgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAgcggc                                                                 >  1:9723085/1‑68 (MQ=255)
           tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                         <  1:41176901/66‑1 (MQ=35)
           tGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                       <  1:41009667/68‑1 (MQ=35)
            ggAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                      <  1:9296996/68‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:42192670/67‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:2409302/67‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:27239886/67‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:27917870/67‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:19145881/67‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:11050196/67‑1 (MQ=35)
              aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:30470102/67‑1 (MQ=35)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTAGACGCCGCTa                                                   >  1:27288189/1‑68 (MQ=11)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTAGAAGCGGCTa                                                   >  1:2493768/1‑68 (MQ=11)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTGGACGCGGCTa                                                   >  1:33265085/1‑68 (MQ=11)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTa                                                   >  1:23779051/1‑68 (MQ=255)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                    <  1:20802062/67‑1 (MQ=35)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                   >  1:20307456/1‑68 (MQ=35)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                   >  1:10009076/1‑68 (MQ=35)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                   <  1:8786003/68‑1 (MQ=35)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTa                                                   >  1:7113710/1‑68 (MQ=255)
               aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTAGAAGCGGCTa                                                   >  1:13913308/1‑68 (MQ=11)
                cAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGCTAAGCGCAAAGACGCCGGCGACGACGCTAt                                                  >  1:26027318/1‑68 (MQ=11)
                 aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTAtt                                                 >  1:42563547/1‑68 (MQ=255)
                 aGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCGGCTAtt                                                 >  1:17884453/1‑68 (MQ=255)
                    aaCAGGCTAATGCGGAACCAGAAGCACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                               <  1:23928249/67‑1 (MQ=25)
                        ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                           <  1:25810879/67‑1 (MQ=35)
                         gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCgcg                                          >  1:22918663/1‑65 (MQ=25)
                         gCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCGGCTATTGCCCgcg                                          >  1:24511871/1‑65 (MQ=11)
                          cTAATTCGGAACCAGAAGAAGAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                           <  1:32136279/65‑1 (MQ=11)
                          cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTAGACGCCGCTATTGACCgacc                                         >  1:8518818/1‑64 (MQ=11)
                          cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTTGAAGCGGCTATTGCCCGAGc                                         >  1:42307341/1‑67 (MQ=11)
                          cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGCAGCCGCTATTGCCCGAGca                                        >  1:17017252/1‑67 (MQ=11)
                          cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGt                                           <  1:25575806/65‑1 (MQ=255)
                          cTAATGCGGAACCAGAAGAACAGGACGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGcc                                        >  1:22124295/1‑68 (MQ=11)
                               gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGTGCCAAAGc                                   >  1:36231612/1‑68 (MQ=25)
                               gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCACGTGCCAAAGc                                   >  1:26881354/1‑68 (MQ=25)
                               gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCGGCTATTGCCCGCGCCaagg                                    >  1:40257328/1‑65 (MQ=11)
                               gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCGaacg                                    >  1:37866618/1‑65 (MQ=11)
                               gCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGACGCCGAAGCAGCCGCTATTGCCCGGGCCAAAGc                                   >  1:34433031/1‑68 (MQ=11)
                               gCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCCGCAGTAGAAGCGGCTATTGCCCGCGCGaagg                                    >  1:1515787/1‑65 (MQ=11)
                                 ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCa                                  <  1:23465298/67‑1 (MQ=37)
                                 ggAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCa                                  <  1:15935749/67‑1 (MQ=37)
                                    aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCa                              >  1:26174563/1‑68 (MQ=40)
                                    aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAcgcacgca                              >  1:34894336/1‑68 (MQ=255)
                                     ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCCGCaa                             >  1:12525722/1‑68 (MQ=25)
                                     ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAGGCCCGCaa                             >  1:6381294/1‑68 (MQ=11)
                                      cacaagacCAGGTAGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:14686787/65‑1 (MQ=14)
                                      cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:43423475/68‑1 (MQ=40)
                                      cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:27055902/68‑1 (MQ=40)
                                      cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:25271479/68‑1 (MQ=40)
                                      cAGAAGAACAGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:39890188/68‑1 (MQ=255)
                                      cAGAAGAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAg                            <  1:18643522/68‑1 (MQ=255)
                                       agaagaACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                     <  1:32923824/42‑1 (MQ=38)
                                         aagaaCAGGTCGATCCGCGCAAAGCCCCCGTACAACCAGCAATAGCCCGTGCCAAAGCACGCAAGCTg                         >  1:21482411/1‑68 (MQ=14)
                                         aacaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                          <  1:21161560/64‑1 (MQ=1)
                                          acaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTgg                        <  1:37229923/66‑1 (MQ=14)
                                           gaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGa                       >  1:23688921/1‑68 (MQ=37)
                                            aaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                      <  1:5211172/68‑1 (MQ=35)
                                            aaCAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGa                       >  1:5028955/1‑67 (MQ=11)
                                            aaCAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGTGCCAAAGCACGAAAGCGGGa                       >  1:5793160/1‑67 (MQ=11)
                                             aCATGTAGATCCGCGCAAAGCCGCCGTAGAAGACGCTATTGACCGTGCCAAAGCACTCAACCTTGacc                     >  1:22845072/1‑66 (MQ=2)
                                             aCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAACGCACGCAAGCTGGAAc                     >  1:32748694/1‑68 (MQ=25)
                                             aCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAAc                     >  1:22389991/1‑68 (MQ=35)
                                             aCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTggca                      >  1:1331257/1‑65 (MQ=255)
                                             aCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAAc                     >  1:31541824/1‑68 (MQ=25)
                                               aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                    >  1:7485392/1‑67 (MQ=35)
                                                  tggATCCGCGCAAAGCCGCCGTGGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGc                  <  1:10826859/64‑1 (MQ=11)
                                                  tCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCa                 >  1:36777523/1‑67 (MQ=11)
                                                  tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCCACGCACGCAAGCGGGCCCAGCa                 >  1:27624246/1‑67 (MQ=11)
                                                  tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCGGGAACAGCa                 >  1:5531145/1‑67 (MQ=255)
                                                  tCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCa                 >  1:38384213/1‑67 (MQ=25)
                                                  gcGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGc                  <  1:11908712/65‑1 (MQ=25)
                                                  gcGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                 <  1:43656674/66‑1 (MQ=25)
                                                      tCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:29412392/68‑1 (MQ=255)
                                                      tCCTCGCAAACCAGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:39305589/68‑1 (MQ=11)
                                                      tCCGCGGAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:41383400/68‑1 (MQ=255)
                                                      tCCGCGCAAAGCGGCCGTCGAAGCCGCTTTTGCCCGTGCCAAAGCACGGAAGCTGGAACAGCAACAgg            <  1:10495849/68‑1 (MQ=11)
                                                      tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:32582510/68‑1 (MQ=35)
                                                      tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:28534530/68‑1 (MQ=35)
                                                      tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:1931057/68‑1 (MQ=35)
                                                      tCCGCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:16610479/68‑1 (MQ=255)
                                                      tCCCCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg            <  1:19693079/68‑1 (MQ=255)
                                                       ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAg             >  1:16265240/1‑66 (MQ=35)
                                                       ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAACGCACGCAAGCGGGAACAGCACCAggg           >  1:5393627/1‑67 (MQ=11)
                                                       ccGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGc           <  1:2317344/68‑1 (MQ=25)
                                                           gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAAt       >  1:31181665/1‑68 (MQ=35)
                                                             aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATg      <  1:6988036/67‑1 (MQ=35)
                                                             aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATg      <  1:15438435/67‑1 (MQ=35)
                                                                gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:35855107/1‑68 (MQ=40)
                                                                gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:20459946/1‑68 (MQ=40)
                                                                gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:41748799/1‑68 (MQ=40)
                                                                gccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa  >  1:37623455/1‑68 (MQ=40)
                                                                |                                                                   
AGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGA  >  CP009273/1702972‑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: