Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,116	C→T	P21L (CCT→CTT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,116	0	C	T	98.3%	188.6 / ‑4.6	60	coding (62/963 nt)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/24);  minor base . (0/1);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACG  >  NC_000913/497871‑498361                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gTGAGGGGAACAGGCCCGCACAAGCAAACTTAACAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑g                                                                                                                                                                                                                                                  <  1:482930/251‑1 (MQ=255)           cAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTGTTCCGCTACAATTAGCAACAAGTTGAATCGAGAAGAGGCGGTAATGCGTCAGACTAAATCCGGTATCCTGCTGGCAGACCGGTGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACg                                                                                                                                                                                                                                        <  1:352844/251‑1 (MQ=255)             ggCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCt                                                                                                                                                                                                                                      <  1:109467/251‑1 (MQ=255)                     caAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTCAAACGCTATCTGaa                                                                                                                                                                                                                               >  1:329884/1‑250 (MQ=255)                          aaaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAAtt                                                                                                                                                                                                                         >  2:454951/1‑251 (MQ=255)                           aaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAAttt                                                                                                                                                                                                                        >  1:132018/1‑251 (MQ=255)                            aCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAAtttt                                                                                                                                                                                                                       <  1:18788/251‑1 (MQ=255)                                 tCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGc                                                                                                                                                                                                                  >  1:277364/1‑251 (MQ=255)                                  cagaaaTCTCAGGCCGGATATTCATTCGGCCTTTTACAACAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGACTCGATAAGAGGCGGTACTGCGTTAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGGTGACCCCACACTTGGA‑GCGGTCAAAAGATATCTGAAACAATTTTTAAGCg                                                                                                                                                                                                                 >  1:368333/5‑251 (MQ=255)                                   aGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCg                                                                                                                                                                                                                 >  2:434889/1‑250 (MQ=255)                                       aTCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCATCACTTTTCCGATACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAACACCGGTATCCTGCTGGCAAACCTGGGTACCCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGCCAgc                                                                                                                                                                                                            >  1:577560/1‑250 (MQ=255)                                        tctcAGGCCGGATATTCATTCGGCCTTTTACAAAAACATCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAc                                                                                                                                                                                                           >  1:306599/1‑251 (MQ=255)                                           cagtccGGATATTCATTCGGCCTGTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTTATGTGACCTCTCACAGCAATTAGTTCGTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAc                                                                                                                                                                                                           <  2:428266/244‑1 (MQ=255)                                                    tatTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGAt                                                                                                                                                                                                >  1:462506/1‑250 (MQ=255)                                                        cattcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTc                                                                                                                                                                                           >  1:64165/1‑251 (MQ=255)                                                        cattcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCAAAGCCATTAGTTCTTCTTCATCACTTTTCCGCTACAATTCTCCACAAGTTGAATCGATAAGAGGCGGTCATGCGTCAGACTAAAACCGGTCTCCTGCTGGCAACCCTGGGCACGCCCGATGCCCCCCCACTTGAA‑GCGGTAACACGCTATCTGAAACAACTTTTAAGCGACAGACGCGTGGTTGATACCTc                                                                                                                                                                                           >  1:375894/1‑251 (MQ=255)                                                         attcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTcc                                                                                                                                                                                          >  1:108101/1‑250 (MQ=255)                                                          ttcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTAACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGAGTGGTTGATACCTCAc                                                                                                                                                                                         >  1:407704/1‑251 (MQ=255)                                                                   ttACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGGTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCTGTAAAACGCTATCTGAAACAATGTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTa                                                                                                                                                                                 <  1:608870/250‑1 (MQ=255)                                                                            ccTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCCAt                                                                                                                                                                       >  2:8130/1‑251 (MQ=255)                                                                            ccTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCAt                                                                                                                                                                       >  1:220051/1‑251 (MQ=255)                                                                              tCAATCTGCGCTATTCATTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCCCTTTTCCGCTACAATTCTCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCAGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtg                                                                                                                                                                     >  1:432798/1‑251 (MQ=255)                                                                                      cgcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCgcgg                                                                                                                                                             >  1:218160/1‑251 (MQ=255)                                                                                       gcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCgcggc                                                                                                                                                            >  1:176387/1‑251 (MQ=255)                                                                                                ttttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAtttt                                                                                                                                                    <  2:427778/250‑1 (MQ=255)                                                                                                  ttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGt                                                                                                                                                                                                       <  1:336388/197‑1 (MQ=255)                                                                                                  ttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGAACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGCCAGACGCGt                                                                                                                                                                                                       >  2:336388/1‑197 (MQ=255)                                                                                                      gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTg                                                                                                                                             <  1:148775/251‑1 (MQ=255)                                                                                                                        aaTTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCg                                                                                                                           <  1:348/251‑1 (MQ=255)                                                                                                                            aGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGc                                                                                                                       >  2:199003/1‑251 (MQ=255)                                                                                                                                       cTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAACCCGGGATACTGCTGGCCAAACTGGGGACGCCCGCTGACCCCACACTTGAA‑GCGGTAAAACGATGTATGAACCAATTTTTAAGCGAAAGAAGCGCGGTTGTGACCTCACGGTTGTTTTGGTGGCCATTGCCGCGGGGTGTGGTTTTTCCGGTGGGCTCGCCGCCGGTGGCGCAGCGGTATGcttcg                                                                                                            >  2:1352/1‑247 (MQ=255)                                                                                                                                         cACTTTTCCGCTGCAATTCTCAACAAGTTGAATCGATAAGAGGCGGTACTGCGTCACACTCAAACCGGTATTCTGCTGGCAAACCTGGGTAAGCCCGATGCCACCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGTCATTGCTGTGCGGCGTGATTTTGCCGGTGCGCGCGCTGGGTGTGGCGCGGCTGTGTGTCtctgt                                                                                                          >  2:226755/1‑251 (MQ=255)                                                                                                                                                     aCAATTATCAACCAGTTGAATCGATAAGAGGCGGTAATGCGTCAGAATAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGCCAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGGCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTGTGGATGGGAgg                                                                                              >  2:575184/1‑251 (MQ=255)                                                                                                                                                       aaTTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCTACAGACGCGTGGTTGATACCTCGCGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAGTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTCTCTGGATGGAAGGTg                                                                                            <  2:407704/251‑1 (MQ=255)                                                                                                                                                       aaTTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                            <  2:64165/251‑1 (MQ=255)                                                                                                                                                           aTCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATCCCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTc                                                                                        <  2:449592/251‑1 (MQ=255)                                                                                                                                                              aacaaGATGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCCGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgcc                                                                                     >  1:74309/1‑251 (MQ=255)                                                                                                                                                                    ttGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGAt                                                                               >  1:233447/1‑251 (MQ=255)                                                                                                                                                                     tGAATCGATAAGAGGCGTTAATGCTTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTGAAGCGACAGACGCGTGGTTGATACCCCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGACTTTGCCGCTGCGCTCGCCGCGTGTGGCTAAGCTGTATGCCTCTGCCTGGATTGAAAGTGGTTCGCAGCTGAt                                                                               <  2:607447/250‑1 (MQ=255)                                                                                                                                                                     tGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGAt                                                                               <  1:454951/250‑1 (MQ=255)                                                                                                                                                                      gAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATgg                                                                             <  2:176387/251‑1 (MQ=255)                                                                                                                                                                          cGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATGTTTAAGCGACAGACGCGTGGTTGATCCCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATGTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGttt                                                                          <  2:329884/250‑1 (MQ=255)                                                                                                                                                                                  ggCGGTAATGCGTCAGACTAAAACCGGTTTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgcc                                                                 <  1:280713/251‑1 (MQ=255)                                                                                                                                                                                     ggTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGc                                                              <  2:462506/251‑1 (MQ=255)                                                                                                                                                                                        aaTGCGTCAGACTAAACCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAc                                                           >  1:225536/1‑251 (MQ=255)                                                                                                                                                                                             gTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCACcagc                                                        >  2:486273/1‑249 (MQ=255)                                                                                                                                                                                               cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                                    >  1:41955/1‑251 (MQ=255)                                                                                                                                                                                                 gACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgcgc                                                  <  1:373947/251‑1 (MQ=255)                                                                                                                                                                                                   atAAAACCGGATTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAGACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGATGTTATGGTGACCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGATGCTGTATGCCTCGGGCTGGAGGGAAGGTGGGTCGCCGCTGACGGTTTACAGCCGACTGTCACAGCAGgcgg                                                  >  2:132439/2‑248 (MQ=255)                                                                                                                                                                                                      aaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGCTGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGCTACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCCGCAACAGCAGGCGCTGGcc                                             >  2:229753/1‑250 (MQ=255)                                                                                                                                                                                                           ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACg                                        <  2:63940/251‑1 (MQ=255)                                                                                                                                                                                                                     tgGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACCGCAGGCGCTGGCACACCGTTTACCgg                                >  2:340052/1‑249 (MQ=255)                                                                                                                                                                                                                     aggCAAACCTGGGTACGCCAGATGCCCCCCCAC‑TGAAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGgag                              <  2:238725/250‑1 (MQ=255)                                                                                                                                                                                                                      ggCAAACCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTGTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTTTATGCCTCTGTCTGGATTGAATGTGGCTCGCCGCTGATGGTTGACATCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGgag                              <  1:304704/250‑1 (MQ=255)                                                                                                                                                                                                                         aaaCCTGGGTACGCCCGATGGCCCCACTCTTGAA‑GAGGCACAACGCTATCTGGAACAATTTTTAAGCGACAGCCGAGTGGTTGATACATCACGGTTGTTATGGTGGCCATTGCTGCGCGGAGTGATTTTGCCGCTGCGCTCGACGCGTGTGGCGCAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGAAGGAGCTGGCACAACGTTTACCGGAGATGc                          >  2:73477/1‑251 (MQ=255)                                                                                                                                                                                                                         aaaCCTGGGTACGCCCGATGCCCCCACACTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTACGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGTGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGTAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCCCACCGTTTACCGGAGATGc                          >  2:375647/1‑251 (MQ=255)                                                                                                                                                                                                                           aCCTGGATAAGCTCGAAGCCCCCACACTTGAA‑GCGGTAAAACACACTCTGAAACAATTTTTAAGAGACAGACGCCTGGTTGATACCTCACCGTTGTTATGGTGGCCATTGCTGCGCGGCGTGCTTTTGCCGCTGCGCTCGCAGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGCTGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGCGATGccc                        >  2:572196/1‑251 (MQ=255)                                                                                                                                                                                                                              tGGGTCCGCCCGATGCCCCCACACTTGAC‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGCCGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGACGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGt                      >  2:304235/1‑250 (MQ=255)                                                                                                                                                                                                                                     gCCCGATGCCCCAACCCTTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGGTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTGTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTgg              <  2:577560/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cacaCTTGAA‑GCGGTAAAACGCTATCTGAACCAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAGCAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACg  <  1:399950/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAA‑GCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACG  >  NC_000913/497871‑498361

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