Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,483,721	G→C	R156R (CGC→CGG)	pabA ←	aminodeoxychorismate synthase component 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,483,721	0	G	C	100.0%	97.1 / NA	30	R156R (CGC→CGG)	pabA	aminodeoxychorismate synthase component 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/16);  total (14/16)

CAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACC  >  NZ_CP009273/3483634‑3483799                                                                                        |                                                                               cAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACtggcgatgccg                                                                              >  1:474197/1‑90 (MQ=255)  aGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACtggcgatgccga                                                                             <  1:204653/90‑1 (MQ=255)   ggAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGaa                                                                            <  2:612607/90‑1 (MQ=255)    gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  2:301661/90‑1 (MQ=255)    gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  2:49364/90‑1 (MQ=255)    gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  2:283502/90‑1 (MQ=255)    gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAAt                                                                           <  1:286714/90‑1 (MQ=255)     aaaTTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATc                                                                          <  2:400941/90‑1 (MQ=255)        ttAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATcccc                                                                       >  1:489555/1‑90 (MQ=255)              agcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAAtc                                                                 <  2:567429/90‑1 (MQ=255)              agcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAAtc                                                                 <  2:541970/90‑1 (MQ=255)               gcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAAtct                                                                >  1:249009/1‑90 (MQ=255)                    tGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCgg                                                           >  2:497144/1‑90 (MQ=255)                         tCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTc                                                      <  1:239896/90‑1 (MQ=255)                         tCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTc                                                      <  1:541140/90‑1 (MQ=255)                            ttGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCt                                                   >  1:319141/1‑90 (MQ=255)                                                 tggatggaACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCa                              >  2:501878/1‑90 (MQ=255)                                                    atggaACTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGc                           <  2:401808/90‑1 (MQ=255)                                                          cTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTaa                     >  2:537211/1‑90 (MQ=255)                                                          cTGCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTaa                     <  2:489555/90‑1 (MQ=255)                                                            gCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGTTAATg                   >  2:44086/1‑90 (MQ=255)                                                            gCACACCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATg                   <  2:603057/90‑1 (MQ=255)                                                                acCTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTc               <  1:78007/90‑1 (MQ=255)                                                                 ccTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGc                               >  2:417428/1‑73 (MQ=255)                                                                 ccTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGc                               <  1:417428/73‑1 (MQ=255)                                                                 ccTTCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCa              >  1:359142/1‑90 (MQ=255)                                                                   ttCCAGATCCCACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAgg            >  1:415342/1‑90 (MQ=255)                                                                           cccACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccacca    >  2:515105/1‑90 (MQ=255)                                                                           cccACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccacca    >  2:165162/1‑90 (MQ=255)                                                                             cACTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccaccacc  >  2:287069/1‑90 (MQ=255)                                                                               cTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTcc                                                 <  2:548021/41‑1 (MQ=255)                                                                               cTGGCGATGCCGAATCCCCATAATCTCTCGGGTTTCGCTcc                                                 >  1:548021/1‑41 (MQ=255)                                                                                        |                                                                               CAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACC  >  NZ_CP009273/3483634‑3483799

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