Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	28,427	C→T	G18G (GGC→GGT)	dapB →	4‑hydroxy‑tetrahydrodipicolinate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	28,427	0	C	T	100.0%	118.7 / NA	36	G18G (GGC→GGT)	dapB	4‑hydroxy‑tetrahydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/22);  total (14/22)

GCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAGGCGTTACCGTGC  >  NC_000913/28306‑28560                                                                                                                          |                                                                                                                                      gctggcGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAggcg                                                                                                                     <  2:7297/140‑1 (MQ=255)    ggcggcGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGgcggcg                                                                                                                  <  2:68789/140‑1 (MQ=255)     gcggcgTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGgcgc                                                                                                                 <  2:114271/140‑1 (MQ=255)     gcggcgTAGCGATGCGCAGGTTTCTCTGAAATCGGACTATGCAAATTATCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGgcgc                                                                                                                 <  2:308883/140‑1 (MQ=255)              cGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTag                                                                                                        <  1:580448/140‑1 (MQ=255)               gATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGgcgc                                                                                                                 >  2:350354/1‑130 (MQ=255)               gATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGgcgc                                                                                                                 <  1:350354/130‑1 (MQ=255)               gATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTaga                                                                                                       >  2:126737/1‑140 (MQ=255)                    gcTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTTATTCAGGCGGCGCTGGCATTAGAGGGCg                                                                                                  <  2:69237/140‑1 (MQ=255)                              ctGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTggg                                                                                        <  2:614676/140‑1 (MQ=255)                                gAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcg                                                                                      <  2:638513/140‑1 (MQ=255)                                  aaaCGGTCTATGCAACTTAACACAACAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgct                                                                                    <  1:317768/140‑1 (MQ=255)                                   aaCGGTCTATGCAAATTAACAACAGTGAATAGCTATGAATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGgcgctg                                                                                   <  1:246817/140‑1 (MQ=255)                                                aaTTAACACAAGAGAATGGCTAGGCATGATTCACAGATCCGCGTTGCCTTCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGa                                                                      <  2:326236/140‑1 (MQ=255)                                                 aTTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGaa                                                                     <  2:115107/140‑1 (MQ=255)                                                 aTTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATCAGAGGGCCTGCCGTCGGGCGCCGCGCTGGAGCGTGaa                                                                     >  1:326236/1‑140 (MQ=255)                                                   tAACAAAAGAGAATAGCTATGCCTGATGCAAACGTCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAgg                                                                   <  2:289366/140‑1 (MQ=255)                                                   tAACAAAAGAGAATAGCCATGCATGATGCAAACATTCGCATTGCCATCCCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAgg                                                                   <  2:225703/140‑1 (MQ=255)                                                    aaCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGt                                                                        <  2:385644/134‑1 (MQ=255)                                                    aaCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGt                                                                        >  1:385644/1‑134 (MQ=255)                                                    aaCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGa                                                                  <  1:360225/140‑1 (MQ=255)                                                     aCAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAt                                                                 <  1:421014/140‑1 (MQ=255)                                                      cAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                                >  2:492180/1‑140 (MQ=255)                                                      cAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtc                                                                >  1:87902/1‑140 (MQ=255)                                                        aaaGAGAATAGCTATGCTTGATGCAAACATCCGCGTTGCTATCGCGGGCGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGAtctt                                                              <  2:226923/140‑1 (MQ=255)                                                                 aGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGAGGCGCTGGCATTAGAGTGCGTACAGTTGGACGCGGCGCTGGAGCGTGTAGGATCTTCTTTACTgg                                                     >  2:174086/1‑140 (MQ=255)                                                                      tGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTT‑GGCGCGTGGCTGGAGCGAGAAAGGTCTTCTTTACTGGGCAACg                                               >  1:215511/1‑140 (MQ=255)                                                                           gatgCAAACATCCGCGTTGCCATCGCGGGAGCTGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGc                                           >  2:381272/1‑140 (MQ=255)                                                                                         cgTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCg                             >  2:643056/1‑140 (MQ=255)                                                                                         cgTTGCCATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTGACTGGGAAGCGACGCCGGTGAGCTGGCCg                             >  1:184778/1‑140 (MQ=255)                                                                                              ccATCGCGGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGcc                        >  2:527844/1‑140 (MQ=255)                                                                                                    cgGGAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGaaa                  <  1:643056/140‑1 (MQ=255)                                                                                                      ggAGCCGGGGGGCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAAc                <  2:87902/140‑1 (MQ=255)                                                                                                             gggggCGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAGGCGtt         >  1:612994/1‑140 (MQ=255)                                                                                                                  cGTATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAGGCGTTACCGt    >  1:483458/1‑140 (MQ=255)                                                                                                                    tATGGGTCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAGGCGTTACCGTGc  <  2:429387/140‑1 (MQ=255)                                                                                                                          |                                                                                                                                      GCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACATCCGCGTTGCCATCGCGGGAGCCGGGGGGCGTATGGGCCGCCAGTTGATTCAGGCGGCGCTGGCATTAGAGGGCGTGCAGTTGGGCGCTGCGCTGGAGCGTGAAGGATCTTCTTTACTGGGCAGCGACGCCGGTGAGCTGGCCGGAGCCGGGAAAACAGGCGTTACCGTGC  >  NC_000913/28306‑28560

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