Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,711,501	A→G	N225S (AAC→AGC)	rsxE →	SoxR iron‑sulfur cluster reduction factor component; electron transport inner membrane NADH‑quinone reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,711,501	0	A	G	100.0%	125.1 / NA	35	N225S (AAC→AGC)	rsxE	SoxR iron‑sulfur cluster reduction factor component; electron transport inner membrane NADH‑quinone reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/19);  total (16/19)

CACCGACTCCCCTTTCCTGCTGGCGATGCTGCCACCAGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAACGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCG  >  NC_000913/1711361‑1711648                                                                                                                                             |                                                                                                                                                    cacCGACTCCCCTTTCCTGCTGGCGATGCTGCCACCAGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACa                                                                                                                                           <  2:299808/151‑1 (MQ=255)      aCTCCCCTTTCCTGCTGGCGATGCTGCCACCAGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGaa                                                                                                                                      <  1:619479/151‑1 (MQ=255)                     tgGCGATGCTGCCACCAGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATaa                                                                                                                        >  2:583089/1‑150 (MQ=255)                      ggCGATGCTGCCACCAGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATaa                                                                                                                        <  1:722855/149‑1 (MQ=255)                                ccaccaGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCt                                                                                                            <  1:472470/151‑1 (MQ=255)                                ccaccaGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCt                                                                                                            <  1:552437/151‑1 (MQ=255)                                 caccaGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATaa                                                                                                                        >  1:452896/1‑138 (MQ=255)                                 caccaGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATaa                                                                                                                        <  2:452947/138‑1 (MQ=255)                                  accaGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTgg                                                                                                          >  2:620888/1‑151 (MQ=255)                                           aTTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCa                                                                                                 <  2:157568/151‑1 (MQ=255)                                             ttATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCAc                                                                                                <  2:459104/150‑1 (MQ=255)                                               aTTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCg                                                                                             <  2:767565/151‑1 (MQ=255)                                                      tGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGt                                                                                      >  2:265344/1‑151 (MQ=255)                                                           cTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTgaga                                                                                  >  1:701741/1‑150 (MQ=255)                                                                       ggAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATcc                                                                     <  2:620144/151‑1 (MQ=255)                                                                       ggAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATcc                                                                     >  2:703639/1‑151 (MQ=255)                                                                        gAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATccc                                                                    >  1:270265/1‑151 (MQ=255)                                                                        gAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATccc                                                                    >  1:642879/1‑151 (MQ=255)                                                                                 tgattgatGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCg                                                             >  1:182819/1‑149 (MQ=255)                                                                                 tgattgatGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCg                                                             >  1:310429/1‑149 (MQ=255)                                                                                   attgatGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCt                                                         <  2:575057/151‑1 (MQ=255)                                                                                     tgatgaAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCtt                                                        >  2:384289/1‑150 (MQ=255)                                                                                         gaAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAAttt                                                   >  1:418932/1‑151 (MQ=255)                                                                                         gaAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAAttt                                                   <  1:479829/151‑1 (MQ=255)                                                                                             gAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGt                                               >  1:575667/1‑151 (MQ=255)                                                                                                       gcgcCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                     <  2:477061/151‑1 (MQ=255)                                                                                                        cgcCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa                                    <  2:415319/151‑1 (MQ=255)                                                                                                         gcCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGaa                                   <  2:236030/151‑1 (MQ=255)                                                                                                            gTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTg                                >  1:692650/1‑151 (MQ=255)                                                                                                                   agcagcTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTg                         <  2:407759/151‑1 (MQ=255)                                                                                                                    gcagcTGCAGAACGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTg                         >  2:215210/1‑150 (MQ=255)                                                                                                                               aCGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCttt             <  2:568733/151‑1 (MQ=255)                                                                                                                                cGTGCATTGCCAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTc            <  2:537815/151‑1 (MQ=255)                                                                                                                                         ccAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGc   <  1:396166/151‑1 (MQ=255)                                                                                                                                          cAAGCGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCg  >  2:192171/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                                    CACCGACTCCCCTTTCCTGCTGGCGATGCTGCCACCAGGTGCATTTATTGGCCTGGGACTGATGCTGGCAGGAAAATACCTGATTGATGAAAGAATGAAAAAGCGCCGTGCTGAAGCAGCTGCAGAACGTGCATTGCCAAACGGTGAAACAGGGAATGTCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCG  >  NC_000913/1711361‑1711648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑