Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,994	A→C	21.8%	A692A (GCA→GCC)	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,994	0	A	C	21.8%	65.1 / 3.7	47	A692A (GCA→GCC)	rsxC	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (26/9);  new base C (9/2);  total (35/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAA  >  NC_000913/1708848‑1709141                                                                                                                                                   |                                                                                                                                                    ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCAcgc                                                                                                                                                  >  1:1523337/1‑150 (MQ=255)  cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCAAAAGCGCGCAAACGGGAACAGCAACAGGCTATTGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGTCAGAAGAAGCGCCTAATGCAcgcc                                                                                                                                                 >  1:1258447/1‑149 (MQ=14)        cccaGGCCGATCCCCGCAAATCTGCCGTCGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAACAACAGGTTGATCCTCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAg                                                                                                                                           <  1:820570/148‑1 (MQ=255)         cccgggcGATCCGCGCAAAGCCCCCGCCCACGCCGCTATTGCCCGAGCCAAAGCGCGCAAACTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGc                                                                                                                                          <  1:920990/144‑1 (MQ=21)            ggTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGAGCCAACGCGCGCAAACGGGAACAGCATCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCAAGGCCCg                                                                                                                                       >  1:1542706/1‑150 (MQ=17)            ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCACCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCGATTGCCCGCGCCAACGCACg                                                                                                                                       >  1:1551220/1‑150 (MQ=14)            ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGACGATCAGGTTGCTCCGCGCATCGCGGCCGTGGCAGAGGCTGTTGCCCGCGCCaccgccag                                                                                                                                       >  1:394430/1‑143 (MQ=14)                   ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGTGCCAAAGCACGCAAGCTg                                                                                                                                >  2:1390264/1‑150 (MQ=255)                    cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAGCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCCAAGCGGCCGTGGCCGCGGCGATTGCCCGCGCACAGGCCCGCAAggagg                                                                                                                               >  1:246041/1‑146 (MQ=14)                    cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAACGCACGCAAGCTgg                                                                                                                               >  1:1452173/1‑150 (MQ=255)                            gccgccGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCGCACGCCAGCAAGCTCGCCCAGCaa                                                                                                                       >  2:109592/1‑150 (MQ=17)                             ccgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCGAAGCTGCCGGCGAAGCGGGTATTGCCCGAGCCGAAGCACGCACGCTGGACCAGCAAc                                                                                                                      >  1:1476164/1‑150 (MQ=14)                              cgccgTTTAAGCTTCTATTCCCCGCTCCAAAGCGAGCAATCGGGAACAGCAACCTTCTAATGCGGTTCCCGATTAAGAGGTTGATCCTCGCAAAGCTGCCGTCGAAGCTGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACa                                                                                                                     <  2:628827/150‑1 (MQ=14)                                      aaGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCTGGAACAGCAGCAGGCTAATg                                                                                                             >  1:1060267/1‑150 (MQ=255)                                      aaGCCGCTATTGCCCGAGACAAATCTCGCAATCGGGAACAGCAACAGGCTAATGCGTATCCAGAATAACAGTTCGATCCGCGCAAAGGTGCCGTCGAAGCCGCTATTGCACGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATg                                                                                                             <  2:957363/150‑1 (MQ=11)                                         ccGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAACGCTGCCGTCGCCGCGGCGATTGCACGCGCCACCGCCCGCCAGCGGGCGCAGCAACAGGCTACTGCgg                                                                                                          >  2:1246263/1‑150 (MQ=21)                                                tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCACCAGGCTAATGCGGTACCaga                                                                                                   >  1:1733196/1‑150 (MQ=255)                                                 gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaa                                                                                                  >  2:1561707/1‑150 (MQ=255)                                                   ccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaaga                                                                                                >  1:1530642/1‑150 (MQ=255)                                                          aaaGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCTCGAACAGCAACAGGCTAATCAGGAACCATAAAAACTGGtt                                                                                         >  2:1368985/1‑150 (MQ=25)                                                            aGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCGGGACCAGCAGCAGGCTAATGCGGTACCAGAAGAACAGGTGGa                                                                                       >  2:1890962/1‑150 (MQ=255)                                                                    aaCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCCCGCGCCAACGCACGCAAGCACGAACAGCAACAGGCTAATGCCGTACCAGAATAACAGGTTGATCagcgaa                                                                               >  2:969375/1‑148 (MQ=17)                                                                      cGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCaaa                                                                             >  1:5565/1‑150 (MQ=255)                                                                        ggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAACGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTCCCCGAAGAACAGGTTGATCCGCGCACAGc                                                                           >  1:955465/1‑150 (MQ=255)                                                                         gAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCGGGAACAGCACCAGGCTAATGCGGTCCCAGAAGAACAGGTTGATCCGCGCAAAGCg                                                                          >  2:1233137/1‑150 (MQ=255)                                                                            cAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCa                                                                       >  1:482395/1‑150 (MQ=255)                                                                             aGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCCGCGGCTATTGCCCGCGCCCAAGCAATCAACCACGAACAGCAACTGGCTTTTGACGAGCCAGAAGACGAGGTATATACGCGAAGCGCGGCAg                                                                      >  2:698618/1‑150 (MQ=9)                                                                                       tAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAGCAGGCTAATGCGGTACCATAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                                                            >  2:1450472/1‑150 (MQ=255)                                                                                                        caCAGGTCGCTCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCGGGAACAGCAACAGGCTAACGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg                                           <  1:1014571/149‑1 (MQ=18)                                                                                                         aCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGc                                          <  1:1390265/150‑1 (MQ=255)                                                                                                         aCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGc                                          >  2:245327/1‑150 (MQ=255)                                                                                                             gTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaa                                      >  2:613857/1‑150 (MQ=255)                                                                                                                gATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaaa                                   >  1:1116882/1‑150 (MQ=255)                                                                                                                 ctCCGCCCAAAGCGGCCTTCGAAGCGTATAATGCAGTCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAg                                  <  2:1502334/149‑1 (MQ=21)                                                                                                                 aTCCGCGCAAAGCTGCCGTCGACGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCATCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCACAGCGGCAGTTGCCGCGGCTATGGCCCGCGCTCAGGCcaacaacg                                  >  1:344665/1‑148 (MQ=25)                                                                                                                  tCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGACCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAACAAgc                                 >  1:878875/1‑150 (MQ=255)                                                                                                                  tCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgc                                 >  2:933715/1‑150 (MQ=255)                                                                                                                   ccGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc                                >  2:178165/1‑150 (MQ=255)                                                                                                                   ccGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc                                >  2:610071/1‑150 (MQ=255)                                                                                                                    cgcgCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccg                               >  1:467434/1‑150 (MQ=255)                                                                                                                     ccgcAATGCTGCCGTCGTCGCTGCTTTTGCTCGCGCATAAGCGCGCAAGGGGGATCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccgc                              <  1:1873981/149‑1 (MQ=18)                                                                                                                                   tCGAAGCGGCTATTGCACGATCCAAATCACGCATGTTTTAAGATCAACAGGCTAATGCGGTTCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTg                <  2:700750/150‑1 (MQ=25)                                                                                                                                    cGACGCGGCTATTGACCGCGCAAAAGCACGCAAGCGGCAACAGCAACCTGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGt               <  2:998589/150‑1 (MQ=21)                                                                                                                                              tttGCACGGTCTAAAACACGCATGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAc     <  1:909528/149‑1 (MQ=37)                                                                                                                                                tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt                  <  1:1104081/135‑1 (MQ=255)                                                                                                                                                tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt                  >  2:1104080/1‑135 (MQ=255)                                                                                                                                                 gCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTaa  >  2:899827/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAA  >  NC_000913/1708848‑1709141

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