Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,599,236	C→A	E175D (GAG→GAT)	dapA ←	4‑hydroxy‑tetrahydrodipicolinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,599,236	0	C	A	100.0%	150.0 / NA	43	E175D (GAG→GAT)	dapA	4‑hydroxy‑tetrahydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/20);  total (23/20)

TTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACA  >  NC_000913/2599104‑2599359                                                                                                                                     |                                                                                                                            tttGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGa                                                                                                                      >  1:621737/1‑140 (MQ=255)           tGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTac                                                                                                            <  2:621737/140‑1 (MQ=255)             ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTacac                                                                                                          >  2:6525/1‑140 (MQ=255)              gCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACAcg                                                                                                         <  2:610403/140‑1 (MQ=255)                  tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACAcg                                                                                                         <  2:671416/136‑1 (MQ=255)                  tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACAcg                                                                                                         >  1:671416/1‑136 (MQ=255)                  tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGtt                                                                                                     <  2:573317/140‑1 (MQ=255)                  tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGtt                                                                                                     >  2:220678/1‑140 (MQ=255)                  tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGtt                                                                                                     >  2:19450/1‑140 (MQ=255)                  tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGtt                                                                                                     <  1:758270/140‑1 (MQ=255)                         cgcTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTcc                                                                                              >  2:713902/1‑140 (MQ=255)                            tGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTg                                                                                           >  2:257912/1‑140 (MQ=255)                              cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTGCCTGtt                                                                                         >  1:788043/1‑140 (MQ=255)                              cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGtt                                                                                         <  2:622379/140‑1 (MQ=255)                                aCGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGc                                                                                       <  1:6525/140‑1 (MQ=255)                                aCGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGc                                                                                       <  1:713902/140‑1 (MQ=255)                                         gTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGAt                                                                              <  1:534581/140‑1 (MQ=255)                                                               cGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTc                                                        >  2:485084/1‑140 (MQ=255)                                                                 cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                      >  1:409602/1‑140 (MQ=255)                                                                 cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                      >  2:458118/1‑140 (MQ=255)                                                                 cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                      >  2:345182/1‑140 (MQ=255)                                                                 cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                      >  2:124206/1‑140 (MQ=255)                                                                 cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                      >  1:540798/1‑140 (MQ=255)                                                                          aTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGcc                                             <  1:485084/140‑1 (MQ=255)                                                                              aGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAcc                                         <  1:220678/140‑1 (MQ=255)                                                                                tCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGcc                                             >  2:688879/1‑134 (MQ=255)                                                                                tCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGcc                                             <  1:688879/134‑1 (MQ=255)                                                                                tCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGt                                       >  1:852748/1‑140 (MQ=255)                                                                                   aGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGc                                              >  1:417911/1‑130 (MQ=255)                                                                                   aGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGc                                              <  2:417911/130‑1 (MQ=255)                                                                                       cgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                <  1:458118/140‑1 (MQ=255)                                                                                       cgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                >  2:570566/1‑140 (MQ=255)                                                                                          tCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGc                             <  1:124206/140‑1 (MQ=255)                                                                                          tCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGc                             >  1:608534/1‑140 (MQ=255)                                                                                                    gccgcTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGc                   <  2:788043/140‑1 (MQ=255)                                                                                                    gccgcTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGc                   >  1:125395/1‑140 (MQ=255)                                                                                                    gccgcTCAGCAGAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGc                   <  2:540798/140‑1 (MQ=255)                                                                                                             cagAACAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACggg          >  2:209706/1‑140 (MQ=255)                                                                                                                aaCAAAATCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACg       <  2:417472/140‑1 (MQ=255)                                                                                                                     aaTCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcaca  <  2:125395/140‑1 (MQ=255)                                                                                                                     aaTCATCTGAAACCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcaca  >  1:522527/1‑140 (MQ=255)                                                                                                                                aCCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGa         <  2:587170/122‑1 (MQ=255)                                                                                                                                aCCAGATCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGa         >  1:587170/1‑122 (MQ=255)                                                                                                                                     |                                                                                                                            TTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACA  >  NC_000913/2599104‑2599359

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