Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,994	A→C	15.5%	A692A (GCA→GCC)	rsxC →	SoxR [2Fe‑2S] reducing system protein RsxC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,708,994	0	A	C	15.5%	69.5 / 3.8	43	A692A (GCA→GCC)	rsxC	SoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base A (23/13);  new base C (4/3);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAC  >  NC_000913/1708847‑1709138                                                                                                                                                    |                                                                                                                                                 gCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAAAGGGAACAGCAACCGGCTAATGAGGAGCCAGAAGAACAGGTTGATACGCGCAAAGCTGCCGTGGCCGCGGCTGTTGcgcg                                                                                                                                                >  1:306375/1‑150 (MQ=17)       acaactGGTCCATCCCCCCCAAGCCCCCGGCGAAGGCTCTCTGCCCCGAGCCAACGCGCGCCAACTGCACCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTGGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaa                                                                                                                                          <  1:251959/150‑1 (MQ=255)             ggTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGAAAACGGGAACAGCAACAGGCTAATGAGGAGCAAGAAGAACAGGATGATCAGCGCAAAGCTGCAGTACACGCGGCTATTGCACGCGCAAAAGCaaa                                                                                                                                    >  2:609504/1‑148 (MQ=17)                        gcAATGCCGCCGTCGAAGCCGGTATTGCCGTAGCCTTAGCTCGCTATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAAc                                                                                                                         <  1:784407/150‑1 (MQ=17)                              ccgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCGAAAGCCCGAAAGCAGGACCAGCAAc                                                                                                                   >  1:1080647/1‑150 (MQ=17)                                      gAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCACGCACAGCAACAGGATaaa                                                                                                           >  2:1165463/1‑149 (MQ=25)                                              tATTGCCCGAGCCAAAGCGCGCAAACGGGAACATCAACATTCTAATGCGGAGCCAGAATAACATTTTGATCCGCTCAAAGCTTCCGTCGAAGCGGCTATTGCAAGCTCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAcc                                                                                                   <  2:756499/150‑1 (MQ=17)                                              tATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAcc                                                                                                   >  2:1403161/1‑150 (MQ=255)                                                  gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaa                                                                                               >  1:935452/1‑150 (MQ=255)                                                  gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAACCGGGAACAGCACCAGGCTAATGCGGAGCCagaa                                                                                               >  1:799190/1‑150 (MQ=37)                                                  gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaa                                                                                               >  2:1189345/1‑150 (MQ=255)                                                  gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCCCGCGCAAACGCACGCAACCGCGCCCAGCAACAGGATCAAGCGGtcgccgaa                                                                                               >  2:1004418/1‑143 (MQ=18)                                                  gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCCAAAGCACGAAAGCAGGAACAGCAACAGGCTAATGCGGTACCagaa                                                                                               >  1:693477/1‑150 (MQ=255)                                                   cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTCGAACAGCAGCAGGCTAATGCGGTACCagaag                                                                                              >  1:272508/1‑150 (MQ=255)                                                           aaaGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCCGCGGCTATTGCCCGCGCCCAAGCACGCAAGCAGGAACAGCAACAGGCTAATGCGGAACCCGAAGAACACGta                                                                                      >  1:1229572/1‑149 (MQ=17)                                                                 cgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCGGGCACAGCAACAGGCTAATGCGGGACCAGAAGAACAGGTTGATcag                                                                                >  2:140042/1‑148 (MQ=18)                                                                                   ccGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGcc                                                              >  1:990424/1‑150 (MQ=255)                                                                                    cGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTATTGCCCGCGCCAAAGCAAGCAAGCTGGAACAGCACCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCAGCGCAAAGCGGCAGTTGccc                                                             >  2:1299593/1‑149 (MQ=38)                                                                                     ggCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCAAGCAAGCTGGAACAGCAAAAGGCTAATGAGGTACCAGAAGAACAGGTTGATCCCCGCAAAGCGGCAGTTGCcgc                                                            >  1:617033/1‑150 (MQ=255)                                                                                            gCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAtt                                                     >  1:390185/1‑150 (MQ=255)                                                                                            gCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAtt                                                     >  1:1060157/1‑150 (MQ=255)                                                                                              ggAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGc                                                   <  1:548580/150‑1 (MQ=255)                                                                                                   cAGAAGAACAGGTTTATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCGAAAGCAAGCAAGCGGGAACAGCAACAGGCTAATGCGGTACAAGAAGAACAGGTTCGTCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgcg                                              >  2:706890/1‑150 (MQ=21)                                                                                                   cAGAACAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCGGGAACAGCAACAGGCTAATGCGGTCCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgcg                                              <  1:765585/150‑1 (MQ=14)                                                                                                         aaCAGGTTGATCCTCGCAAAGCCGCCTTCGAAGCTTCTATTGCCCGTGCCATCGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg                                        <  2:388053/150‑1 (MQ=17)                                                                                                             ggTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCCAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaa                                    >  1:812976/1‑150 (MQ=255)                                                                                                             ggTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaa                                    >  2:818170/1‑150 (MQ=255)                                                                                                                tGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaa                                 >  1:590573/1‑150 (MQ=255)                                                                                                                 gATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaaa                                >  2:1150023/1‑150 (MQ=255)                                                                                                                  aTCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAg                               >  1:1427266/1‑150 (MQ=255)                                                                                                                               ggCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGaa                  <  1:609504/149‑1 (MQ=255)                                                                                                                               agCCGTCGCATCGGCGACTGCAGGCGCCAACGCTCGCAATCCGCAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGaa                  <  1:748369/149‑1 (MQ=25)                                                                                                                                  cGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt               >  1:674099/1‑150 (MQ=255)                                                                                                                                   gTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCTAAAAAAGCCGCCCAGCAGAAGGtt              >  2:9461/1‑150 (MQ=255)                                                                                                                                    tCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTg             <  1:1189345/150‑1 (MQ=255)                                                                                                                                      gAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTa           <  1:140042/150‑1 (MQ=255)                                                                                                                                        aGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTaaa         <  2:1080647/150‑1 (MQ=255)                                                                                                                                           ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACgc      >  1:615662/1‑149 (MQ=255)                                                                                                                                            gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGa                                                                                    <  1:356600/71‑1 (MQ=255)                                                                                                                                            gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGa                                                                                    >  2:356600/1‑71 (MQ=255)                                                                                                                                            gCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAg     <  2:14336/150‑1 (MQ=255)                                                                                                                                             cTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg    <  2:935452/150‑1 (MQ=255)                                                                                                                                              tATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCTAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGa   <  1:9461/150‑1 (MQ=255)                                                                                                                                               aTTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAc  <  2:1172478/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAC  >  NC_000913/1708847‑1709138

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