Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,709,038	T→A	8.9%	V707E (GTA→GAA)	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,709,038	0	T	A	8.9%	129.0 / 2.0	50	V707E (GTA→GAA)	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 (24/20);  new base A (3/2);  total (28/22)
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

TGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAC  >  NC_000913/1708895‑1709180                                                                                                                                                |                                                                                                                                               tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGACGCGGCTATTGCCCGCGCAAAAGCCAGCAAGCTCGGAAGGGAACAGGGTAATGCCGtcg                                                                                                                                              >  8:196587/1‑144 (MQ=18) tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAc                                                                                                                                              >  8:53525/1‑146 (MQ=255)   cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAACGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCag                                                                                                                                           >  1:9774/1‑147 (MQ=255)         ccAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACa                                                                                                                                    <  7:53525/148‑1 (MQ=255)             aGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAACGCACGCAAGCGGGAACAACAACAGGCTAATGCGGGACCAGCAGAACAGGtg                                                                                                                                >  8:137429/1‑147 (MQ=37)                          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                                              <  2:112114/147‑1 (MQ=255)                                                                                                ttGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAgg                                              <  7:162013/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)                                                                                                     cgcgCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGCTCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTaaa                                        >  2:212873/1‑148 (MQ=255)                                                                                                         ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGt                                    <  5:370754/148‑1 (MQ=255)                                                                                                         ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGATGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATgg                                     >  3:298284/1‑147 (MQ=255)                                                                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTg                                                       <  4:323664/127‑1 (MQ=255)                                                                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTg                                                       >  3:323664/1‑127 (MQ=255)                                                                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTAt                                  >  5:252831/1‑148 (MQ=255)                                                                                                           aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTAt                                  >  3:300708/1‑148 (MQ=255)                                                                                                             aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAg                                                                                                             <  8:114687/71‑1 (MQ=255)                                                                                                             aGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAg                                                                                                             >  7:114687/1‑71 (MQ=255)                                                                                                                aggctaGCTGGAACTGCAACAGGCTAATGCGGTCCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGGTGTAAACGAGGACTAAATGGTATTCAGa                             <  1:98715/143‑1 (MQ=255)                                                                                                                  gCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAt                           <  6:252831/148‑1 (MQ=255)                                                                                                                     aGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAAGGGTATTCAGAAtagc                        >  3:236977/1‑148 (MQ=255)                                                                                                                             cAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCa                                                                                    >  3:133799/1‑80 (MQ=255)                                                                                                                             cAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCa                                                                                    <  4:133799/80‑1 (MQ=255)                                                                                                                              aGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagct                       <  6:282015/140‑1 (MQ=255)                                                                                                                              aGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagct                       >  5:282015/1‑140 (MQ=255)                                                                                                                                       cTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCa      >  5:117902/1‑148 (MQ=255)                                                                                                                                           tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAc  <  4:16958/148‑1 (MQ=255)                                                                                                                                            gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagctag                     <  7:172157/128‑1 (MQ=255)                                                                                                                                            gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAAtagctag                     >  8:172157/1‑128 (MQ=255)                                                                                                                                            gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAc  >  7:173030/1‑147 (MQ=255)                                                                                                                                                |                                                                                                                                               TGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAGAAGGTTGTAAACGAGGACTAAATGGTATTCAGAATAGCTAGCTCCCCTTATACCCATAAC  >  NC_000913/1708895‑1709180

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