Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,709	T→A	13.5%	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
*	NC_000913	1,708,709	0	T	A	13.5%	64.7 / 3.1	52	R597R (CGT→CGA)	rsxC	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (5/36);  new base A (6/2);  total (11/41)
Fisher's exact test for biased strand distribution p-value = 7.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAC  >  NC_000913/1708561‑1708854                                                                                                                                                     |                                                                                                                                                  cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                                                                                                                  <  1:925711/150‑1 (MQ=255)   gcaACCGTTGCTGATCCACGTAAAACTTCCGTTGAAGCAGCTATCGCCCGCTCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                                                                                                <  1:1727466/150‑1 (MQ=255)   gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGAGcc                                                                                                                                                >  1:856559/1‑150 (MQ=255)   gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                                                                                                <  2:1339432/150‑1 (MQ=255)          ttgttgATCCACGCAAAGCAGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                                                                                                         <  1:570116/150‑1 (MQ=11)           tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                                                                        <  1:431707/150‑1 (MQ=255)            gCTGATCCACGTATAACTGCAGTTGAAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                                                                       <  2:283230/150‑1 (MQ=255)            gCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                                                                       <  2:38767/150‑1 (MQ=255)                 tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAGTGCGGAACCGGACCCGGCGGACGAACCGCGCCAGGCCGCCGTCGAAGCCGCTCTTGCCCGTGCCAAAGCACGCAGGcg                                                                                                                                  >  1:17076/1‑149 (MQ=17)                  ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                                 <  2:748071/150‑1 (MQ=255)                   cACTTAAAACTGCCGTTGAAGCAGACATCGCCCGCGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGGCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTgg                                                                                                                                <  2:1712250/150‑1 (MQ=17)                    cctttaaaCTGCCCTTCTACCAGCTAACGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                               <  1:611573/145‑1 (MQ=14)                    ccGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                               <  1:771066/149‑1 (MQ=255)                    aCGTAAAACTGCAGTTGAAGCAGCTATCTACCGCGCAAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAGCAGAACAAGAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                               <  2:685049/150‑1 (MQ=17)                      gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                             <  1:590321/150‑1 (MQ=255)                      ctAAAACTGCCGTTGAAGCAGCTATCGACCGCTCCAAAGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                                                                             <  2:822698/149‑1 (MQ=255)                         aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                          <  2:461868/150‑1 (MQ=255)                          taCTGCCGGAGTAGCAGCTATGTCCCGCGCCAACTCTCGCAATCTGGAACATCAACAGGAAACTCCGGAACCAGCACACCAGGCCGATCCGCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                                                                                                         <  2:1151136/149‑1 (MQ=19)                          aaCTGCCGTAGAAGCAGCTATTTCCCGCGCCTAAGCGCGCAAGCTGGAACACCAACAGGCTAATCAGGAACCAGAACCCCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                                                                                                                         <  2:40444/150‑1 (MQ=11)                                       gcagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                                                                            <  1:639889/150‑1 (MQ=255)                                           cTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                        <  1:1500528/150‑1 (MQ=255)                                           cTATCGCCCGCGCCAAAGCGCGGAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                        <  1:1814140/150‑1 (MQ=255)                                             aTCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCGCAGGCGCGCAAGCTGGAACAGCACCCGGCTAATGCGGAACCa                                                                                                      >  1:632185/1‑150 (MQ=14)                                                     cgcCAAATCGCGCAATCTGGAACAGCATCAGGCTAATGCGGAACCAGAACACCTGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAcg                                                                                              <  1:1799872/150‑2 (MQ=14)                                                     cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                                                                              <  1:1058046/150‑1 (MQ=255)                                                         aaatggcgcAATCTGGAAGATCTACAGGATAATGACGAACCAGAAAAACAGGTCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                          <  1:1037723/145‑1 (MQ=14)                                                         aaaGCTCGCAAGCTGGAACATCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                          <  1:1802796/150‑1 (MQ=255)                                                            tgccgcAAGCTGGATCATCAACAGGATAATTAGTAACCATTACAAATTTTAGATCCTCGTAAATCCTCCTTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACACCAACTGGCTAATGCGGAACCAGAACAACAGGTCGAt                                                                                       <  2:1161740/147‑1 (MQ=2)                                                            gctcgcAATCTGGAACATCAACAGGATAATTCGGAACCAGCACAACAGGTCGATCCGCGCAAATCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt                                                                                       <  1:911692/147‑1 (MQ=14)                                                             ggcggAATCTGGAACATCAACAGGAAAATTAGGAACCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                      <  1:388762/149‑1 (MQ=14)                                                             ggcgcAATCTGGAACATCATCAGGCAAATTCGGAACCAGAAACACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                                                                      <  1:395639/149‑1 (MQ=14)                                                                       ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                              <  1:562783/84‑1 (MQ=18)                                                                       ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                              >  2:562783/1‑84 (MQ=18)                                                                                   ggCTAATTCGGAACCAGAACCACAGGCAGATCCTCGAAAATCTGCCGTAGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc                                                                <  2:526279/150‑1 (MQ=14)                                                                                   ggCTAATTCGGAACCAGAACAAGAGGTCGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc                                                                <  2:1460507/150‑1 (MQ=255)                                                                                   ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaac                                                                                                                                   >  2:1719390/1‑81 (MQ=9)                                                                                           cGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACTCGAACAGCAACAGGCTAATGCCGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTTGAAGCGGCTATTg                                                        >  2:578448/1‑150 (MQ=11)                                                                                                  aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGc                                                 >  1:1057478/1‑150 (MQ=255)                                                                                                    aacaacAGGTCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                               <  2:627507/150‑1 (MQ=255)                                                                                                          ccgacGATCCTCGCAAATCCGCCGTCCAAGCAGCTATAGCCCGTCCCAAAGCGCGCAAGCTGGAACACCAACAGGCTAATGCGGAACCAGAACCACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCAc                                         <  1:680096/146‑1 (MQ=11)                                                                                                             tagATCCTCGCAAAACAGCCGTCGAAGCCGCTATTGCCCGGGCCAAATCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCa                                      <  2:1492609/148‑1 (MQ=14)                                                                                                             tCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCa                                      <  2:1734876/150‑1 (MQ=255)                                                                                                             ctgATCCTCGCAAACCCGCCGGCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAATCTGGAAGAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCa                                      <  2:1808616/148‑1 (MQ=11)                                                                                                                    gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAATCGGAACCAGAAGAACAGGTGGATCAGCGCAAAGAGGACGTCGAAGCGGGTATTGCCCGAGCCAAAACAAGAAAGATggg                               >  2:1383374/1‑149 (MQ=11)                                                                                                                     cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAATCGGAACCAGATCAACAGGTCGATCCGAGAAAAGCCGACGTCAAAGCCGCTATTGCCCGAGCAAAATCAAGAAAGCTGGaa                              >  2:10915/1‑150 (MQ=11)                                                                                                                       cAAAACGGCCGTGGCAGCCGCGATCCGCCGTGCCATATCGCGCAAGCTGGAACATCATCAGGATAATTCGGAACCAGAACACCATGCCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                            <  1:1136267/150‑1 (MQ=2)                                                                                                                       aaaacctgccgTCGAAGCAGCTAGTCCACGGGCCAAATCTCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                            <  1:1039195/143‑1 (MQ=1)                                                                                                                               ccgTCGACGCCGCTATTGCCCGGGCCACAGCGCGCAAGCTGGAAGAGCAACAGGCTACTGCGGAACCAGAACAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAgg                    <  1:490719/150‑1 (MQ=11)                                                                                                                                         cGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGCACAGCAACAGGCTAATGCGGACCCAGACGAACAGGTAGCTCCGCGCAAAGCCGCCGTCGCCGCGGCTATTGCCCGCGCGCACGCCCGCAAGCGGGCAGAGCCGCAGGCGaaaggggt          >  1:123241/1‑144 (MQ=2)                                                                                                                                                 cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGtt                                                                                          <  1:1167851/62‑2 (MQ=255)                                                                                                                                                 cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGtt                                                                                          >  2:1167851/1‑61 (MQ=255)                                                                                                                                                 cccGTGCAAAAGGTCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac  <  2:1722346/150‑1 (MQ=11)                                                                                                                                                     |                                                                                                                                                  CAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAC  >  NC_000913/1708561‑1708854

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: ‑