Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,402,336	Δ1 bp	coding (920/1446 nt)	abgB ←	p‑aminobenzoyl‑glutamate hydrolase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,402,333	0	T	.	100.0%	85.1 / NA	21	coding (923/1446 nt)	abgB	p‑aminobenzoyl‑glutamate hydrolase subunit B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (11/10);  total (11/10)

AATTTGTTTCGCAAAAGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAA  >  NC_000913/1402257‑1402409                                                                             |                                                                             aaTTTGTTTCGCAAAAGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCAtttc                                                                           >  1:268680/1‑79 (MQ=255)   tttGTTTCGCAAAAGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCAtttct                                                                          >  1:268683/1‑77 (MQ=255)         tCGCAAAAGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGc                                                                  <  1:268658/80‑1 (MQ=255)             aaaaGCCAGTTCTTCGGAGTCCCATTCCGGGGTACCAAAATGGGTTAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTc                                                             >  2:268690/1‑81 (MQ=255)               aaGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCgg                                                           <  1:268659/81‑1 (MQ=255)               aaGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCgg                                                           >  2:268697/1‑81 (MQ=255)                       tcttcGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAAc                                                   >  2:268699/1‑81 (MQ=255)                               gTTCCATTCCGGGGTCTCAAAATGGGTTAGGTCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAAGTGGAACAg                                           >  2:268691/1‑81 (MQ=255)                                ttCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAgg                                          >  2:268687/1‑81 (MQ=255)                                  ccATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCt                                        >  2:268693/1‑81 (MQ=255)                                   cATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCtt                                       <  2:268660/81‑1 (MQ=255)                                    aTTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCttt                                      <  1:268661/81‑1 (MQ=255)                                       ccGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGt                                    <  1:268662/80‑1 (MQ=255)                                         ggggTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGa                                 <  1:268663/81‑1 (MQ=255)                                          gggTACCAAAATGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCg                                  >  2:268685/1‑79 (MQ=255)                                                    aTGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTc                       >  2:268701/1‑80 (MQ=255)                                                    aTGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCa                      <  1:268664/81‑1 (MQ=255)                                                     tGGGATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCaa                     <  1:268665/81‑1 (MQ=255)                                                      gggATAGGGCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAAc                    >  1:268686/1‑81 (MQ=255)                                                              gCCTGGTACATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGttt            <  1:268666/81‑1 (MQ=255)                                                                     aCATGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGtcat     >  2:268703/1‑81 (MQ=255)                                                                       aTGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGtcatc    >  2:268692/1‑80 (MQ=255)                                                                        tGGCA‑TTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCaa  <  2:268667/81‑1 (MQ=255)                                                                             |                                                                             AATTTGTTTCGCAAAAGCCAGTTCTTCGGAGTTCCATTCCGGGGTACCAAAATGGGATAGGGCCTGGTACATGGCATTTTCTAAGGTGCGATTCGGGAGATAACTGGAACAGGCTTTGTCGAAGCGGCATTCAACCGTGGTTTCGGTCATCAA  >  NC_000913/1402257‑1402409

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