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.2%	175.5 / NA	55	P21L (CCT→CTT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/22);  minor base G (0/1);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTAC  >  NC_000913/497883‑498360                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      ggCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCt                                                                                                                                                                                                                                      >  1:50244/1‑251 (MQ=255)          aaGCAAACTTATCCGCAATCCCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCCATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGAATAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCTCCACACTTGAAGCGGTaaaaagagtgctgaaa                                                                                                                                                                                                                              >  2:242620/1‑238 (MQ=255)           aGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACa                                                                                                                                                                                                                            <  1:338574/251‑1 (MQ=255)                        gCAATCTCAGGCCGGATATTCAGTCGGCCGTTTACAAAAACCTCAATCTGCGCTATTCCTTGTTCTGATTTGACCTCTCACAGCAATTAGTTCGTCTTCCTCACTGTTCCGCTACAATTAGCAACAAGTTGAATCGACAAGAGGCGGTAATGCGTCGTACTAAAACCTTTATCCTGCTGGCAAGCCTGGGTTCGCCCGATGGCCCCACACTTGAAGCGGTAAAACGCTGTCTGAAGCAA‑GTTTCAAGCGAc                                                                                                                                                                                                              <  2:182054/251‑1 (MQ=255)                          aagctcAGGCCGGATATTCATTCGGGCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGCCCTCTCCCAGCAATTAGTTCTTCTTCCTCACTGTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAg                                                                                                                                                                                                            <  2:254773/248‑1 (MQ=255)                            tctcAGGCCGCATATTCATTCGGCCTTTTACAAAACCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAAGGCCTCAGACTAAAACCGGTATACTGCTGGCAAACCTGGGAACGCCCGATGCCCCAACACTTGAAGCGGTAAAACGCTAGCTGAAACAA‑TTTTTAAGCGACAGAc                                                                                                                                                                                                          >  2:316300/1‑251 (MQ=255)                                   ccGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTCAGCGACAGACGCGTGGt                                                                                                                                                                                                   >  1:38321/1‑251 (MQ=255)                                           tcattcGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCt                                                                                                                                                                                           >  1:95318/1‑251 (MQ=255)                                                  gCCTTCTACAAAAACCTCAATCTGCGCTGTTCCTTTTTCTGCTTTGACATCTAACAGCAATTAGTTCTAATTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATACGAGGCGGTAATGCGACAGAATAAAACCGGTATCCTGCTCGCAAACCTGGGGACGCCCGACGCCCCCCCACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGAAGCGTGGTTGATACCTCACGgtt                                                                                                                                                                                    >  2:456434/1‑251 (MQ=255)                                                       ttACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATATGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGAAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGGGGTTGCTACCTCACGGTTGTTAt                                                                                                                                                                               >  1:200461/1‑251 (MQ=255)                                                           aaaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAAAAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTTTTAtggtg                                                                                                                                                                           >  1:385903/1‑251 (MQ=255)                                                           aaaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGGTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAATCAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggtg                                                                                                                                                                           <  2:385903/251‑1 (MQ=255)                                                              aaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACCAGTTGAATCGATAAGAGGCGGTAATGAGTCAGACTCAAACCGGTATCCTGATGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGTTACCTCACGGTTGTTATGGTGGcc                                                                                                                                                                        >  2:172146/1‑251 (MQ=255)                                                               aCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAAATTTTTAAGCGACAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGcc                                                                                                                                                                        >  1:99231/1‑251 (MQ=255)                                                                   cAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCAtt                                                                                                                                                                     <  1:178272/249‑1 (MQ=255)                                                                         gcgcTATTCCTTTTTCTGATTTGACCTCTCTCAGCAATTTGTTCTTCTTCCTCACTTTTCCGCTTCAATTATCAACAAGTTGAATCGATAAGAGGCGGTAGTGCGTCAGACTAAACCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATTCCCCCACACGTAAAGCTGTAATGCGCTATCTGAAACAA‑TTTTTAAGCGACAGTCGCGTGGTTGATACCTCACTGTTGGTATGGTGGCCATTGCTgcg                                                                                                                                                               <  2:268880/249‑1 (MQ=255)                                                                         gcgcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGGGGCCATTGCTGcgcg                                                                                                                                                             >  2:104931/1‑251 (MQ=255)                                                                             tATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGGTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCgcggcgg                                                                                                                                                         >  2:65405/1‑250 (MQ=255)                                                                               ttCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACCGACGCGTGGTTGATACCTCACGGTTGTTTTGGTGGCCCTTGCTGCGCGGCGTgt                                                                                                                                                       >  1:47185/1‑250 (MQ=255)                                                                                    ttttCTGATTTGACCTCTCTCAGCAATTAGTTCTTCTTCCTCACTTGTCCGCTACAATTAACAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAATCGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGGTGTTATGGTGGCCATTGCTGCGCGGCGTGAtttt                                                                                                                                                   <  2:256378/250‑1 (MQ=255)                                                                                      ttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGc                                                                                                                                                 <  1:327203/250‑1 (MQ=255)                                                                                          gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTg                                                                                                                                            >  1:358360/1‑251 (MQ=255)                                                                                           aTTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgc                                                                                                                                           <  1:296009/251‑1 (MQ=255)                                                                                           aTTTGACCTCTAACAGCCATTAGCTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATACGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTCCTGGAAAACCTGGGTACGCCCGATGCCCCCAAACTTGAAGCGGATAAACGCTATCTGAAAAAA‑TTTTTAAGCGACAGAAGCGTGGTTGATCCCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTTCCGcctg                                                                                                                                           >  2:222994/1‑248 (MQ=255)                                                                                               gACCTTTCACAGGAATTAGATCTTTTTCCTCACTTTTCCGCTACAATTATCACCAAGTTGACTCGATAAGAGGCGGTAATGCGTAAGACTAAAACCGGTATCCTGCTGGCAAACCTGTGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGCGAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCCATTGCTGTCCGGCGGGTTTTTGCCGccgatccc                                                                                                                                       >  2:151846/1‑244 (MQ=255)                                                                                               gACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTc                                                                                                                                                                                          >  1:433063/1‑200 (MQ=255)                                                                                               gACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAAGCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTc                                                                                                                                                                                          <  2:433063/200‑1 (MQ=255)                                                                                                    ctcACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAACCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgcc                                                                                                                                    <  2:102571/249‑1 (MQ=255)                                                                                                          gCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTACTGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGTTTTTGCCGCTGCGCTCGCCGCGTGTgg                                                                                                                            >  2:27571/1‑251 (MQ=255)                                                                                                            aaTTAGTTCTTCTTCCTCACTTTTACGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCATACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGGAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACACGTGGTTGATACCTCACGGTCGTTTTGGTGGCCCTTGCTGAGCGGCGGGATTTTGCCGCTGCGCTCGCCGCGGGTgggg                                                                                                                          >  2:381245/1‑249 (MQ=255)                                                                                                                aGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGc                                                                                                                      >  1:160214/1‑251 (MQ=255)                                                                                                                           cTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATACGAGGCGGGACTGCGTCAGAATCACAACGGTATCCTTCAGGCAAACCTGGGTACGCCCGATGCCCCCACCATTGAAGCGGGCAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTTGTTGATACCTCACGGTTGTTATGGCGGCCCTTGGTGCGCGGCGGGTATTTGCCTCTGCGCTCGCCGCGTTTGGCGTAGCttcatgtctct                                                                                                           >  2:432129/1‑241 (MQ=255)                                                                                                                              ccgtttCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCTCACTGGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACTGTTGTTATGGTTGCCATTGCTGCGCGGTGTGATTTTGCCGCTGCGCTCGCCGCGTTTGGCGAATCTGTATGCCtctgtc                                                                                                        <  2:181136/248‑1 (MQ=255)                                                                                                                                    ccGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACACGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCtggatg                                                                                                  >  1:252048/1‑251 (MQ=255)                                                                                                                                      gCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCtggatgg                                                                                                 >  1:218438/1‑250 (MQ=255)                                                                                                                                              gaTCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCGATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAGGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTTTCTGGATGGAAGGTGGc                                                                                         <  1:321549/249‑1 (MQ=255)                                                                                                                                                     aaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCt                                                                                 >  1:164245/1‑251 (MQ=255)                                                                                                                                                         tGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGTTAc                                                                                                                                                                                             >  2:5107/1‑139 (MQ=255)                                                                                                                                                         tGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATAc                                                                                                                                                                                             <  1:5107/139‑1 (MQ=255)                                                                                                                                                         tGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGAt                                                                              <  1:212632/250‑1 (MQ=255)                                                                                                                                                                aTAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTCTAAGCGACAGACGCGTGGTTGATACCTCACGTGTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTTCCGCGGGGCTCGCCGCGTGTGGTTAAGCTGGATGCCGCGGTCTGGATGGAAGGTGGCTCGCCGCTGATTGTTTACa                                                                      <  1:456434/251‑1 (MQ=255)                                                                                                                                                                      ggCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCTGATGGTTTACAgccgcc                                                                >  2:336049/1‑251 (MQ=255)                                                                                                                                                                          gTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCa                                                            <  2:160214/251‑1 (MQ=255)                                                                                                                                                                          gTAATGCGTCAGAATAAACCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACCCTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTCAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGAGGCTCGCCGCTGATGGTTTACAGCCGCCAGc                                                             >  1:303584/1‑250 (MQ=255)                                                                                                                                                                                   cAGACTAAAACCGGTATCCTGCCGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCGCACGGTTGTTATGGTGGCCATTGCGGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                                   <  2:50244/251‑1 (MQ=255)                                                                                                                                                                                       cTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTg                                               >  1:196452/1‑251 (MQ=255)                                                                                                                                                                                           aaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcac                                           <  1:184352/251‑1 (MQ=255)                                                                                                                                                                                                  aTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGtt                                     <  1:7471/250‑1 (MQ=255)                                                                                                                                                                                                   tCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGttt                                    <  2:252048/250‑1 (MQ=255)                                                                                                                                                                                                             aaaCCTGGGTACGCACGAAGCACCCACACTTGAAGCGGTAAAACGCTATCCGAAACAC‑TTTTTAAGCGACAGACGCGTGGTTGCTACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGCAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGc                         >  2:76291/1‑251 (MQ=255)                                                                                                                                                                                                               aCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGccc                       >  2:343571/1‑251 (MQ=255)                                                                                                                                                                                                                       aCGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAgcgc                <  1:343571/250‑1 (MQ=255)                                                                                                                                                                                                                       aCGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCt               >  2:97881/1‑251 (MQ=255)                                                                                                                                                                                                                                gCCCCCACACTTGAAGCGGTACACCGCTATCTGAAACCC‑TTTTTAACCGACAGACGCGTGGTAGACACCCCCCGCTTGTTATGGTGCACATTGCCACGAGGCGTGATTTTGCCGCTGCCCTCGCCGCGTGTGGCGAAGCTGTCTGCCTCTGACTGGGTGGCAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCAGGCACAACGTTTACCGGCGATGCCAGTAGCGCTGGGAATg        >  1:194408/1‑249 (MQ=255)                                                                                                                                                                                                                                    ccACACTTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCTCCAGCAACAGCAGGCGCTGTCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTAc  <  1:268844/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      GGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTGAAACAA‑TTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTAC  >  NC_000913/497883‑498360

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