Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,709,163 T→C S8P (TCC→CCC)  rsxD → SoxR iron‑sulfur cluster reduction factor component, putative membrane protein of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,709,1630TC94.7% 125.0 / ‑0.2 38coding (22/1059 nt)rsxDSoxR iron‑sulfur cluster reduction factor component, putative membrane protein of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/19);  minor base . (1/1);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

GGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCG  >  NC_000913/1709027‑1709294
                                                                                                                                        |                                                                                                                                   
ggCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGcccc                                                                                                                                   >  1:1012968/1‑139 (MQ=255)
             ccAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATa                                                                                                                      >  2:93386/1‑139 (MQ=255)
               agaagaACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                    <  1:1654371/139‑1 (MQ=255)
               agaagaACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                    <  1:613551/139‑1 (MQ=255)
               agaagaACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                    <  1:1551721/139‑1 (MQ=255)
                          ttGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGAc                                                                                                         >  1:499356/1‑139 (MQ=255)
                            gATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACAt                                                                                                       >  2:8904/1‑139 (MQ=255)
                               ccGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCCAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                    >  2:837028/1‑123 (MQ=255)
                               ccGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAAc                                                                                                                    <  1:837028/123‑1 (MQ=255)
                                   gcAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCAt                                                                                                <  2:1326898/139‑1 (MQ=255)
                                    cAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATc                                                                                                      <  1:1356973/132‑1 (MQ=255)
                                    cAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATc                                                                                                      >  2:1356973/1‑132 (MQ=255)
                                                      ggCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTc                                                                             <  1:132454/139‑1 (MQ=255)
                                                          aTTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAg                                                                         <  1:93386/139‑1 (MQ=255)
                                                               ccGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTg                                                                    <  2:194093/139‑1 (MQ=255)
                                                                  cgcTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCa                                                                 >  1:1546190/1‑139 (MQ=255)
                                                                  cgcTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAG‑‑‑CCC‑‑‑‑‑CCCATAACCAGCGCCAGACATCGCGCATTATGCTg                                                                                         >  1:106941/1‑107 (MQ=255)
                                                                  cgcTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAG‑‑‑CCC‑‑‑‑‑CCCATAACCAGCGCCAGACATCGCGCATTATGCTg                                                                                         <  2:106941/107‑1 (MQ=255)
                                                                             aaaaaGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgc                                                      >  1:481782/1‑139 (MQ=255)
                                                                                   ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                              <  2:1528746/61‑1 (MQ=255)
                                                                                   ccgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTta                                                                                                                              >  1:1528746/1‑61 (MQ=255)
                                                                                    cgccCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACtg                                               <  1:948618/139‑1 (MQ=255)
                                                                                                  ttGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTttggttg                                 <  1:8904/139‑1 (MQ=255)
                                                                                                     tAAACGAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTgggg                              <  2:291773/139‑1 (MQ=255)
                                                                                                          gAGGACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTActc                         >  1:297/1‑139 (MQ=255)
                                                                                                            ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGcc                                                                  >  2:1433266/1‑96 (MQ=255)
                                                                                                            ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGcc                                                                  <  1:1433266/96‑1 (MQ=255)
                                                                                                            ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACtctc                       >  2:1111927/1‑139 (MQ=255)
                                                                                                            ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACtctc                       >  2:957567/1‑139 (MQ=255)
                                                                                                            ggACTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACtctc                       >  1:1435017/1‑139 (MQ=255)
                                                                                                               cTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCa                                                                                                 <  1:67518/62‑1 (MQ=255)
                                                                                                               cTAAATGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCa                                                                                                 >  2:67518/1‑62 (MQ=255)
                                                                                                                   aTGGTATTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGa                <  2:297/139‑1 (MQ=255)
                                                                                                                        aTTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATcgc                                                                                                    >  2:945040/1‑50 (MQ=255)
                                                                                                                        aTTCAGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATcgc                                                                                                    <  1:945040/50‑1 (MQ=255)
                                                                                                                            aGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCa                                                                                                 <  1:964563/49‑1 (MQ=255)
                                                                                                                            aGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCa                                                                                                 >  2:964563/1‑49 (MQ=255)
                                                                                                                            aGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcgc                                                    <  2:649253/94‑1 (MQ=255)
                                                                                                                            aGAATAGCTAGCCCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcgc                                                    >  1:649253/1‑94 (MQ=255)
                                                                                                                                 agctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGtttt                                        <  2:261484/101‑1 (MQ=255)
                                                                                                                                 agctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGtttt                                        >  1:261484/1‑101 (MQ=255)
                                                                                                                                 agctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCg  >  1:982957/1‑139 (MQ=255)
                                                                                                                                   ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                     >  1:990860/1‑86 (MQ=255)
                                                                                                                                   ctagccCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAgcg                                                     <  2:990860/86‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
GGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAACCAGCGCCAGACATCGCGCATTATGCTGTTGGTGTTGCTCGCAGCCGTGCCAGGAATCGCAGCGCAACTGTGGTTTTTTGGTTGGGGTACTCTCGTTCAGATCCTGTTGGCATCG  >  NC_000913/1709027‑1709294

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