Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,424	A→T	A208A (GCA→GCT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,424	0	A	T	97.1%	216.1 / NA	70	A208A (GCA→GCT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (34/34);  minor base G (0/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGG‑TCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCA  >  NC_000913/3181‑3666                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gggCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttc                                                                                                                                                                                                                                                >  2:621000/1‑250 (MQ=255)   gCGAGCTGGAAGGACGAATCTCCGGGAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGGTTGATGAGTGGCTGTGGGTGCTGGCGCATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTAGCGCCGCCAGGATTGCATTGCGCACGGGCTACATCTGGCTGGCgtcatt                                                                                                                                                                                                                                             <  2:310401/251‑1 (MQ=255)     gAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttcattc                                                                                                                                                                                                                                            >  2:83417/1‑250 (MQ=255)         tGGAAGGGCGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGGCGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACCGGCGACATCTGGCTGGCTTCATTCACGcc                                                                                                                                                                                                                                       >  2:145299/1‑251 (MQ=255)               gacgaaTCTCCGGCAGCATTCATTACGACAACGTGGCCCCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTAt                                                                                                                                                                                                                                 <  1:422116/246‑1 (MQ=255)                  gaaTCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGAGCACGGGCGACATCTGGCTGGCTTCCTTCACGCCTGCTATTcc                                                                                                                                                                                                                              >  2:499070/3‑251 (MQ=255)                          ggCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGcc                                                                                                                                                                                                                      >  1:528377/1‑251 (MQ=255)                           gcagcaTTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCt                                                                                                                                                                                                                     >  1:346461/1‑251 (MQ=255)                                attcattACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCt                                                                                                                                                                                                                >  2:164701/1‑251 (MQ=255)                                         gACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCTCCATGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGaa                                                                                                                                                                                                       <  1:257455/251‑1 (MQ=255)                                          aCAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCCCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGa                                                                                                                                                                                                        >  2:165209/1‑249 (MQ=255)                                            aaCGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCt                                                                                                                                                                                                    <  1:600536/251‑1 (MQ=255)                                              cGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGGTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCtga                                                                                                                                                                                                  <  1:95315/251‑1 (MQ=255)                                                  gCACCGTGTTTTCGTGGTGGTATGCAGTGGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGAGTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGaa                                                                                                                                                                                              <  2:618930/251‑1 (MQ=255)                                                        tgtTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCCGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACCTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGCGCTTGCCGCGAAGCTGATGCAAGATGt                                                                                                                                                                                        >  2:494094/1‑251 (MQ=255)                                                         gtTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGtt                                                                                                                                                                                       >  1:343246/1‑251 (MQ=255)                                                                         gCAGTTGATGATCGAAGAAAACGACAGCATCAGCCAGGAAGGGCCAGGGTTGGATGAGTGGCTGTGGGGGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCGAGGATTGCATTGCGCGCGGGCGACATCTGGCGGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTGATCGCTGagctgtacc                                                                                                                                                                       <  2:315092/251‑9 (MQ=255)                                                                          gaGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATGGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGATGCTGATGAAAGATGTTATCGCTGATCCCTAcct                                                                                                                                                                      <  2:40727/250‑2 (MQ=255)                                                                            gTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTg                                                                                                                                                                    <  1:434993/251‑1 (MQ=255)                                                                               gatgatCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAAc                                                                                                                                                                 <  1:375604/251‑1 (MQ=255)                                                                               gatgatCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAAc                                                                                                                                                                 <  2:24657/251‑1 (MQ=255)                                                                                  gatCGAAGAAAACGCCATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGt                                                                                                                                                              >  1:94810/1‑251 (MQ=255)                                                                                   ggcGAAGAAAACGACATAATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGACGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTTCCGGCGCAGTATCGCCGCTAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCTGCATGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGtt                                                                                                                                                             <  2:282642/249‑1 (MQ=255)                                                                                   atCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGtt                                                                                                                                                             <  1:608831/251‑1 (MQ=255)                                                                                         gaaAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGcc                                                                                                                                                       >  1:587062/1‑251 (MQ=255)                                                                                             aCGACATCAGCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGc                                                                                                                                                   <  1:129437/251‑1 (MQ=255)                                                                                                 catcatCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGCCGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTcc                                                                                                                                               >  1:94390/1‑251 (MQ=255)                                                                                                 catcatCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTcc                                                                                                                                               >  1:498990/1‑251 (MQ=255)                                                                                                       cagccagcTAGTGCCAGGGTCTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCACGACGGCAGAAGCCAG‑GACTATTTTGCCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCAAGGGGGGCATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTAcggcagg                                                                                                                                         <  1:288156/251‑1 (MQ=255)                                                                                                       cagccagcAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcagg                                                                                                                                         <  2:265381/251‑1 (MQ=255)                                                                                                             gcAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggc                                                                                                                                   >  1:36562/1‑251 (MQ=255)                                                                                                              cAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggc                                                                                                                                   >  1:300670/1‑250 (MQ=255)                                                                                                                aGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggcagg                                                                                                                                >  1:156231/1‑251 (MQ=255)                                                                                                                  ggCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGGATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCGACCGTGAACGGGTACGGCCTGGCTTCcggcaggctcggcaggcg                                                                                                                              <  2:449074/250‑1 (MQ=255)                                                                                                                  agCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTACAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGCCCTCTGGCTGGCTTCCTTCACGCCTGCTATTCCCGTCAGCCTGCGCTTGCCGCGTAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcgggcgcggccggcg                                                                                                                              >  2:83134/2‑251 (MQ=255)                                                                                                                       gggTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑ttgc                                                                                                                        <  1:572450/251‑4 (MQ=255)                                                                                                                              atgaGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGTTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATc                                                                                                                 <  1:337627/251‑1 (MQ=255)                                                                                                                                      cTGTGTGTGCTGGCGTCTCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCAGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGCTTG‑GGAAATAGGCGCGGt                                                                                                         >  2:224215/1‑251 (MQ=255)                                                                                                                                         tgGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGCTTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAAACGGCGCGGTagg                                                                                                      >  1:415338/1‑250 (MQ=255)                                                                                                                                            gTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTagcgag                                                                                                   <  1:183171/251‑1 (MQ=255)                                                                                                                                                    gTATCCGGGGATTAAAGTTGCGCGGGGAGAAGCCAGTGGGT‑GTTTACCGGCGCAGTATCGCCGCCAGGATTTCATTGCGCACGGGCGACATCTGGCGGTCTTCACTCACGCCTGGTATTCCCGTCAGCCTGATCTTGCCGCGATGCTGATGAAATATGTTATCGCTGATCCCGACCTTGAACGGGTACTGCCTGGGTTGCGGCAGGCGCGGCCGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGGGGTAtct                                                                                           <  2:384400/251‑1 (MQ=255)                                                                                                                                                    gTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGGAAACGGCGCGGTagcgggggggatct                                                                                           >  2:473555/1‑241 (MQ=255)                                                                                                                                                      aTCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATctcc                                                                                         <  2:36562/251‑1 (MQ=255)                                                                                                                                                        ccGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATctccgg                                                                                       <  1:387706/251‑1 (MQ=255)                                                                                                                                                                        cgtcgtgaGAAGCCTG‑GGCTATTTGAGCGGCGCGGTATCGCCGCCAGGATTGCATTGCGCACGGGCGGCTTCTGGCTGGCTTCTTTCAAGCCTGCTAGTCCCGTCAGACTGAGCTTGCCGCGAAGATGATGAAAGATGTTATCGCTGATCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCTGCAATCCGCGCGGTAGCGAGCGGTATCTCCGTCTCCGGCCCGACCtt                                                                        <  2:156756/243‑1 (MQ=255)                                                                                                                                                                          aCGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGtt                                                                     >  1:590276/1‑251 (MQ=255)                                                                                                                                                                           cGGCAGAAGCCAG‑GGCTATGTTACCGGCTCAGTATCTCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTAGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCAGGCTCCGGCCCGACCTTGTTc                                                                    <  2:499857/251‑1 (MQ=255)                                                                                                                                                                                 aaGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCtg                                                              >  1:197491/1‑251 (MQ=255)                                                                                                                                                                                  aGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGGAAACGGCGCGGCAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGTTCGCTCtgg                                                             >  1:183010/1‑250 (MQ=255)                                                                                                                                                                                      aG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCtgtgtg                                                          >  1:285957/1‑250 (MQ=255)                                                                                                                                                                                       g‑ggCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCCTTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGCTGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATAGGCGCGGTAGCGAGCGGTATCCCCGGGTCCGGCCCGACCTTGTTCGGTCTGTGTGAc                                                        >  2:207776/1‑251 (MQ=255)                                                                                                                                                                                         ggCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCAGTGCGCACGGGCGACATCTGGCTGGCTGCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTGACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTTTGACa                                                       <  2:325684/251‑1 (MQ=255)                                                                                                                                                                                             aTTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGCCATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGc                                                    <  2:346461/250‑1 (MQ=255)                                                                                                                                                                                             aTTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCGTGTTCGCTCTGTGTGACAAGc                                                    <  1:207776/250‑1 (MQ=255)                                                                                                                                                                                                   ccGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGatac                                             <  1:401876/251‑4 (MQ=255)                                                                                                                                                                                                         cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGATGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGCACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑GTGTGGAACTCGGCGCGGTAGCCAGCGGTATCTCCTGCTCCGGCCCCACCTTGTTTGCTCTGTGTGACAAGCCGGATACCGcccc                                       >  1:353812/1‑250 (MQ=255)                                                                                                                                                                                                                gccgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTTCGGGAATCGGCGCGGTTGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTTTTCGCTCTGTGTGACAAGCCGGaatcccccccgcccgtt                                >  2:615222/1‑236 (MQ=255)                                                                                                                                                                                                                   gccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTCGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGcc                            >  1:453054/1‑251 (MQ=255)                                                                                                                                                                                                                       ggATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACt                        >  1:511198/1‑251 (MQ=255)                                                                                                                                                                                                                        gcttgaattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGGAAACGGCGCGGTAGCGAGAGGGATCTCCGGCTCCGGCCCGACCTTGTTTGCCCTTTGTGACAAACCGGGTAACGCCCAG‑GCGCTTGCCGACTg                       >  2:408943/7‑251 (MQ=255)                                                                                                                                                                                                                        gATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTg                       <  2:285957/251‑1 (MQ=255)                                                                                                                                                                                                                         attgcattgcGCACGGGCGACCTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGCCCTTGTTCGCTCTGTGTGCCAAGCCGGATACCGCCCCGCGCG‑TTGCCGACTgg                      >  1:207131/1‑251 (MQ=255)                                                                                                                                                                                                                             cattgcGCACGGGCGACATCGTGCTGGCTTCATTCACGCCTGACATTCCCGTCAGCCTGAGCTTGCCGCTAAGCTGATGAAAGATGTTAGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCACGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTgg                  <  2:165440/251‑1 (MQ=255)                                                                                                                                                                                                                              attgcGCACGGGCGACTGCTGGCTGGCTTCATTCACGCCTGTTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGAGGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCAGCTCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCA‑TTGCCGACTGGTTggg                 <  2:459118/251‑1 (MQ=255)                                                                                                                                                                                                                               ttgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCCGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGt                >  1:297839/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggCGACATCTTGCTGGCTTCATTCACTCCTGCTATTCCCGGCATCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTGCTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATTCCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTAc       <  2:503356/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGGTCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTAc       <  1:60601/251‑1 (MQ=255)                                                                                                                                                                                                                                           gACATCTGGTTGGCTTCATTCACGCCTGCTATTCCCTGCAGCCTGAGCGTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTTACAAGCCGGATACCGCCCAGCGCG‑TTGTCGACTGGTTGGGTAAGAACTACCTg    <  2:580093/251‑1 (MQ=255)                                                                                                                                                                                                                                             cATCGGGCTGGCTTCATTCACGCCTGCTAGTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACGGTTTGGGTAAGAACTACCTGCa  <  2:207131/251‑1 (MQ=255)                                                                                                                                                                                                                                              aTCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGG‑TTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCG‑TTGCCGACTGGGTGGGTAAGAACTACCTGCa  >  1:493145/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAG‑GGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGG‑TCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCG‑TTGCCGACTGGTTGGGTAAGAACTACCTGCA  >  NC_000913/3181‑3666

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