Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,662 C→G 11.4% Q582E (CAA→GAA)  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
*NC_0009131,708,6620CG11.4% 103.7 / 2.6 48Q582E (CAA→GAA) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (7/35);  new base G (3/3);  total (10/38)
Fisher's exact test for biased strand distribution p-value = 9.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCT  >  NC_000913/1708519‑1708796
                                                                                                                                               |                                                                                                                                      
aaGCCAGAGCGAAACAGGCAGAACTGCAGCAGACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAgg                                                                                                                                  <  1:1015637/150‑1 (MQ=255)
      gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                                            <  2:904768/150‑1 (MQ=255)
      gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                                            <  2:891373/150‑1 (MQ=255)
      gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                                            <  2:223314/150‑1 (MQ=255)
                   aGAACTGCAGCAAACTAACGACGCAGCAACGGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgccgc                                                                                                               <  1:881563/150‑1 (MQ=255)
                        tGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                                                                          <  2:897312/150‑1 (MQ=255)
                        tGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGCCAAAGGAAGAACCAGGACTACAGGACGATCCGAGAAACGACGCCGTCg                                                                                                          >  2:69113/1‑150 (MQ=25)
                         gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGa                                                                                                         <  1:1103721/150‑1 (MQ=255)
                         gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGa                                                                                                         <  1:908913/150‑1 (MQ=255)
                          aagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                                                        <  2:1192998/149‑1 (MQ=255)
                            gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTAGATCAGAGCGCATCCGCCCTCgtcgc                                                                                                      >  1:900690/1‑146 (MQ=21)
                            gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGAGACAGGCTAATGCGGCCCCAGAACACCAGGGCGCTCCGCGCCACGCCGCCGTCGAAGc                                                                                                      >  1:955900/1‑150 (MQ=21)
                             tcatctcacgaagCAGCACCCCTTGCTGATCCTCGTAAATCTGCCGTCGAAGCAGCTATAGCCCGATCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                     <  2:504443/143‑1 (MQ=11)
                             caAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                     <  1:813478/150‑1 (MQ=255)
                                 cTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCggat                                                                                                 >  1:504443/1‑147 (MQ=255)
                                  tAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                                                                <  2:920/150‑1 (MQ=255)
                                   aaCGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACACGCTAATGCGGAGCCAGAAGAAAAGGTCGATACGCGAAAAGACGCCGTCGAAGCCGCTAtt                                                                                               >  2:1187564/1‑150 (MQ=18)
                                      gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                                            <  1:831308/150‑1 (MQ=255)
                                            gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                                      <  1:7979/150‑1 (MQ=255)
                                              aaCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATGCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                                                    <  2:1128807/150‑1 (MQ=255)
                                                ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCAAAAGCGCGCAAGCTTTAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAg                                                                                  <  2:400759/150‑1 (MQ=255)
                                                    tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                              <  1:476303/150‑1 (MQ=255)
                                                    cgCTGATCCACGTAAACCTGCCGTACAAGCAGCTAGCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                              <  2:893340/149‑1 (MQ=17)
                                                           ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGGCAAAGCACGCaacagg                                                                       >  1:1123090/1‑146 (MQ=255)
                                                               gTATAACTGCAGATGTTGCAGCGATCGCCCGGGCCTAAGCGGGCAAGCTGGAACAGCTTCAGGATAATTCGGAACCAGAACACCAGGCCGATCCGCGTAAACCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                   <  2:127755/150‑1 (MQ=14)
                                                                  aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCATCAGGCTAATGCGGAACCAGAACAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                <  2:520692/150‑1 (MQ=255)
                                                                  aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTg                                                                                              >  1:1353498/1‑120 (MQ=14)
                                                                  aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTg                                                                                              <  2:1353498/120‑1 (MQ=14)
                                                                     cTGCCGTTGAAGCAGCTCGCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAAc                                                             <  1:1212662/150‑1 (MQ=255)
                                                                     cTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg                                                                                                             <  1:197121/102‑1 (MQ=255)
                                                                     cTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg                                                                                                             >  2:197121/1‑102 (MQ=255)
                                                                     cTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAAc                                                             <  1:1218266/150‑1 (MQ=255)
                                                                       gCCGTTGAAGCAGCTCTGGCCCGCGACCACGCGCGCAAGCTGGACCTGCAACAGGCTACTCCGGACCCAGAACACCAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg                                                           <  1:69113/150‑1 (MQ=14)
                                                                         cGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGccgccg                                                                                                             <  2:1328024/98‑1 (MQ=255)
                                                                         cGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGccgccg                                                                                                             >  1:1328024/1‑98 (MQ=255)
                                                                          ggtGAAGCAGCTATCGCCCGCGCCAACTCGCGCAATCTGGAACAGCAACAGGATAATGCGGAACCAGAACCCCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCt                                                        <  2:811202/148‑1 (MQ=14)
                                                                           ttGCAGCAGCTATCGCCGGCGACAAATCGCGCAAGCTGGCACACCATCAGGCTACTGCGGAACCAGAACAACAGGTCGATCCCCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                                                       <  2:783413/150‑1 (MQ=14)
                                                                                        cGCCCGCGCAAAATCGCGCAAGCTGGAACATCATCAGGCTAATGCGGAACCAGCACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga                                          <  1:955610/150‑1 (MQ=14)
                                                                                          cccGGGCCAAAGCGCGCAAGCTGGAACATCTACAGGCTACTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                                        <  1:1062036/150‑1 (MQ=14)
                                                                                              cgcCAAATCGCGGAAGCTGGAACAGCAACAGGCTAATTCGGAACCAGAACAACAGGTAGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                    <  2:553872/150‑1 (MQ=14)
                                                                                              cgcCAAAGCGCGCAAGCTGGAACATCAACAGGATAATGCGGAAGCAGAACAAGAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                    <  2:1185258/150‑1 (MQ=14)
                                                                                                  aaaGCTCGCAATCTGGAACATCATCAGCATAATTAGGAACCAGAACAAGAGGTCGATCCTCGCAAAGCAGCCGTCGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                <  1:1281030/150‑1 (MQ=14)
                                                                                                          cAAGCAGGAACAGCATCAGGATAATCCGGAACCAGAACAACAGGTCGATCCGCGTAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg                        <  2:735197/150‑1 (MQ=14)
                                                                                                                      tcAACAGGCTAATGCGGAAGCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGt            <  1:626289/149‑1 (MQ=255)
                                                                                                                        aaCAGGATAATTCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCg          <  2:570275/150‑1 (MQ=255)
                                                                                                                           aGGCTAATGCCGAACCAGAACACCAGGACGATCCGCGAAAAGCCGCCGGCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg       <  2:496749/150‑1 (MQ=14)
                                                                                                                            ggCTAATGCGGAACCAGAACAACAGGTCGCTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGAATAATACGGAACCACAAGAACCGGTCGAACATCGCAAACACGCCGTCGACGc      >  1:608470/1‑150 (MQ=14)
                                                                                                                                aaTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt  <  1:620078/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                      
AAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCT  >  NC_000913/1708519‑1708796

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