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	93.5%	125.7 / ‑3.5	47	A208A (GCA→GCT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (25/18);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTT  >  NC_000913/3186‑3648                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     aGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCGTCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCAc                                                                                                                                                                                                                          <  1:163148/251‑1 (MQ=255)    tGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTCTGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGcc                                                                                                                                                                                                                       >  1:121678/1‑251 (MQ=255)        aggacgaaTCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCt                                                                                                                                                                                                                   >  1:302839/8‑251 (MQ=255)           acgaaTCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTAtt                                                                                                                                                                                                                >  1:9166/5‑251 (MQ=255)            cgaaTCTCCGGCAGAATTCATTACGACAACGTGGCACCGTGTTTTCTTGGCGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGAGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGCCGGCACAAGCCAGGGATATTTTCCCGGCGCCGTATCGCCCCCAGGGTTGCAATGCGCACGGGCGACATCTGGCTGGCTTCATTCAACCCCGCTATTc                                                                                                                                                                                                               >  2:178858/4‑251 (MQ=255)             gaaTCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCAGCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTcc                                                                                                                                                                                                              <  1:32527/249‑1 (MQ=255)              aaTCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTccc                                                                                                                                                                                                             >  1:40322/2‑251 (MQ=255)                               attACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑cc                                                                                                                                                                                          <  1:50074/251‑1 (MQ=255)                                            ggCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCAGCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCtgatga                                                                                                                                                                             <  1:282069/251‑1 (MQ=255)                                                  gtgtTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGCTGAAAGATg                                                                                                                                                                       >  1:195284/1‑251 (MQ=255)                                                       ttCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTAt                                                                                                                                                                   <  1:185142/250‑1 (MQ=255)                                                          ttggtggtATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCt                                                                                                                                                               >  1:251038/3‑251 (MQ=255)                                                           tggtggtATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCTCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCt                                                                                                                                                               <  1:355662/249‑1 (MQ=255)                                                           tggtggtATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGc                                                                                                                                                                <  1:378205/248‑1 (MQ=255)                                                              gggtATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAAc                                                                                                                                                           <  1:269587/250‑1 (MQ=255)                                                                     cAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCg                                                                                                                                                    >  1:4234/1‑251 (MQ=255)                                                                             gatCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGCGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCCGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGATGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGt                                                                                                                                            >  1:16273/1‑251 (MQ=255)                                                                             gaacGAAGAAAACGACATCATCAGCCAGCAAGTGACAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGTCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCGGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGt                                                                                                                                            <  1:40696/248‑1 (MQ=255)                                                                              atCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTACTGGCGTATCCGGGGATTAAAGTATCGACGGCAGAAGCCAGGGCTATTTTAACGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGtt                                                                                                                                           >  1:217170/1‑251 (MQ=255)                                                                                      aaaCGACATCATCAGCCAGCAAGTGCCAGGGTATGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTATCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATAGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAAGCTTCATTCACGCCTGCTATGCCCGGCAGCCT‑AGACTTG‑CCGAGAAGGTGATGACAGATGTTATCGGTGATACCTACCGTGAACGGGTAATGcc                                                                                                                                     >  2:208635/1‑249 (MQ=255)                                                                                       aacgccagcatcaGCCATCAAGTGCCAGGTGTTGATGAGTGGCTGTGGGTGGTGGAGTTGTCGGGGATTAAAGTCTCGACGGCAGAAGTCAGGGCTATGTTACCGGCGCAGGATCGACGCCAGGATTGCATTGCGCACGTGCGACATCTGGCTGGTTTCATTCACGCCTGCTATTCTCGTCAGCCTGAG‑CTTG‑CCGGGAAGCTGATGAAAGATGTTATCGCTGAACCTTACCTTGAACGGTTCCTGcctgg                                                                                                                                  <  2:60349/243‑4 (MQ=255)                                                                                       aaCGACATCCTCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTATCGCAGGCAGAAGCCAGGGCTATTTTACCGGCGCGGTATCGCCGCAAGGATTGCATTGCGCAAGGGCGACATCTGGCTGGCTTCATTCAAGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGAGAAGCTGTTGCTAGATCTTATCGATGAACCCTACCGGGAGCGGTTACTGcctcg                                                                                                                                  >  2:68650/1‑248 (MQ=255)                                                                                       aaCGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGCTTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACACTACCGTGAACGGTTACTGcctgg                                                                                                                                  >  1:39608/1‑248 (MQ=255)                                                                                         cGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCTGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCt                                                                                                                                <  1:75637/251‑1 (MQ=255)                                                                                           aCATCATCGGCCAGCTAGTGCCATGGTTTGATGATTGGCTGTGGGGGCTGGCGTTTCCGGGTATTAAAGTCTCGACGGTATAAGCCAGTGTTATTTTACCGGCGTAGTATCGCCGTCAGGATTGCCTTTCGCACGAGCGATATCTGGCTGGCTTCAGTCGCGCCTGCTATTCCCGTCCGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTc                                                                                                                              <  2:237543/251‑1 (MQ=255)                                                                                            catcatCAACCCGCAAGAGCCAGGGTTTTATGAGTTGCTGTGGGAGATGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGAGACATCTGTCAGGCTTCATTACCGCCTGCTATACCCGTAAGCCTGAG‑CTTG‑CCGCGAAGCTGATGCAAGATTTTATCGCTGAACCCTCCCGTAAACGGTTACTGCATGGCTTcc                                                                                                                             >  2:194377/1‑251 (MQ=255)                                                                                                   agccagcAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggc                                                                                                                      >  1:340357/1‑251 (MQ=255)                                                                                                               ccAGGGTTTGATGAGTGGCTGTGGGTGCTGGAGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATATGGATTGCTTCATTACCGCATGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGCAGCTGTTGAAAGATGTTATCGCTGAACCCTACCGTGAACAGTTACTGCCTCGCTTCcggcaggcagggaatgcggg                                                                                                          >  2:353408/1‑239 (MQ=255)                                                                                                                     tttGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCggg                                                                                                    >  1:258904/1‑250 (MQ=255)                                                                                                                               ggCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGcgcg                                                                                          <  1:280007/251‑1 (MQ=255)                                                                                                                                   gtgGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGt                                                                                        >  1:154166/1‑249 (MQ=255)                                                                                                                                      ggTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGACGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGgcgggcgggtggggggataggcgcggtagcgg                                                                                   >  1:370125/1‑227 (MQ=255)                                                                                                                                                 aTCCGGGGATTAAAGGCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TAtctc                                                                        <  1:35068/250‑1 (MQ=255)                                                                                                                                                           tAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATAGGCGCGGGAGCGAGCGG‑GATCCCCGGCTCCGGcc                                                             >  1:207984/1‑251 (MQ=255)                                                                                                                                                                             aGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGAAC‑TTGTTCGCTCtgt                                          >  1:144861/1‑251 (MQ=255)                                                                                                                                                                                         tgtcccggtgcagttccgccgcCATGATTGAATTGCGCACGGGTGCCATCTGGCTGGCTTCATTCTCGCCTGCTATCCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTAGCGCTGATCCTTACCGTGAACGGTTACTGCCTGGTTTCCGGCAGGAGCGGCCGTCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TAGCTCCGGCTCCGGCCCGACC‑TTGTTCGCGATGTGTGACAAgacg                               <  2:5604/234‑4 (MQ=255)                                                                                                                                                                                            aCCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGATAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGat                            <  1:207840/250‑2 (MQ=255)                                                                                                                                                                                                 cgcAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGAGGGCGG‑GATCTCCGGATCCGGCCCGACC‑TTGGTCGCTCTGGGTGTCAAGCCGGGAACCGcc                      >  1:265378/1‑251 (MQ=255)                                                                                                                                                                                                   cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑GATCTCCGGCCCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATAACGcccc                    >  1:338857/1‑250 (MQ=255)                                                                                                                                                                                                       aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTTCGGAAATCGGCGCGGGAGCGAGCGG‑TATCTCCGGCTCCGGCCCGGCC‑TTGTTCGCTCTGTGTGACAAGCCGGAATACGCCCAgcgc                >  1:196253/1‑251 (MQ=255)                                                                                                                                                                                                       aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATAGGGGCGGGAACGAGAGG‑AATCTCCGGCTCCGGCCCGACC‑ATGTTCGCTCTGTGTCACAAGACGGATACCGCCCAgcgc                >  1:261877/1‑251 (MQ=255)                                                                                                                                                                                                          gccgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGtt             <  1:178858/251‑1 (MQ=255)                                                                                                                                                                                                           cctccAGGATTGCATTGCGCTCGGGCGACAGTTGGCTGGCTTCATTCACGCCGGCGATTCCCGTCAGCCTGAG‑CTTG‑CCGGGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCGGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCTGTAGCGAGCGG‑TATCTCCGGCTCCGTCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTg            <  2:368578/248‑1 (MQ=255)                                                                                                                                                                                                            cgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCAGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTCTCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGGAACGGTGCGGTAGAGCGCGGGCATC‑CCGGCTCCGGGCCGACCATTGTT‑GGTCCGGGTGACAAGCCGGagcacgcccaacgccgtgg           >  2:77237/1‑233 (MQ=255)                                                                                                                                                                                                                   attgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCCGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGTTGCCGGCTgg    >  1:122527/1‑251 (MQ=255)                                                                                                                                                                                                                      gcatggcGCACGGGTGACCTCTGGCGGGCTTCAGTCACGCCTGCTTTTCCCGGCCGCCTGAG‑‑TTGTCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCGCCGGCCCGACC‑TGGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGtt  <  2:338857/250‑1 (MQ=255)                                                                                                                                                                                                                               aCGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑tt                                                     <  1:60431/190‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     AGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAG‑CTTG‑CCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTT  >  NC_000913/3186‑3648

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