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%	166.6 / NA	55	P21L (CCT→CTT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/24);  minor base A (0/1);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATG  >  NC_000913/497873‑498354                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                gAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCGACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGc                                                                                                                                                                                                                                          <  2:352413/250‑1 (MQ=255)    gggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAATCCTCAATCTGCGCTATTCCTTTTTCGGATTTGACCGCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATCATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGTGATCCTGCTGGCAAACCTGGGTACGCCTGATGCCCCCACACTTG‑AAGCGGTa                                                                                                                                                                                                                                      <  1:338714/251‑1 (MQ=255)     ggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTaa                                                                                                                                                                                                                                     <  1:389376/251‑1 (MQ=255)     ggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTGCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTGG‑AAGTGGTaa                                                                                                                                                                                                                                     <  2:256809/251‑1 (MQ=255)            gCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTa                                                                                                                                                                                                                              <  1:17286/251‑1 (MQ=255)              ccGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCGCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATc                                                                                                                                                                                                                            <  2:393101/251‑1 (MQ=255)                gCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTg                                                                                                                                                                                                                          <  1:183070/251‑1 (MQ=255)                     aGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAAc                                                                                                                                                                                                                      >  1:332391/1‑250 (MQ=255)                          aCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAAtttt                                                                                                                                                                                                                <  1:426968/251‑1 (MQ=255)                                               ggATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCCCTTTTCCGCTACAATTATCACCAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATc                                                                                                                                                                                                                            >  1:525293/1‑218 (MQ=255)                                               ggATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTAAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTg                                                                                                                                                                                           >  1:110568/1‑251 (MQ=255)                                               ggATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAAGTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCGGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATc                                                                                                                                                                                                                            <  2:525293/218‑1 (MQ=255)                                               ggATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTGTTCTGATTTGACCTCGCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCTGCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACATG‑AAGCGGTAAAAAGCTAGCTGAAACAATGTGTAAGCGACAGACGCGTGGTTg                                                                                                                                                                                           <  2:107908/251‑1 (MQ=255)                                                            gCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGGTGATACCTTACGgtt                                                                                                                                                                              >  1:494000/1‑251 (MQ=255)                                                                   aCAAAAACATCAATCTGCCCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCAATTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTCAAACCGGAATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGATATCTGAAACCATTTTTAAGCGACAGACGCGTGGTTGCTACCTCAAGGTTGTTTtgg                                                                                                                                                                       >  2:254012/1‑251 (MQ=255)                                                                       aaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGc                                                                                                                                                                   >  1:510557/1‑251 (MQ=255)                                                                        aaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGcc                                                                                                                                                                  <  2:9752/251‑1 (MQ=255)                                                                          ccTCAATCTGCGCTATTCCTTTTTCTGATTGGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCGATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTGGTTATGGTGGCCAt                                                                                                                                                                <  2:211889/251‑1 (MQ=255)                                                                              aaTCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGCTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTGTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtgct                                                                                                                                                            <  2:366535/251‑1 (MQ=255)                                                                                tCTGCGCTATTCTTTTTTCTGATTTGACCTCTCACAGCCATTCGTTCTTCTTCCTCATTTTTCCGCTACAATTCTCAACACGTTGAATCGCTAAGAGGCGGTAATGCGTCAGACTAAACCCGGTCTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCACCCACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTCTGGTGGCCATtgctgc                                                                                                                                                          >  2:265608/1‑251 (MQ=255)                                                                                        aTTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTg                                                                                                                                                  >  2:43294/1‑251 (MQ=255)                                                                                                           ccTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgc                                                                                                                               <  1:212672/251‑1 (MQ=255)                                                                                                            ctctCACAGCACTTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgcc                                                                                                                              >  2:42559/1‑251 (MQ=255)                                                                                                            ctctCACAGCAATTCGTTCTTCTTACTCACTTTTCCGCTACACTTATCAACAAGATGAATCGCTCAGAGGCGGTAATGCGTCAGACTCAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGAGCTcgcc                                                                                                                              >  2:376217/1‑251 (MQ=255)                                                                                                            ctctCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgcc                                                                                                                              >  2:426177/1‑251 (MQ=255)                                                                                                                     cAATTAGTTCTTCTTCCTCACGTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGGAAGGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTGG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGTTTGTTATGGTGTCCATTGCGACGCTGCGTGATTTTGCCGCGGCGCTCGCCGCGTGTGGc                                                                                                                     <  2:506325/251‑1 (MQ=255)                                                                                                                              cttcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGTCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGCCAGACGCGTGGTTGATACCTCCCTGTTGTTATGGTGGCCCTTGCTGCGCGGCGGGATTTTGCCCCTGCGCTCGCCGCGCGTGGCGACGCttgg                                                                                                            >  2:2826/1‑248 (MQ=255)                                                                                                                                tcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATg                                                                                                          >  1:243074/1‑251 (MQ=255)                                                                                                                                                 cTACACTTATCAACAAGTTGACTCGATAAGAGGCGGTAATGCGTCCGACTAACACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTCTCTGAAACAATTTTTCAGCGACCGACGCGTGGTTGATACCTCACGGTTGGTTTGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGGGCTCGCCGCGAGGGGCGAAGCTGGATGCCTCTGTCtggatggc                                                                                          >  2:273848/1‑249 (MQ=255)                                                                                                                                                   aCAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACCATTTTTAAGCGACAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCCCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTTGATGGAAgg                                                                                       >  2:176493/1‑251 (MQ=255)                                                                                                                                                     aaTTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                     <  2:332391/251‑1 (MQ=255)                                                                                                                                                         aTCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACAC‑TGTAAGCGGTAAAGCGCTATCTGAAACAATTTTTAAGCTACAGGCGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCGTTGCTGCTCGGCGTGATGTTGCCGCTGCGGGCGTCGCGTTTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTc                                                                                 <  1:526019/251‑1 (MQ=255)                                                                                                                                                               aaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCt                                                                           <  2:145947/251‑1 (MQ=255)                                                                                                                                                               aaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGGATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCAAACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGGGGCCCTTGCTGCGCGGCGTGATTTTTCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATAGAAGGGGGCTCGCCGCt                                                                           >  1:907/1‑251 (MQ=255)                                                                                                                                                                aGTTGAATCGATAAGAGGCGGTGATGAGGCAGACTAAACCCGGTATCCTGATGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑ACGCGGAAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGGGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgccgcgg                                                                          >  2:408219/1‑249 (MQ=255)                                                                                                                                                                    gAATCGATAAGCGGCGGTAATGCGTCAGACTAAAAACGGTATCATGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGCAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCGGTCTGGGTGGAAGGGGGCTCGCCGCTGaagg                                                                      >  2:93574/1‑248 (MQ=255)                                                                                                                                                                     aaTCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGt                                                                     <  1:426177/251‑1 (MQ=255)                                                                                                                                                                                 gCGGTAAGGCGGCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCa                                                         <  2:184045/251‑1 (MQ=255)                                                                                                                                                                                         gCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagc                                                 >  1:545434/1‑251 (MQ=255)                                                                                                                                                                                              aGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgcg                                            >  1:434693/1‑251 (MQ=255)                                                                                                                                                                                                aCTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGGCGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCt                                          >  2:132612/1‑251 (MQ=255)                                                                                                                                                                                                   aaaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTgg                                        <  1:43294/250‑1 (MQ=255)                                                                                                                                                                                                      aCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGcac                                     >  2:471536/1‑250 (MQ=255)                                                                                                                                                                                                       ccGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGGCGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACaa                                   <  1:132612/251‑1 (MQ=255)                                                                                                                                                                                                       ccGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTGG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACaa                                   <  1:59849/251‑1 (MQ=255)                                                                                                                                                                                                         ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACg                                 >  1:152552/1‑251 (MQ=255)                                                                                                                                                                                                          gTATCCTGCTGGCAAACCTGGGTACACCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGAGGTTGATACCTCACGGTTGTTATGCTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGACGGGGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCAAAACGt                                >  1:512250/1‑251 (MQ=255)                                                                                                                                                                                                                       aaaCCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAACACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTCATAACTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCATCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGc                   >  1:258389/1‑251 (MQ=255)                                                                                                                                                                                                                       aaaCCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGc                   >  1:207601/1‑251 (MQ=255)                                                                                                                                                                                                                       aaaCCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTGTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTGGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGTAAGTTGGCTCGCCGCTGATGGGTTAGATTCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGc                   <  1:376217/251‑1 (MQ=255)                                                                                                                                                                                                                        aaCCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGcc                  <  1:150127/251‑1 (MQ=255)                                                                                                                                                                                                                        aaCCTGGGTACGCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCAATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGcc                  >  2:28353/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCCCGATGCCCCCACACTTG‑AAGCGGTAAACCGCTATCGGAAACAATTTTTAAGCGACAGACGCGTGGTTGATCCCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTgg       <  2:60838/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCCCGATGCCCCCACACTTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTCCCGGAGAGGCCCGGAGCGCTgg       >  2:59841/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTGCCCCCACACTTG‑AATCGGTAAACCGCTATCTGAAACAATTTTTACGCGACAGACTCGTGGTTGATACCTCACGGTTGTTATGGTAGCCATTGCTGCGCGGCGCGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTCTATGCCTCTGTCTGGATGGAAGGCGGCTCGCCGCTGATGGTTTACAACCGCCAGCAACCGCAGGCGCTGGCACAACGTTTACCGGACATGCTCTTAGCGCTGCGAATg  >  1:515197/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                GAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTG‑AAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATG  >  NC_000913/497873‑498354

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