Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,508	A→G	Q724Q (CAA→CAG)	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	2,508	0	A	G	98.6%	222.8 / NA	71	Q724Q (CAA→CAG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/38);  minor base T (0/1);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGT  >  NC_000913/2267‑2753                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      tCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTcgac                                                                                                                                                                                                                                                 <  1:398032/248‑1 (MQ=255) tCGGCAAGCTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGAt                                                                                                                                                                                                                                              >  1:164705/1‑251 (MQ=255)     cAAGTTAGACGAAGGCATTAGTTTCTCCGAAGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGGGATGTTGCCGCTTTTATGGCGAATCTGTCCCAGCTCGACGATctct                                                                                                                                                                                                                                          >  1:493971/1‑251 (MQ=255)       aGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCttt                                                                                                                                                                                                                                        <  1:297099/251‑1 (MQ=255)         ttAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGACATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACATGTGCTGCCCGCAGAGATTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCGCGACGGTCTCTTTGc                                                                                                                                                                                                                                      >  1:367858/1‑251 (MQ=255)                  ggCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCcgcgcgcgg                                                                                                                                                                                                                             >  2:94889/1‑250 (MQ=255)                   gCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTGTCTTGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCGCGTTAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGATGCTGTGCTGCCCGCAGTGTTTAACGCTGAGGGTGATGTTGCCGCTTTGATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTg                                                                                                                                                                                                                            <  1:287129/251‑1 (MQ=255)                   gCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATAACGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTg                                                                                                                                                                                                                            >  1:358085/1‑251 (MQ=255)                                               ggCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAg                                                                                                                                                                                                   <  1:169551/248‑1 (MQ=255)                                                     ggAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCt                                                                                                                                                                                          <  2:141144/251‑1 (MQ=255)                                                      gAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGCTGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTTTGTCGAATCTGTCACAGATCGGCGATATCTTTGCCGCGCGCGTGTCGAAGTCCCGTGTTGAAGGAAAAGTTTTGCGCtt                                                                                                                                                                                         >  2:136296/1‑250 (MQ=255)                                                           gggTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGCAACGGGACGTGAACTGGAGCTGGCGGCTATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTTACAGCTCGACGATCTCTTTGCCCCGCGCGTGGCGAAGGCCCGTGTTGAAGGAAAAGTTTTGCGCTATGTTg                                                                                                                                                                                    >  2:617304/1‑251 (MQ=255)                                                                   ccGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATAtt                                                                                                                                                                            <  2:179693/251‑1 (MQ=255)                                                                   ccGAACCGGATCCGCGAGATGATCTTTCGGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGCCGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCGGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTGTTATGGCGAAGCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGTCTCGTGATGAAGGAAAAGGTTTGCGCTATGTTGGCAATAtt                                                                                                                                                                            <  2:350956/251‑1 (MQ=255)                                                                        ccGGATCCGCGCGATGATCTTTCTGGTATGGATGTAGCGCGTCAGCTATTGCTTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGGAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGCGATGTTGCCGCTTTTCTGGCGCATCTGTCACAGCTCGCCGAACTCTTTGCCCCGCGCGTGGCGAAGGACCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATtgatg                                                                                                                                                                        >  2:572124/1‑250 (MQ=255)                                                                                  gagaTGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATggtgcc                                                                                                                                                             >  1:35346/1‑247 (MQ=255)                                                                                   agaTGATCTTTCTGGTATGGATTTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATggtgcct                                                                                                                                                            <  1:517883/251‑6 (MQ=255)                                                                                   agaTGATCTTTCTGGTATGGATGTAGCGCGGAAGCTATTGATTCTCGCTCGTGAAACGGGACGGGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCGCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGGGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGTCt                                                                                                                                                            <  2:610390/251‑1 (MQ=255)                                                                                                ggTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGGGAACTGGAGCGGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGGTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGttaa                                                                                                                                                  <  2:493971/248‑4 (MQ=255)                                                                                                     tgatttagcgcagaagcTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTTAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGcc                                                                                                                                           <  1:601329/234‑1 (MQ=255)                                                                                                           agcgcGTAAGCTATTGAGTCTCGCTCGTGAAACGGGACGTTAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTgg                                                                                                                                    <  2:378185/250‑1 (MQ=255)                                                                                                              gcgTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGAAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATg                                                                                                                                 >  1:520194/1‑251 (MQ=255)                                                                                                              gcgTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCAGCTTTTATGGCGAATCTGTCACCGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGCAAAGTTTTGCGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGACGTGGATg                                                                                                                                 >  1:305033/1‑251 (MQ=255)                                                                                                               cgTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATgg                                                                                                                                >  2:243420/1‑251 (MQ=255)                                                                                                                  aagcTATTGCTTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGtg                                                                                                                                                                                                               >  1:506814/4‑169 (MQ=255)                                                                                                                  aagcTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGtg                                                                                                                                                                                                               <  2:506814/166‑1 (MQ=255)                                                                                                                  aagcTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTa                                                                                                                              >  1:483197/4‑250 (MQ=255)                                                                                                                  aagcTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGAGATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTTGATGGTa                                                                                                                              <  1:528503/247‑1 (MQ=255)                                                                                                                   agcTATTGATTCTCGCTCGTGAAACGGGACGTGAACGGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAAt                                                                                                                            <  1:545917/249‑1 (MQ=255)                                                                                                                             tctcGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTg                                                                                                                   <  1:129073/250‑1 (MQ=255)                                                                                                                                    cgTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCaaa                                                                                                             <  1:314684/249‑1 (MQ=255)                                                                                                                                                   gAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGaaa                                                                                             >  1:498391/1‑250 (MQ=255)                                                                                                                                                   gAACTGGAGCTGGCGGATAGTGAAATTGAACCTGTGCTGCCCGCAGAGTGTAACGCTGAGGGTGATGATGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTGTGCGCGATGTTGGCAATATTGATGAAGAGGGTGCCTGCCGCGTGAAGATTGCTGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGaaaa                                                                                            <  1:241683/251‑1 (MQ=255)                                                                                                                                                       gggAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGGTCAAAGTGAAAAATGGCGAAAACGcc                                                                                        <  2:393869/250‑1 (MQ=255)                                                                                                                                                        ggAGCTGGCGGATATTTAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGTGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGc                                                                                         <  1:105358/249‑1 (MQ=255)                                                                                                                                                         gAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCt                                                                                      <  1:546960/251‑1 (MQ=255)                                                                                                                                                           gCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGTTCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTgg                                                                                    <  2:438327/251‑1 (MQ=255)                                                                                                                                                             tGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGCTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAAAGATCCGCTGTTCAAAGGGAAAAATGGCGAAAACGCCCTGGcc                                                                                  >  2:246162/1‑251 (MQ=255)                                                                                                                                                              ggCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGAGCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATTGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGcc                                                                                  <  1:621111/250‑1 (MQ=255)                                                                                                                                                               gCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCtt                                                                                <  2:151005/251‑1 (MQ=255)                                                                                                                                                               gCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCtt                                                                                <  1:112968/251‑1 (MQ=255)                                                                                                                                                                 ggATATTGAAATTGAACCTGTGCTGCCAGCAGAGTTTAACGCTGAGGGTGATGTGGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTGTTGCGGTATGTAGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCtttt                                                                              <  2:215431/251‑3 (MQ=255)                                                                                                                                                                  gATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCtttta                                                                             <  1:94889/251‑4 (MQ=255)                                                                                                                                                                  gATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCtttt                                                                              >  1:604689/1‑248 (MQ=255)                                                                                                                                                                   atatTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTtat                                                                            >  2:494236/1‑251 (MQ=255)                                                                                                                                                                     attgaaattgaaCCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTtata                                                                           <  1:289231/250‑1 (MQ=255)                                                                                                                                                                         aaattgaaCCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCAc                                                                      <  1:284120/251‑1 (MQ=255)                                                                                                                                                                         aaattgaaCCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACACCCTGGCCTTTTATAGCCAc                                                                      >  2:574117/1‑251 (MQ=255)                                                                                                                                                                             tgaaCCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACtatt                                                                  >  1:175768/1‑251 (MQ=255)                                                                                                                                                                                  cTGTGCTGCCCGCAGAGTTGAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAg                                                                          <  2:499196/238‑1 (MQ=255)                                                                                                                                                                                  cTGTGATGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAg                                                                          >  1:499196/1‑238 (MQ=255)                                                                                                                                                                                   tgtgCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGc                                                            >  1:511183/1‑251 (MQ=255)                                                                                                                                                                                      gctgcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGc                                                         >  1:349653/1‑251 (MQ=255)                                                                                                                                                                                        tgcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTg                                                       >  2:506026/1‑251 (MQ=255)                                                                                                                                                                                           ccGCAGAGTTTAACGCTGAGGGTGATGTTTCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGGTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCg                                                    <  1:349463/251‑1 (MQ=255)                                                                                                                                                                                               agagTTTAACGCTGAGGGCGATGTTGCCGATTTTATGGCGAATCTGTCACATCGCGACGATCTCTTTGCCGCGCGCGTGGGGAAGGCCCGGGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGAGGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTAGCTGCCGCTGCCGTTgg                                                <  2:531451/251‑1 (MQ=255)                                                                                                                                                                                                   tttAACGCTGAGGGTGATGATGCCGCTTTTATGGCGAATCTGTCACAGCACGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGCAAAGTTTTGCGCTATGTCGGCCATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGACAACGCCATGGCCTTTTATAGACACTATTATCAGCCGCTGCCGTTggcgct                                            >  2:218248/1‑247 (MQ=255)                                                                                                                                                                                                      aaCGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAAGGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTgcg                                         <  1:67944/251‑1 (MQ=255)                                                                                                                                                                                                             aGGGTGATGTTGCCGATTTTATGGCGAATCTGTCAAAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGTCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGatat                                  >  2:246568/1‑251 (MQ=255)                                                                                                                                                                                                             aGGGTGATGTTCCCGCTTTTATGGCGAATCTGTCACAGCTCGCCGATCTCTTTGCCGCGCGCGTGGCGTAGGCCCGTGATGAAGGACACGTTTTGCGCTACGTTGGCCATATTGATGCCGATGGTGCCTGCCGCGTGAAGATTGCCGAGGGGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAAAGCCCTGGCCTTTTATAGCCACTATTATCCGCCGCTGCCGTTGGTGCTGCGCggcgtt                                  >  2:231539/1‑247 (MQ=255)                                                                                                                                                                                                               ggTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATgg                                <  1:134612/251‑1 (MQ=255)                                                                                                                                                                                                                   aTGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGCTGAAGATGGTGCCTGCCGCGTGAAGATTGACGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCCCTATTATCAGCCGCTGCCGTTGGTGCCGCGCGGATATGGAGCg                            >  2:204407/1‑251 (MQ=255)                                                                                                                                                                                                                          cGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATg                     <  1:107150/251‑1 (MQ=255)                                                                                                                                                                                                                           gCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGCAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGCAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATg                     >  1:29563/1‑250 (MQ=255)                                                                                                                                                                                                                             ttttATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGCTGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTCTTATCAGCCGCTGCCGTTGGGGCTGCGCGGATATGGTGCTGGCAATGACg                  >  2:168367/1‑251 (MQ=255)                                                                                                                                                                                                                              tttATGGCGAATCTGTCACAGCTCGACTATCTCTTTGCCGTGCGCGGGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTACCGGAGTGGATGGTAATGATCCGCTGTTCAAAGTGAATAATTGCGAAAACGCCCTGGCCGTTTATAGCGACTAGTAGAAGCCGCTGCTGGTGGTGCTGCGCGGGGATAGTGCTGGCAATGAc                   <  2:118569/249‑1 (MQ=255)                                                                                                                                                                                                                              tttATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGAGCTCTGTTGGCAATATTGTTGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGGGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATAGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCCACCGATGCCGTGGTTGCTGAGCGGCTATGGGGCgggacatgaagg                 >  2:502024/1‑242 (MQ=255)                                                                                                                                                                                                                                tATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGCAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGTCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGtt                <  2:29563/250‑1 (MQ=255)                                                                                                                                                                                                                                    gggAATCTGTCACAGCTCGACGACCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGTAAAAGTTTTGCGCTATGTTGGCAAGAGTGATGAAGATGGGGCCTGCCGCGTGAAGAGTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAGCGCCCTGGGATTGTGTAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGc           <  2:305033/249‑1 (MQ=255)                                                                                                                                                                                                                                      gAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTg         <  2:519942/251‑1 (MQ=255)                                                                                                                                                                                                                                             tCACAGCTCGACGATATCTTAGCCGCGCGCGTGGCGAAGGCCCGCGATGAAGGACAAGTTTTGCGCTATGTTGGCAATTTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAACCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTCCAGCTGCCGgtgt  >  2:62021/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      TCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGT  >  NC_000913/2267‑2753

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: ‑