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	98.2%	199.8 / NA	54	K252R (AAA→AGA)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (33/20);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGA  >  NC_000913/498570‑499050                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  ggCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCgg                                                                                                                                                                                                                                          <  2:182624/249‑1 (MQ=255)        aGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAgg                                                                                                                                                                                                                                 >  1:14191/1‑251 (MQ=255)                    tcttTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAGCCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTCGGTCatat                                                                                                                                                                                                                     >  2:481653/1‑251 (MQ=255)                         tGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAgg                                                                                                                                                                                                                >  2:252805/1‑251 (MQ=255)                                cATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGAtgtg                                                                                                                                                                                                         >  2:128454/1‑251 (MQ=255)                                cATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGAtgtg                                                                                                                                                                                                         <  1:328586/251‑1 (MQ=255)                                      gAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGGTGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGACTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGCTGTGCCCggg                                                                                                                                                                                                   >  2:3757/1‑251 (MQ=255)                                      gAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCggg                                                                                                                                                                                                   >  1:17508/1‑251 (MQ=255)                                       aaCCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCAAGCGTTATGCAGATGAAGGAGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGCCGAAACGCTGAGAATGCTCGGCGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGc                                                                                                                                                                                                  >  2:383433/1‑251 (MQ=255)                                         ccGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑GCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCTGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGTATAAAAATGCGTAGGTCATGTGCAGGTGATGTGCCCGGGCtt                                                                                                                                                                                                <  2:537753/251‑1 (MQ=255)                                            gATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTCCCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGCGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTc                                                                                                                                                                                                                         >  1:503894/1‑223 (MQ=255)                                            gATCTGCTACTGCTCGCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGAGGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTGGGTCGGGAACCCTGGCGGATGCCGTATACCTACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTc                                                                                                                                                                                                                         <  2:503894/223‑1 (MQ=255)                                                        ctctctTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAACGGCGTAGGTCATATTCAGGCGATGTGCCCGGGCTTTGCTGCGGATTGTCt                                                                                                                                                                                 >  1:567420/1‑251 (MQ=255)                                                                 cATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCt                                                                                                                                                                        <  1:323562/251‑1 (MQ=255)                                                                   gggCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGAGGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCAGATGCAGGTGAGGTGCCCGGGCGTGGCGGCTGATTGTCTGGAGACGCTgg                                                                                                                                                                      <  2:503626/250‑1 (MQ=255)                                                                    ggCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTgg                                                                                                                                                                      >  2:74738/1‑250 (MQ=255)                                                                         tCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGGCCTTTCAGTCGCGCTGTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCGGGAagag                                                                                                                                                                 <  2:555547/250‑1 (MQ=255)                                                                           cccAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGCGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAtt                                                                                                                                                              >  1:102932/1‑251 (MQ=255)                                                                             cAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGGCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGc                                                                                                                                                            <  1:537583/251‑1 (MQ=255)                                                                                cGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGACGAGCTTGCCgg                                                                                                                                                         >  2:173498/1‑250 (MQ=255)                                                                                                gCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTc                                                                                                                                         <  2:136553/251‑1 (MQ=255)                                                                                                gCGATGATTACCCGCAACGTTGCCGCA‑AACCACATCGTGCACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGCTGACCTTTCAGTCGCGCTTTGGTCGGGCACCCTGGCTGATGCCTTATACCGACGCACCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAATGACGAGCAAAAACGTGCGGTc                                                                                                                                         >  2:332636/1‑251 (MQ=255)                                                                                                 cGATGATTACACGCAACGTTGCCGCACCACGAC‑TCGCGAACGGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCCGTCGCGCTTTGGTCGGGCACCATGGCTGATGCCTTATACCGACGACACGCTGAGACTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGCGATGTGCCCGGGCTTTGCTGCGGCTTTTCTGGAGACGCTGGGCGAGATTGCCGAGCAAAAccctgcgggct                                                                                                                                        >  2:427858/1‑242 (MQ=255)                                                                                                  gatgatTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGTTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCtt                                                                                                                                       >  1:493443/1‑251 (MQ=255)                                                                                                     gatTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCt                                                                                                                                    >  2:115255/1‑251 (MQ=255)                                                                                                         aCCCGCAACGTTGCCGCACAACGAC‑TCGTTAACTGGCTTCCGCATTGGGGATGTCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTGTGCTGCGGATTGTCTGTAGTCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGt                                                                                                                                <  1:279280/251‑1 (MQ=255)                                                                                                          cccGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGc                                                                                                                                                                         <  1:310434/209‑1 (MQ=255)                                                                                                                  gTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCgggggg                                                                                                                       >  2:182965/1‑251 (MQ=255)                                                                                                                      ccGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCGCGGTGCCGGCGGGaaaa                                                                                                                   <  1:450205/251‑1 (MQ=255)                                                                                                                        gCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAACGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTc                                                                                                                                                                                                                                            >  2:169336/1‑128 (MQ=255)                                                                                                                        gCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaaa                                                                                                                  >  2:565743/1‑250 (MQ=255)                                                                                                                        gCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTATACCGACGAAACGCTGAGAATGCTc                                                                                                                                                                                                                                            <  1:169336/128‑1 (MQ=255)                                                                                                                            aaCGAC‑TCGTGAACTGGCTTCGGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAAACCTGGCTGATGACTTATAACGACGAAACGCTGAGAATGCTAGGGGAAAAAGACGTCGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGCGACGCTGGCAGAGATTGCCGAGCACGACCGTGAGGACATCCTCGGTGCCGGCGGGAAAAAAGATg                                                                                                             >  1:234212/1‑251 (MQ=255)                                                                                                                                            ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGCCGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATTTTCCGgcgc                                                                                              >  1:459754/1‑251 (MQ=255)                                                                                                                                            ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                                                                                               <  2:110027/250‑1 (MQ=255)                                                                                                                                               ttCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCttt                                                                                           >  1:395502/1‑250 (MQ=255)                                                                                                                                                  cGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTTCGGGGCTTTATg                                                                                        >  2:261975/1‑251 (MQ=255)                                                                                                                                                       tGGGGATGGCACCGGAAAAAGTGGTGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCAc                                                                                    <  1:3757/250‑1 (MQ=255)                                                                                                                                                        ggggATGGCACCGGAAAACGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCGCTTAATGCCACGc                                                                                  >  1:166152/1‑251 (MQ=255)                                                                                                                                                                 cccGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACata                                                                         <  2:12240/250‑1 (MQ=255)                                                                                                                                                                 aCCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTAATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGCAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATATATTCCGGCGCTTAAAGCCAAGCCgggaaat                                                                          >  1:276071/1‑245 (MQ=255)                                                                                                                                                                       aaaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCtt                                                                                            <  1:208488/226‑1 (MQ=255)                                                                                                                                                                       aaaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGaaa                                                                   <  1:128454/251‑1 (MQ=255)                                                                                                                                                                        aaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTTCCGCGCTTTATTCCACCCCgggaactaactgaa                                                                   >  2:138551/1‑238 (MQ=255)                                                                                                                                                                        aaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTTCGGCGCTTTATGCCACCCCGGAACatattgcaa                                                                   >  2:113799/1‑245 (MQ=255)                                                                                                                                                                                    aCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATAAGAATATATTACGGGGCTTAAATCCACGCCGGGACCCAACGAAATGATGGCTaaact                                                      >  2:25044/1‑248 (MQ=255)                                                                                                                                                                                     ccTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCtt                                                     >  1:483815/1‑251 (MQ=255)                                                                                                                                                                                     ccTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATGGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCTGGAGGTCGTCCTCGGTGCCGGCGGGAAAAAAGATGAATTTATTCCGGCGCGTTATGCCACGCCGGAACATATCGAAATTATGTCTAATCtt                                                     <  2:71463/251‑1 (MQ=255)                                                                                                                                                                                        ttCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgt                                                   <  2:473199/250‑1 (MQ=255)                                                                                                                                                                                            gTCGCGCTTTGGTCGGGAACCATGGCTGATGCCTTATACCGACGACACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATTTTCACGCGCTTAATGCCACGCCGGAACAAATCGAAATTATGGCTAATCTTTTTGCcg                                              >  1:48612/1‑251 (MQ=255)                                                                                                                                                                                                             aaCCCTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCTGCGCTTTATGTCACGCCGGGACATATCGAAAAGATGGCCTATATTGTTGGCGCGTAACGccacagctc                              >  2:72304/1‑242 (MQ=255)                                                                                                                                                                                                                 cTGGCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATTGAAATTATGGGTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCgg                         >  1:442247/1‑251 (MQ=255)                                                                                                                                                                                                                   ggCTGATGCCTTATACCGACGAAACGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTa                       <  1:115255/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaaCGCTGAGAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCGGCGCTTAATGCCACGCCGGAACATATAGGAATGATGGCTCATCTTGTTTCCCCGTCTTGCGAAAGCAGAGCGGGAAAAAACCGAGCGCCGTcgcgg  >  2:454221/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGAC‑TCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGA  >  NC_000913/498570‑499050

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