Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,994	A→C	19.3%	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	19.3%	76.3 / 2.2	43	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 (18/16);  new base C (6/3);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAC  >  NC_000913/1708855‑1709138                                                                                                                                            |                                                                                                                                                 gctctaccGATCTGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGCTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaa                                                                                                                                          <  7:192890/141‑1 (MQ=255) aacaGGTCGATCCGCGCAAAGCAGACGTCGAAGCCGCTATTGCCCTATCAAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaa                                                                                                                                          <  1:39484/148‑1 (MQ=17)      gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGACGCGGCTATTGCACGCGCGAAGGCAc                                                                                                                                     >  1:169225/1‑148 (MQ=255)            ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCAAAAGCGAGCAAACGGGAACAGAAACCGGCTAATGCGGATCCAGAAGAAGAGGTTGCTCCGCGCAAAGCTGCAGTAGAAGAGGCTATTGCCCGCGCGAAGgcgcgcgggc                                                                                                                               >  4:330685/1‑144 (MQ=14)            ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGCAGCGGCTATTGCCCGCGCAAACGCCCGCAAGc                                                                                                                               >  1:42411/1‑148 (MQ=14)               cgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGTGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGACGCGGCTATTGCACGCGCCATCGCACGCAAggggg                                                                                                                            >  8:240758/1‑144 (MQ=17)                                    gCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCt                                                                                                                              <  1:271103/125‑1 (MQ=255)                                    gCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCt                                                                                                                              >  2:271103/1‑125 (MQ=255)                                         tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGACGCGGCTATTGCCCGCGCAAAAGCCAGCAAGCTCGGAAGGGAACAGGGTAATGCCGtcg                                                                                                    >  8:196587/1‑144 (MQ=18)                                         tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAc                                                                                                    >  8:53525/1‑146 (MQ=255)                                           cccGAGCCAACGCGCGCAAACTTGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGAGGCTCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTaa                                                                                                            <  3:274901/136‑1 (MQ=14)                                           cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAACGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCag                                                                                                 >  1:9774/1‑147 (MQ=255)                                                 ccAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACa                                                                                          <  7:53525/148‑1 (MQ=255)                                                    aaGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAAc                                                                                                                   <  4:333117/120‑1 (MQ=255)                                                    aaGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAAc                                                                                                                   >  3:333117/1‑120 (MQ=255)                                                     aGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAACGCACGCAAGCGGGAACAACAACAGGCTAATGCGGGACCAGCAGAACAGGtg                                                                                      >  8:137429/1‑147 (MQ=37)                                                            aaaCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                                                                                                                                      >  5:147473/1‑93 (MQ=255)                                                            aaaCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                                                                                                                                      <  6:147473/93‑1 (MQ=255)                                                                  gAACAGAAACCGGCTTATGCGGAGCCAGATGAACAGTTTCATCAGCTCAAAGCTGCGGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCaa                                                                           <  7:256170/146‑1 (MQ=17)                                                                    acacCAACCGGCTAATGCCGAGCCAGAAGAACAGGTAGCTCCGCGCAAAGCGGCCGTCGAAGCGTCTATTGCACGCGCCAACGCTCGCATCCGGGATCAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTTGATCCGCGCAAAGCg                                                                       <  2:389461/148‑1 (MQ=11)                                                                         aaCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGt                                                                  >  3:16958/1‑148 (MQ=255)                                                                               cTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCAGGAACAGCAACAGGCTAATGCGGTACCAGAACAACAGGTTGAACCGCGCACAGCGGCAGTTGCcgc                                                            >  3:264318/1‑148 (MQ=255)                                                                                   tcaggAGCCAGAAGAACAGGTTTATCCGCGCCCCGCGGACGTCGCAGCGGCTCTTGCCCGCGCCAAAGCACGCAACCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCt                                                        <  4:301482/145‑1 (MQ=18)                                                                                    gCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTAATCACCTCGCCAAAGCACGCAAGCTGGACCAGCAACAAGCGACTGCGGTACCAGAAGAACAGGTTGAGCCGCGCAAAGCGGCa                                                                    >  7:314676/1‑135 (MQ=18)                                                                                         gCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTAGCAGCGGCTATTGCCCGCGCCAAGGCAAGAAAGCGGGCCCAGCAACAGGATAAAGCGGTGCCAGAAGAGGAGGTTGCACCGCGCAAAGCGGCCGTTGCCGCGGCTATTga                                                   >  2:294780/1‑146 (MQ=16)                                                                                           cagcagAACTGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcg                                                <  7:240758/148‑1 (MQ=255)                                                                                              aagaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGCTCCGCGCAAAGCGGCAGTTGCCGCGGCTCTTGCCCgcgc                                             >  1:287061/1‑148 (MQ=255)                                                                                              aagaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCTGCTATTGCCCgcgc                                             >  3:292384/1‑148 (MQ=255)                                                                                                  aCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCAATTGCCCGCGCTCAg                                         >  6:266822/1‑148 (MQ=255)                                                                                                       atGATCCGCGCAAACCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaa                                    <  7:196587/147‑1 (MQ=255)                                                                                                                     gccgCCGGCGAAGCGGCTATTGCCCGCGCCAAAGCACGCAATCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagc                      <  5:193538/145‑1 (MQ=255)                                                                                                                            tCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGt               <  1:69325/148‑1 (MQ=255)                                                                                                                             cGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTc                                           <  6:233417/119‑1 (MQ=255)                                                                                                                             cGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTc                                           >  5:233417/1‑119 (MQ=255)                                                                                                                               aaGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTg             >  1:112114/1‑147 (MQ=255)                                                                                                                               aaGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTGGATCCGCGCAAAGCCGCAGTGGACGCGGCTATTGCCCGCGCGAAGGCCAAAAAAGCCGCCCAGCAGCAgggt              >  2:69325/1‑144 (MQ=14)                                                                                                                                 gCGGCTATAGCACGCGCCACAGCACGAAAGCTGGAACCGCTTCAGGCCAATGCGGTTCCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTaa          <  1:387723/148‑1 (MQ=37)                                                                                                                                       aTTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGCCAGTTGCCGCGGCTATTGCCCGCGCTCCGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAg     >  1:217071/1‑147 (MQ=255)                                                                                                                                        ttGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg    <  7:162013/147‑1 (MQ=255)                                                                                                                                        ttGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg    <  2:112114/147‑1 (MQ=255)                                                                                                                                        ttGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg    <  3:95832/147‑1 (MQ=255)                                                                                                                                        ttGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGa   >  6:370754/1‑148 (MQ=255)                                                                                                                                          gCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAc  >  5:141004/1‑147 (MQ=255)                                                                                                                                            |                                                                                                                                                 AACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGAC  >  NC_000913/1708855‑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: ‑