Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,593	A→G	A419A (GCA→GCG)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,593	0	A	G	91.7%	121.3 / ‑2.5	49	A419A (GCA→GCG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (24/20);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTG  >  NC_000913/1359‑1834                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   cgcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCCGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAata                                                                                                                                                                                                                                    <  1:405549/251‑1 (MQ=255)  gcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAatat                                                                                                                                                                                                                                   >  1:375113/1‑251 (MQ=255)      gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCGCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAc                                                                                                                                                                                                                               <  2:92439/251‑1 (MQ=255)        aTTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACa                                                                                                                                                                                                                              <  1:121379/250‑1 (MQ=255)         tttCCGTGGTGCTGATTACGCAATCATCTTCCGAATCCAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGTAGGGGGTGGGATGCGCACCT‑TGCGAGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCCATATCAACAtt                                                                                                                                                                                                                            >  2:61588/1‑251 (MQ=255)                    cTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGTGGGAGGGGATGGTAATAGCAACTGTG‑GTGGGATATCGGCGCAATACTTTTACGCCGTGGCCCGCGCCCATATCCACATTGTTGccctagt                                                                                                                                                                                                                 >  2:21150/1‑246 (MQ=255)                    cTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGGGGCAATATCAACATTTTCGCCATTGc                                                                                                                                                                                                                 >  1:81001/1‑251 (MQ=255)                      gATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTAGTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGc                                                                                                                                                                                                                 <  1:144553/249‑1 (MQ=255)                       gtttcGCAATCATCTTTCGTATACTGGATCAGTTTCGGCGTTCCGCATAGCGATGGGGTGCGAGCTAAACGGGCAATGCAGGAAGATTTCGACGTGGAAGTGACAGAAGGCTGACTGGAGGCGTAGGCGGTGACGGAACGGCTGGCCATGATCGCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAGTATCAACATTGTCGCCATTGCTCa                                                                                                                                                                                                              <  2:239051/247‑1 (MQ=255)                            gCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAt                                                                                                                                                                                                         <  1:49764/251‑1 (MQ=255)                               atcatcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAtctt                                                                                                                                                                                                      >  1:146347/1‑251 (MQ=255)                                tcatcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTTCTCAGGGGtcttc                                                                                                                                                                                                     >  1:481598/1‑251 (MQ=255)                                       ccGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCTGGCCATTATCTCGGTGGTAGGCGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTTCCGCGCTGGCCCGCGCCAATATCAAAATTGTCGCCATTGCTCAGGGATCTTC‑TGaaagc                                                                                                                                                                                             >  1:215168/1‑248 (MQ=255)                                            tACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTTAACGGGCAATGCAGGCAGAGTTCTACCTGGAACTGAAAGAATGCTTACTGGAGCCCCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGCGGCAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCTCCAATATCAACATTGCCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAAt                                                                                                                                                                                        >  1:143812/1‑251 (MQ=255)                                                 taTCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCATTATCTCGGTGGTAGGTGCTGGTATGCGAACCT‑TGCGTGGGGTCTCGGCGAAATTCTTTGCCGCGGTTGCCCGCGCCCATATCAACATTGTCGCCATTTCTCAGGGGTCTTT‑TGCACGCGCCATCTCTg                                                                                                                                                                                   >  2:395493/2‑251 (MQ=255)                                                         tCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCCCAATCTCTGTCGGGGt                                                                                                                                                                            >  2:397278/1‑250 (MQ=255)                                                              gTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATaa                                                                                                                                                                      <  1:184590/251‑1 (MQ=255)                                                                    caAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGAAGGAACGGATGGCCATTATCTCGGGGGTAGGGGATGGTAAGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGGTGGCCCGCGCCAATATCAACATTGCCGCCATTGCACCGGGATCTTC‑TGAACGCTCATTCTCTGTCGTGGTAAATAACcatg                                                                                                                                                                 >  1:363077/1‑250 (MQ=255)                                                                     aaaGCGACTGTGTGCGAGCTGAACGGGAAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGGAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCAGCGATGGCCCGCGCCAAAATCAAACTTGTCGACATTGCTCAGGGATCTTC‑TTAACGCCCAATATCTGTCGAGGAAAATAAAgatgat                                                                                                                                                               >  1:156817/1‑251 (MQ=255)                                                                       aGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAACGCAGGCTTACTGGAGCCGCTGGCGGTGACGGGACGGCTGGCCATTATCTCGGTGGTAGGCGAAGGTATGCGAAACT‑TGCGTGGGATCTCGGCGAAGTTCTTTGCCGCACTGGCCCGCGCCCATATCAAAATTGTCGCCCTTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGGAAATAACGATGATGc                                                                                                                                                             >  2:460917/1‑251 (MQ=255)                                                                           aCTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGAcc                                                                                                                                                         <  1:61588/251‑1 (MQ=255)                                                                                   cgcgCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTATCTCGGTGGTCGGTGGTGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTT‑TTCACGCTCAATCCCTGTCGTGGTAAATAACGGTGATGCGACCACCGGCgg                                                                                                                                                 >  2:289490/4‑250 (MQ=255)                                                                                    gAGCTGAACGGGCAATGCAGGAAGAGTTCTCCCTGGACCTGAAAGAAGGCTTACTGGAGCCGCTGGCGGAGACGGCACGGCTGGCCATTATCTCGGGGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGACGCGCTGGCCCACGCCACTATCAACATTTTACCCATTACTCAGGGAACTTCCTGAATG‑TCAATCTCTGTCATGGTAAATAACgatattgccaccaatagcagg                                                                                                                                                >  2:309720/1‑233 (MQ=255)                                                                                    gAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAAGGAACGGCTGGCAATTATCTCGGTGGGAGGGGATGGTATGGGCACCT‑TGCGTGGGATCTCGGCGAAATTATTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCCTTGCTCAGGGGTCTTC‑TGAACGCTCAATCTCTGTCAGGGGAAATAACGATGTGGCGACCACTGGCgcg                                                                                                                                                >  1:451336/1‑249 (MQ=255)                                                                                             gggCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGTCGGAACGGCTGGCCATTATCTCGGTGGTAGGAGTTGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGAACTTC‑TGAACGCCCAATCTCTGTCGTGGTAAATGACGATGGTGCAAACACTGGCGTGCGCGTTACt                                                                                                                                       >  2:119246/1‑251 (MQ=255)                                                                                               gCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGGGTTACtc                                                                                                                                      >  1:98430/1‑250 (MQ=255)                                                                                                      gggaagaggctggccTGGAACTGAAAGAAGGCTTACTGGTGCCGCTGGCGGTGACGGAACGGTTGGCCATTATCACGGTGGTAGGTGGTGGTATGCGCACCT‑TGCGTGGGATGTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAACTTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAACCTCTGCCGTGGTAAATAACGAGGATGCGACCACTGGCGTGCGCGTTACTCATCAGATg                                                                                                                              <  2:143812/238‑1 (MQ=255)                                                                                                      aGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGCACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATg                                                                                                                              >  1:238868/1‑251 (MQ=255)                                                                                                       ggAAGAGTTCTACCTGGAACTGAAAGCAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCa                                                                                                                                                                                          >  1:364290/1‑190 (MQ=255)                                                                                                       ggAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGGAGGGGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGc                                                                                                                             >  1:219127/1‑251 (MQ=255)                                                                                                       ggAAGAGGTCTCCCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCTGTCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCa                                                                                                                                                                                          <  2:364290/190‑1 (MQ=255)                                                                                                        gAAGAGTTTTACCTGGAACGGAAAGAAGGCTTACTGGAGGCGCGGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCt                                                                                                                            <  2:151720/251‑1 (MQ=255)                                                                                                                       gAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCa                                                                                                             <  1:224728/251‑1 (MQ=255)                                                                                                                           tgatagcagGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGCTTGTATGCGCACCT‑TGCGTGGGTTCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCCACATTGTCGCCATTGCTCAGGGCTCTTC‑TGAACGCTCATTCTCTGTCGTGTTAAATACCGCTGATGCTACCACTGGCGTGCGCGTTTCTCTTCCGATGCTTTTCAATAcctcttcggtt                                                                                                         >  2:192260/8‑242 (MQ=255)                                                                                                                                              gCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTg                                                                                      >  1:425470/1‑251 (MQ=255)                                                                                                                                                                             atcgcGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTgg                                                        <  2:376125/246‑1 (MQ=255)                                                                                                                                                                             aTCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTggg                                                       >  2:292933/1‑250 (MQ=255)                                                                                                                                                                                     ggtAGGTGATGGTATGCGCACCC‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑GGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTGACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGa                                               <  2:471764/251‑1 (MQ=255)                                                                                                                                                                                               ggTATGCGCACCT‑TGCGGGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATACGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACAGATCAGGTTATCGAAGTGTTTGTGATTGGCGACGGTGGCGTTGGCGGTGGGCTGCTGGGACAACTGAAACGTCaaa                                      >  2:253685/1‑248 (MQ=255)                                                                                                                                                                                                 tATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCa                                                                                                                                                        >  1:146296/1‑134 (MQ=255)                                                                                                                                                                                                 tATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCa                                                                                                                                                        <  2:146296/134‑1 (MQ=255)                                                                                                                                                                                                 tATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGATTATCGAAGTGGTTGTGCTTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGGGCCAATGAAGCGTcaacaa                                    >  2:129869/1‑250 (MQ=255)                                                                                                                                                                                                      gcACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAgctgg                              >  1:94100/1‑251 (MQ=255)                                                                                                                                                                                                       cACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAgctggc                             <  2:375113/251‑1 (MQ=255)                                                                                                                                                                                                              gCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaaga                       <  2:94100/251‑1 (MQ=255)                                                                                                                                                                                                                   ggaactcGGCGAAATTCTTTGCCGCGCGGGCCCGCGCCACTATCAACATTGTCGCCATTGCTCAGGGAACTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATaaa                  <  2:124887/247‑1 (MQ=255)                                                                                                                                                                                                                               aaTTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTa      <  2:219127/251‑1 (MQ=255)                                                                                                                                                                                                                                aTTCTTTGCCGCGCTGGCCCGCGCCAATAGCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCGCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCAGCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTAc     <  1:441356/251‑1 (MQ=255)                                                                                                                                                                                                                                     ttGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACgtg  <  1:398521/249‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   CGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCT‑TGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTC‑TGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTG  >  NC_000913/1359‑1834

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