Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,895,120	G→A	A106T (GCA→ACA)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,895,120	0	G	A	100.0%	72.9 / NA	21	A106T (GCA→ACA)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/12);  total (9/12)

TACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCT  >  NC_000913/1895043‑1895196                                                                              |                                                                             tACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCacac                                                                           >  2:301874/1‑81 (MQ=255)   cAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCacac                                                                           <  2:301841/79‑1 (MQ=255)   cAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTg                                                                         <  1:301840/81‑1 (MQ=255)         cTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGgtt                                                                    <  1:301842/80‑1 (MQ=255)            cagcagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGttt                                                                <  2:301843/81‑1 (MQ=255)                agGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCt                                                            <  1:301844/81‑1 (MQ=255)                agGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCt                                                            <  2:301845/81‑1 (MQ=255)                  gTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTAc                                                          <  2:301846/81‑1 (MQ=255)                             cgcAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGgccg                                               <  1:301847/81‑1 (MQ=255)                                  aCATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtt                                          <  2:301848/81‑1 (MQ=255)                                    aTTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtttt                                        >  2:301877/1‑81 (MQ=255)                                    aTTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtttt                                        >  2:301870/1‑81 (MQ=255)                                       cGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                     >  1:301872/1‑81 (MQ=255)                                        gCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGt                                    >  2:301878/1‑81 (MQ=255)                                           cAACGCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCAc                                 <  2:301849/81‑1 (MQ=255)                                               gCATAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGc                              >  1:301871/1‑80 (MQ=255)                                                  tAACGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAg                           >  2:301876/1‑80 (MQ=255)                                                    aCGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGaa                         <  1:301850/80‑1 (MQ=255)                                                     cGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAt                       <  2:301852/81‑1 (MQ=255)                                                     cGAAGATTTGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAt                       <  1:301851/81‑1 (MQ=255)                                                            ttGCCATTTCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGaa                >  2:301861/1‑81 (MQ=255)                                                                  tttCAGGGCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGa           >  1:301873/1‑80 (MQ=255)                                                                       gggCGGCACACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGt     >  2:301875/1‑81 (MQ=255)                                                                         gcggcaCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTc   <  2:301853/81‑1 (MQ=255)                                                                          cggcaCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCt  <  2:301854/81‑1 (MQ=255)                                                                              |                                                                             TACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCT  >  NC_000913/1895043‑1895196

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