Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,154,109	G→A	*270* (TAG→TAA)	pabC →	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,154,109	0	G	A	100.0%	185.1 / NA	59	*270* (TAG→TAA)	pabC	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/28);  total (31/28)

TCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGATATTT  >  NC_000913/1153975‑1154237                                                                                                                                       |                                                                                                                                 tCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCAt                                                                                                                              >  1:816283/1‑139 (MQ=255) tCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCAt                                                                                                                              >  1:275756/1‑139 (MQ=255)           cagATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAgtg                                                                                                                    <  2:275756/139‑1 (MQ=255)                  gatgGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGa                                                                                                             <  2:847224/139‑1 (MQ=255)                      gTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAAt                                                                                                         <  1:62825/139‑1 (MQ=255)                               aaTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGTATTAATTATGAAAAAAGTGTTATTGATAATCTTGTTAtt                                                                                                <  2:1229061/139‑1 (MQ=255)                               aaTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTAtt                                                                                                <  2:816283/139‑1 (MQ=255)                                               tGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGt                                                                                >  1:552297/1‑139 (MQ=255)                                                  tGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTAt                                                                              <  2:39126/138‑1 (MQ=255)                                                  tGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTAt                                                                              >  1:39126/1‑138 (MQ=255)                                                   gCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           >  2:1146688/1‑124 (MQ=255)                                                   gCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           <  1:1146688/124‑1 (MQ=255)                                                        tATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTa                                                                                      <  2:757892/124‑1 (MQ=255)                                                        tATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTa                                                                                      >  1:757892/1‑124 (MQ=255)                                                             gCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           <  1:327374/114‑1 (MQ=255)                                                             gCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           >  2:327374/1‑114 (MQ=255)                                                                  tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGAt                                                                                                            <  1:412736/92‑1 (MQ=255)                                                                  tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGAt                                                                                                            >  2:412736/1‑92 (MQ=255)                                                                  tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCg                                                             >  2:897779/1‑139 (MQ=255)                                                                   ggCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGt                                                            >  2:470125/1‑139 (MQ=255)                                                                                 ttCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTAt                                                                                                                >  1:663954/1‑73 (MQ=255)                                                                                 ttCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTAt                                                                                                                <  2:663954/73‑1 (MQ=255)                                                                                 ttCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGcc                                              <  1:662438/139‑1 (MQ=255)                                                                                 ttCGTCAGCAACGTTATATGAATATTTAGCCCCACTCTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGcc                                              <  1:897779/139‑1 (MQ=255)                                                                                   cGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCAt                                            >  2:1080371/1‑139 (MQ=255)                                                                                       aGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                               <  1:386787/116‑1 (MQ=255)                                                                                       aGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                               >  2:386787/1‑116 (MQ=255)                                                                                           aCGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGa                                    <  1:1080371/139‑1 (MQ=255)                                                                                                      atatTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGc                         >  2:275689/1‑139 (MQ=255)                                                                                                      atatTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGGTATTGATAATCTTGTTAtt                                                                                                <  1:954529/68‑1 (MQ=255)                                                                                                      atatTTAGCCCCACTGTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTAtt                                                                                                >  2:954529/1‑68 (MQ=255)                                                                                                       tatTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGtgtg                                                                 <  2:998091/98‑1 (MQ=255)                                                                                                       tatTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGtgtg                                                                 >  1:998091/1‑98 (MQ=255)                                                                                                       tatTTAGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCt                        >  2:380576/1‑139 (MQ=255)                                                                                                            aGCCCCACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCa                   >  2:544864/1‑139 (MQ=255)                                                                                                              ccccACTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGtt                                                                                                                  <  1:1141332/42‑1 (MQ=255)                                                                                                              ccccACTTTGTGAGCGCCCGAATTAATAATGAAAACAGTGtt                                                                                                                  >  2:1141332/1‑42 (MQ=255)                                                                                                                  aCTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           >  1:91575/1‑61 (MQ=255)                                                                                                                  aCTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           <  2:91575/61‑1 (MQ=255)                                                                                                                   cTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           >  1:892043/1‑60 (MQ=255)                                                                                                                   cTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgc                                                                                           <  2:892043/60‑2 (MQ=255)                                                                                                                   cTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTg                                                                      >  2:885487/1‑81 (MQ=255)                                                                                                                   cTTTGTGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTg                                                                      <  1:885487/81‑1 (MQ=255)                                                                                                                        tgAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGc                                                                                              >  1:218982/1‑52 (MQ=255)                                                                                                                        tgAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGc                                                                                              <  2:218982/52‑1 (MQ=255)                                                                                                                        tgAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGtgtg                                                                 <  1:683618/81‑1 (MQ=255)                                                                                                                        tgAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGtgtg                                                                 >  2:683618/1‑81 (MQ=255)                                                                                                                        tgAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGa       >  1:274472/1‑139 (MQ=255)                                                                                                                            cgcCCGAATTAATCATGAAAAAAGTGTTATTGATAAt                                                                                                         >  2:1087601/1‑37 (MQ=255)                                                                                                                            cgcCCGAATTAATCATGAAAAAAGTGTTATTGATAAt                                                                                                         <  1:1087601/37‑1 (MQ=255)                                                                                                                            cgcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGc                                       <  1:273167/103‑1 (MQ=255)                                                                                                                            cgcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGc                                       >  2:273167/1‑103 (MQ=255)                                                                                                                            cgcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCa                               <  1:1212933/111‑1 (MQ=255)                                                                                                                            cgcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCa                               >  2:1212933/1‑111 (MQ=255)                                                                                                                             gcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           <  1:436349/50‑1 (MQ=255)                                                                                                                             gcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                           >  2:436349/1‑50 (MQ=255)                                                                                                                             gcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTAt                                                                              <  1:124058/63‑1 (MQ=255)                                                                                                                             gcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTAt                                                                              >  2:124058/1‑63 (MQ=255)                                                                                                                             gcCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGATAttt  >  1:861746/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                 TCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGATATTT  >  NC_000913/1153975‑1154237

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