Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,809	A→G	K252R (AAA→AGA)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,809	0	A	G	97.3%	263.4 / NA	73	K252R (AAA→AGA)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/41);  minor base T (1/1);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCG  >  NC_000913/498562‑499049                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  aaTGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAAt                                                                                                                                                                                                                                                >  2:471242/1‑251 (MQ=255) aaTGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAAt                                                                                                                                                                                                                                                <  2:345262/251‑1 (MQ=255)         ggCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCCCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCTCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGga                                                                                                                                                                                                                                         <  2:532127/250‑1 (MQ=255)             aaCAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCTGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGaaaa                                                                                                                                                                                                                                    <  2:320359/251‑1 (MQ=255)                aGCGTATGCGCATCTTCTGAAAAAAATGGCGAACCGGATCTGCTTCTGCTCTCTTATCATGGCATTCGCCAGCGTTCTGCAGATGAAGGCGATGATTACCCGCACCGTTGACGCACAACGACTCGTGAACTGGCTTCCGGCTTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGCCGAGCTTTGGTCGGGAACCCTGGCTGCTGCCTTATACCGACGAAACGATGATAATGCTCGGAGAAAAAgg                                                                                                                                                                                                                                 >  2:415512/1‑251 (MQ=255)                  cGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCg                                                                                                                                                                                                                               <  2:109492/251‑1 (MQ=255)                     aCGCGCTTCTTTTGCCAAACATGGCGAACCGGTGCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTGATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCGTCCGCATTGGAGATGGCACCGGAAAAAGTGATGATGACCTTTCAGGCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTATACCGACGAAACGCTGAGAATGCTAGGAGAAAAAGGCGTAg                                                                                                                                                                                                                            <  2:577119/251‑1 (MQ=255)                      cgcgCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAgg                                                                                                                                                                                                                           <  1:510870/251‑1 (MQ=255)                      cgcgCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCAATCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAgg                                                                                                                                                                                                                           >  1:499265/1‑251 (MQ=255)                         gcTTCTTTTGCCAAACCTGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCAAGCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGTCGTAGGTCa                                                                                                                                                                                                                        <  2:499265/251‑1 (MQ=255)                             cttTTGCCAAACACGGCGCACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAACAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATAtt                                                                                                                                                                                                                    >  2:452638/1‑251 (MQ=255)                             cgtttGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTAGCAGGGCAGTCCCCAGCGTTATGCAGAGGAAGGCGATGATTACCCGCAAGGTTGCCGCACAACGACTCGTGAGCTGGCTGCCGCATTGGGGATGGCAGCGGAAAAAGGGATGATGACCTTTCAGGCGCGCTTTGGTCGGGAACCCTGGCGGTTGCCTTATACCGGCGAAACGCGGAGAATGCTCGGACAAAAAGGCGTAGGTCatat                                                                                                                                                                                                                     <  2:109151/248‑1 (MQ=255)                               tttGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTc                                                                                                                                                                                                                   <  2:122824/250‑1 (MQ=255)                                ttGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCa                                                                                                                                                                                                                  <  2:332365/250‑1 (MQ=255)                                  gCCAAACATGGCGAACCGGCTCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGAGGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGATAATGCTCGGAGAAAAAGGCGTAGGTCTTATTCAGGt                                                                                                                                                                                                               <  2:630113/251‑1 (MQ=255)                                   ccAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTGTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTg                                                                                                                                                                                                              <  2:133461/251‑1 (MQ=255)                                     aaaCATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGAt                                                                                                                                                                                                            <  1:402198/251‑1 (MQ=255)                                                  cGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCttt                                                                                                                                                                                               <  1:642772/251‑1 (MQ=255)                                                  cGGATATGCTACTGATCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGCAGGAGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCATTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCttt                                                                                                                                                                                               >  2:271187/1‑251 (MQ=255)                                                            aCTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGAGATCATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGACGGCACCGGCAAACGTGATGATGACCTTTCCGTCGCGCTTTGGTCGGGAACCATGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGCAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGaat                                                                                                                                                                                     >  1:324329/1‑249 (MQ=255)                                                            aCTGCGCTCGTAGCATGGCATTCCCCAGCGTTATGCAGAGGAAGGCGATGAGTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCGTCCGCATTGGGGAGGGCGCCGGAAAAAGTGATGATGACCTTTGAGGCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCTGAGAAAAATGCGTAGGTCATATTCAGGAGATTTGCCCTGGCGTTGCTGCGGAtt                                                                                                                                                                                     <  2:310411/251‑1 (MQ=255)                                                                  ctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTgg                                                                                                                                                                               >  1:374287/1‑251 (MQ=255)                                                                      tATCATTGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGGCTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTAGACCGACGAAACGCTGAGAAGGCTCGGAGAAAAAGGCGTAGGTCATAGTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGgaga                                                                                                                                                                            <  1:656047/250‑1 (MQ=255)                                                                      tATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTGCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTGTGCTGCGGATTGTCTGgag                                                                                                                                                                             <  2:448945/249‑1 (MQ=255)                                                                            ggCATTCCCCAGCGTTATGCAGATGAAGGCGATGCTTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTggg                                                                                                                                                                     >  2:608877/1‑250 (MQ=255)                                                                               aTTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCCGCGGATTGTCTGGAGACGCTGGAaga                                                                                                                                                                  <  1:347755/251‑1 (MQ=255)                                                                                    ccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTg                                                                                                                                                             <  1:325598/251‑1 (MQ=255)                                                                                                aGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCGGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGATGGATGAGATTGCCGAGCAAAAcc                                                                                                                                                 >  1:324298/1‑251 (MQ=255)                                                                                                   tGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTg                                                                                                                                              >  1:252547/1‑251 (MQ=255)                                                                                                      aGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTAGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAgg                                                                                                                                           >  1:125342/1‑251 (MQ=255)                                                                                                              atTACCCGCAACGTTGCCGCACAACGTCTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGAGTCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAATGCGTAGGTCAGATTCAGGGGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTc                                                                                                                                   <  2:354819/251‑1 (MQ=255)                                                                                                              atTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTc                                                                                                                                   <  1:188110/251‑1 (MQ=255)                                                                                                              atTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTc                                                                                                                                   >  1:453759/1‑251 (MQ=255)                                                                                                                tACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGAGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCgg                                                                                                                                 >  1:236064/1‑251 (MQ=255)                                                                                                                  cccGCAACGTTGCCGCACAAAGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATCGCCGAGCAAAACCGTGAGGTCTTCCTCGGTg                                                                                                                               >  1:244099/1‑251 (MQ=255)                                                                                                                   ccGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGc                                                                                                                              >  1:211048/1‑251 (MQ=255)                                                                                                                       aaCGGTGCCGCACAACGACGCGTGAACTGGCTGCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGAGTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCCTCCTCGGTGCcggc                                                                                                                          <  2:238331/251‑1 (MQ=255)                                                                                                                              ccGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCggggaa                                                                                                                    >  1:125454/1‑247 (MQ=255)                                                                                                                                   caACGACTCGTGAACTGGCTTCCGCATTGGAGATGGCACCGGAAAAATTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAatat                                                                                                              <  1:667128/251‑1 (MQ=255)                                                                                                                                      cGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGaa                                                                                                           <  1:471242/251‑1 (MQ=255)                                                                                                                                           cGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTGCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATAGGAatatat                                                                                                      <  2:475989/251‑1 (MQ=255)                                                                                                                                            gTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGCTGATGAACTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATAtt                                                                                                     >  1:642861/1‑251 (MQ=255)                                                                                                                                                   ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcgc                                                                                              <  1:195158/251‑1 (MQ=255)                                                                                                                                                   ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgc                                                                                                <  1:271187/249‑1 (MQ=255)                                                                                                                                                        ccGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAt                                                                                         <  1:28454/251‑1 (MQ=255)                                                                                                                                                        ccGCAGTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCAGATGCCTTATACCGACGAAACGCGGAGAATGCTCGGAGAAAAAGGCGTGGGTCATATTCAGGTGAGGTTCCCTGGCGTTGCTGCGGATTGTCTGGAGACGCTGGTAGAGATTGCCGAGCAAAACCTTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATAGTCCGGCGCTTAAt                                                                                         <  2:352756/251‑1 (MQ=255)                                                                                                                                                              tGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCCGCGCTTAATGCCACg                                                                                   >  1:281234/1‑251 (MQ=255)                                                                                                                                                                gggATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGcc                                                                                 <  2:281234/251‑1 (MQ=255)                                                                                                                                                                   aTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATCTATTCCGGCGCTTAATACCACGCCggg                                                                              >  1:167300/1‑250 (MQ=255)                                                                                                                                                                   aTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCGGAGCCTAATGCCACGCCggg                                                                              >  2:448888/1‑250 (MQ=255)                                                                                                                                                                            gAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGCCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGa                                                                    <  1:120133/251‑1 (MQ=255)                                                                                                                                                                             aaaaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGaa                                                                   <  2:6783/251‑1 (MQ=255)                                                                                                                                                                               aaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAat                                                                 <  1:584602/251‑1 (MQ=255)                                                                                                                                                                                   ggatgatgaCCTTTCAGTCGCGCGTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAatgatg                                                             <  1:325373/250‑1 (MQ=255)                                                                                                                                                                                       gatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCt                                                          <  1:485771/250‑1 (MQ=255)                                                                                                                                                                                            ccTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCtt                                                    <  1:29627/251‑1 (MQ=255)                                                                                                                                                                                               ttCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAAAAATTCGGCGCTTAATTCCTCCCCGGCAAC‑TATCGAAAAGATGGCTAATCttgtt                                                 >  2:466464/1‑251 (MQ=255)                                                                                                                                                                                                   gTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTAATGCCACGCCGG‑GACATATCGAAATGATGGGTAATCTTGTTTCcg                                             >  1:555020/1‑251 (MQ=255)                                                                                                                                                                                                    gagcgcTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCcgc                                            <  2:125342/249‑1 (MQ=255)                                                                                                                                                                                                        gcgttGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCGCGGAGAAAAAGGCTTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCGGGAAGAGATTGCCGAGCAAAGCCGTGAGGGCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACTCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCCGCGTAt                                        <  1:415512/248‑1 (MQ=255)                                                                                                                                                                                                             ggTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAAAGATGGCTAATCTTGTTGCCGCGTATCGCt                                    >  1:16523/1‑250 (MQ=255)                                                                                                                                                                                                             ggTCGGGAACCCTGGCTGATGCCTTATAACGACGAAACGCCGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTTATGCCCCGCCGG‑GACATATCGAAATGATGGCTAATTTTGTTGCCGCGTATCGCt                                    >  1:220766/1‑250 (MQ=255)                                                                                                                                                                                                                  ggAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAACGGCGTAGGCCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGCATATATTTCGGGGCGTTATGGCCCGCCGG‑AACCTATTGAAATGATTGATAATCTTTTTTCcgcggtttcccaaagcc                              >  2:600411/1‑238 (MQ=255)                                                                                                                                                                                                                   gAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTg                             <  2:468329/251‑1 (MQ=255)                                                                                                                                                                                                                      cccTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGGGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAAAATATTACGGCGCTTAATGCCACGCCGG‑AATATATCGACATGATGGCTCAACTTTTTTGCGCGTATCGCTCAAGCTGAgg                          >  1:588778/1‑250 (MQ=255)                                                                                                                                                                                                                          ggCTGATGCCTTATACCGACGAAACGCTGCGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcgg                                                                                                                           >  2:630655/1‑151 (MQ=255)                                                                                                                                                                                                                          ggCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGTCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCcgg                                                                                                                           <  1:630655/151‑1 (MQ=255)                                                                                                                                                                                                                              gATGCCTTAGACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGGGCGGTAAAg                   <  2:225021/250‑1 (MQ=255)                                                                                                                                                                                                                                   cTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGGCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTg              >  1:346875/1‑250 (MQ=255)                                                                                                                                                                                                                                   cTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGAAGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTgg             >  1:58559/1‑250 (MQ=255)                                                                                                                                                                                                                                        aCCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGcc        <  2:549676/251‑1 (MQ=255)                                                                                                                                                                                                                                           gacgaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAAAATATTTCGGCGCTTAATGGCACGCCGG‑AAACTATCGAAATGATGGGTAATTTTGTTTCCGCGTGTCGCTAAAACTGAGCGGGAAAGAAATGAGCGCCGTc     >  1:221156/1‑251 (MQ=255)                                                                                                                                                                                                                                              gaAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTTTTGCCGCGTATCGCTAAAGCTGAGCGGTGAAGAACAGAGCGCCGTcgcg  >  2:163240/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  AATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGG‑AACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCG  >  NC_000913/498562‑499049

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