Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,149,506	T→C	S251P (TCC→CCC)	pabC →	aminodeoxychorismate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,149,506	0	T	C	100.0%	107.0 / NA	29	S251P (TCC→CCC)	pabC	aminodeoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/16);  total (13/16)

TCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATC  >  NZ_CP009273/1149421‑1149591                                                                                      |                                                                                      tctGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCttt                                                                                  <  1:25756/91‑1 (MQ=255) tctGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCttt                                                                                  <  1:284629/91‑1 (MQ=255)    ggAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCg                                                                               <  2:704307/91‑1 (MQ=255)       agagTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCa                                                                            <  1:525976/91‑1 (MQ=255)        gagTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAg                                                                           >  2:838922/1‑91 (MQ=255)            cGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAAc                                                                       >  2:806844/1‑91 (MQ=255)            cGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAAc                                                                       <  2:143577/91‑1 (MQ=255)              ttGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGt                                                                     <  1:182066/91‑1 (MQ=255)                  aggcagATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTtata                                                                 >  2:13833/1‑91 (MQ=255)                     cagATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGa                                                              <  2:225512/91‑1 (MQ=255)                       gatgagatgGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAat                                                            >  2:423501/1‑91 (MQ=255)                         tgagatgGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAatat                                                          <  1:468165/91‑1 (MQ=255)                             atgGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAg                                                      <  2:530783/91‑1 (MQ=255)                               ggTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGcc                                                    <  1:332886/91‑1 (MQ=255)                                 ttATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGcccc                                                  <  1:110516/91‑1 (MQ=255)                                 ttATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGcccc                                                  <  1:160068/91‑1 (MQ=255)                                         aaTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTtgtg                                          >  1:139761/1‑91 (MQ=255)                                          aTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGa                                         >  2:309073/1‑91 (MQ=255)                                          aTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGa                                         >  2:350/1‑91 (MQ=255)                                                      cAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTa                             >  1:252978/1‑91 (MQ=255)                                                             gCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATAttt                                                        >  2:163627/1‑57 (MQ=255)                                                             gCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATAttt                                                        <  1:163627/57‑1 (MQ=255)                                                                   aTGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAg                >  1:400025/1‑91 (MQ=255)                                                                       gCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGtt            >  1:764529/1‑91 (MQ=255)                                                                       gCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGtt            <  1:853205/91‑1 (MQ=255)                                                                           gtgGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTg        >  1:578819/1‑91 (MQ=255)                                                                           gtgGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTg        >  2:413196/1‑91 (MQ=255)                                                                                 aTGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATc  <  1:568407/91‑1 (MQ=255)                                                                                 aTGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATc  <  1:106821/91‑1 (MQ=255)                                                                                      |                                                                                      TCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATC  >  NZ_CP009273/1149421‑1149591

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