Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,317,129	G→A	Q32* (CAG→TAG)	trpA ←	tryptophan synthase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,317,129	0	G	A	100.0%	100.1 / NA	29	Q32* (CAG→TAG)	trpA	tryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/19);  total (10/19)

CATCAATGCCTTTGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCG  >  NC_000913/1316886‑1317373                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     catcaATGCCTTTGTTAAACACCAGATGGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATg                                                                                                                                                                                                                                               <  1:301581/251‑1 (MQ=255)    caaggCCTTTGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCGGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGccc                                                                                                                                                                                                                                            <  1:24279/247‑1 (MQ=255)             tGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGGGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCAGCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCATTTAGCGTATCGATAAGTTTCAATGACTACTCAATGCCCGGATCAcc                                                                                                                                                                                                                                    <  1:36725/250‑1 (MQ=255)                           ttGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAg                                                                                                                                                                                                                      <  1:36166/250‑1 (MQ=255)                                      cAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATgcg                                                                                                                                                                                                           <  1:150228/250‑1 (MQ=255)                                          aGGCCAATGGGCATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGCTCACCGAGCGTGACGAAAGGAACGAATGCGCcttc                                                                                                                                                                                                      >  1:294898/1‑251 (MQ=255)                                                                     tGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCAGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAAc                                                                                                                                                                             >  2:270340/1‑249 (MQ=255)                                                                                 gcaaGCATTTAAAACCATTGGGCCGGATACCCCCCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACGGGGCAAACAGAGATTCGTAGCg                                                                                                                                                               <  2:164074/248‑1 (MQ=255)                                                                                 gCCAGCAGTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCAGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTAGCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCGTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCg                                                                                                                                                               <  1:270340/251‑1 (MQ=255)                                                                                     gCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTcc                                                                                                                                                           <  1:334764/251‑1 (MQ=255)                                                                                            aaaaCATTGGGCCGGAGCCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGAt                                                                                                                                                    <  2:64983/251‑1 (MQ=255)                                                                                                                 cccTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAAAGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAat                                                                                                                               <  2:282611/250‑1 (MQ=255)                                                                                                                                 cgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGa                                                                                                                 >  1:112120/1‑249 (MQ=255)                                                                                                                                      agTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTc                                                                                                          <  2:93422/251‑1 (MQ=255)                                                                                                                                         ggCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCtt                                                                                                        <  1:220662/250‑1 (MQ=255)                                                                                                                                                 gAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATcgccgc                                                                                               >  1:50158/1‑251 (MQ=255)                                                                                                                                                                   ccAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAacca                                                                             <  1:119639/251‑1 (MQ=255)                                                                                                                                                                    cAGTGGGTCGGCGCAGGGGCTACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGAGCTCATTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAaccac                                                                            >  2:374830/1‑251 (MQ=255)                                                                                                                                                                             ggAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCt                                                                                                                                             >  1:288367/1‑177 (MQ=255)                                                                                                                                                                             ggAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTAACGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCt                                                                                                                                             <  2:288367/177‑1 (MQ=255)                                                                                                                                                                                              aCTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTATCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGtttt                                                  <  1:247305/251‑1 (MQ=255)                                                                                                                                                                                                        gTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAgctgct                                                                 >  2:171438/1‑226 (MQ=255)                                                                                                                                                                                                        gTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAgctgct                                                                 <  1:171438/226‑1 (MQ=255)                                                                                                                                                                                                             cACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCa                                   <  2:50158/251‑1 (MQ=255)                                                                                                                                                                                                              aCCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCGCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGTCCGGAAAGGTTAACCACCAGTAGCTGCGCTTTATCCGGGTTTTCGCGCATCATTTTCaa                                  <  2:143380/251‑1 (MQ=255)                                                                                                                                                                                                                           aGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTGTAGCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAgg                     <  1:382776/251‑1 (MQ=255)                                                                                                                                                                                                                                   gATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTggg              >  1:377862/1‑250 (MQ=255)                                                                                                                                                                                                                                       aTTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGAt         >  1:399743/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGACAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAgcg  >  2:224973/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CATCAATGCCTTTGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCG  >  NC_000913/1316886‑1317373

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