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.1%	195.9 / NA	54	P21L (CCT→CTT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/18);  minor base G (0/1);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGA  >  NC_000913/497874‑498351                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              aGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCgg                                                                                                                                                                                                                                       >  1:126574/1‑251 (MQ=255)                   aaGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAAc                                                                                                                                                                                                                    <  1:236183/251‑1 (MQ=255)                       aaaCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAAtt                                                                                                                                                                                                                >  1:451535/1‑251 (MQ=255)                           ttATCCGCAATCTCAGGCCGGATATTCCTTCGGCCGTTTAAAAACACCTCAATCTGAGCTATTCCTTTTTCTGATTTGACATCTCACAGCAATTAGTTCTTCTTCATCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGGCATCAGGCAGACTAAAACCGGTCTCCTGCTGGCAACACTGGGTACGCCCGAGCCCCCCACACTTGAAAAGGTACAACGCTATCTG‑AACCAAttttt                                                                                                                                                                                                             >  2:479736/1‑250 (MQ=255)                             aTCAGCAATCTCAGGCAGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGGACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTa                                                                                                                                                                                                            >  1:409537/1‑249 (MQ=255)                               cAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGc                                                                                                                                                                                                         <  1:178568/250‑1 (MQ=255)                                      ctcAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGAcg                                                                                                                                                                                                 >  1:331205/1‑251 (MQ=255)                                        cAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGAc                                                                                                                                                                                                  <  2:60068/248‑1 (MQ=255)                                            ccGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGt                                                                                                                                                                                           <  1:359235/251‑1 (MQ=255)                                            ccGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACGTGAAGCGGTAAAACGCGATCGG‑AAACAATTTTGAAGCGACAGACGCGTGGt                                                                                                                                                                                           <  2:47175/251‑1 (MQ=255)                                                               tttACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgttt                                                                                                                                                                        >  2:165786/1‑250 (MQ=255)                                                                  aCAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACACGTGGTTGATACCTCACGGTTGTTAtgg                                                                                                                                                                     >  1:207745/1‑251 (MQ=255)                                                                     aaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACCGACGCGTGGGTTATACCTCCCCGTTGTTAtggtgg                                                                                                                                                                  >  2:378125/1‑251 (MQ=255)                                                                       aaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGGATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGcc                                                                                                                                                                >  1:74773/1‑251 (MQ=255)                                                                        aCCTCAATATGCGCTATTCCTTTTTCTGATTTGACCTCTCCCAGCAATTAGTTCTTCTTCCTCACTTTTCAGCTCCAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGAAGCGTGGTTGATACCTCACGGTTGTTATGGTGGccc                                                                                                                                                               >  1:63126/1‑250 (MQ=255)                                                                                  gcgcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtgctg                                                                                                                                                         <  2:58486/247‑1 (MQ=255)                                                                                       aTTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTg                                                                                                                                                <  1:272145/251‑1 (MQ=255)                                                                                           cTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTCTCCTGCTGGCAAACCCGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGGTGATACCTCACGGTTGTTATGGTGGCCATTTCTGCGCGGCGTGAttt                                                                                                                                            >  2:272631/1‑251 (MQ=255)                                                                                                tCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGcc                                                                                                                                        <  1:239519/250‑1 (MQ=255)                                                                                                   gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCACCAAGTTGAATCGATAAGAGGCGGTAATGCGTCCGACTAAACCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCCATCTG‑AAACAATTTTTAAGCGACATACGCGTGGTTGATACCTCACGGGTGTTATGTTGGCCATTTCTGCCCCGGGGGATTTTGCCGCTg                                                                                                                                    >  2:438187/1‑251 (MQ=255)                                                                                                        gACCTCTCACAGCAATTAGTTCTTCTTCCTTACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAAACGGTATCCTGCTGGCAAACCTGGGTACGGCCACTGCCCACACACATGAAGCGGTTAACCGCTATCTGCAAACA‑TTTTTAAGCGACAGAAGCGTGGTTGATACCTCCCGGGTGTTTTTGTTGCCATTGCTGGGCGGCGTGTTTTTGcccctttcctc                                                                                                                               >  2:419536/1‑242 (MQ=255)                                                                                                              tcACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCcgcg                                                                                                                         >  1:290184/1‑251 (MQ=255)                                                                                                                    cAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGc                                                                                                                   <  1:261813/251‑1 (MQ=255)                                                                                                                      aTTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGa                                                                                                                 <  2:325220/251‑1 (MQ=255)                                                                                                                             cttcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTAAGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTCTCTG‑AAACAATTTGTCAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGCTTTTGCCGATGCGCCCGCCGCGTGTGGGGAAGcggct                                                                                                          >  2:167781/1‑247 (MQ=255)                                                                                                                                          tttCCGCTACAATTATCAACCAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCAGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGGGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCACTTTCTGCGCGGCGTGATTTTGCCGCTGCGCCTGCCGCGGGGGGCGGAGCTGGATGCCTctggctgg                                                                                             >  2:417300/1‑251 (MQ=255)                                                                                                                                               gCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCtggatgg                                                                                         >  1:283911/1‑250 (MQ=255)                                                                                                                                                    aaTTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                   >  1:116335/1‑251 (MQ=255)                                                                                                                                                        aTCAACAAGTTGAAACGATAAGAGGCGGTCATGCGTCAGCCTAAAACCGGTATCCTGCTGGCAAACCTGGGCACGCCCGATGCCCCCCCACTTGAAGCGGTAAAAAGATATCTG‑AAACAATTTTTCAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTTGCCATTGCTGCGCGGCGTGATTTTGCCGCCGCGCTCGCCGCGGGTGGCGAAGCTGTATGACCCTGTCTGGGTGGAAGGGGGCTc                                                                               >  2:80031/1‑251 (MQ=255)                                                                                                                                                             caaGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgccgc                                                                          <  2:355828/251‑1 (MQ=255)                                                                                                                                                                           aaGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGACGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGCCGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGGAGGGGGCTCGCCGCTGATGGTT‑TACAgc                                                           >  2:108666/1‑251 (MQ=255)                                                                                                                                                                             gagGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAgccg                                                         >  1:174230/1‑251 (MQ=255)                                                                                                                                                                              agGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCt                                                                                                                                     >  2:70301/1‑175 (MQ=255)                                                                                                                                                                              agGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTGGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCt                                                                                                                                     <  1:70301/175‑1 (MQ=255)                                                                                                                                                                               ggCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCTGATGGTT‑TACAgccgcc                                                       >  2:376296/1‑251 (MQ=255)                                                                                                                                                                                 cGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGcc                                                                                                      <  2:156436/203‑1 (MQ=255)                                                                                                                                                                                 cGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGcc                                                                                                      >  1:156436/1‑203 (MQ=255)                                                                                                                                                                                    tAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCa                                                   >  1:104392/1‑250 (MQ=255)                                                                                                                                                                                         cGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAAcagc                                              >  1:192954/1‑250 (MQ=255)                                                                                                                                                                                          gTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAAcagc                                              >  2:364379/1‑249 (MQ=255)                                                                                                                                                                                                cTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGGGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTg                                      >  2:23907/1‑251 (MQ=255)                                                                                                                                                                                                  aaaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGc                                    >  1:338257/1‑251 (MQ=255)                                                                                                                                                                                                  aaaaCCGGTACCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCTGAGGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGc                                    <  2:74773/251‑1 (MQ=255)                                                                                                                                                                                                   aaaCCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGGGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGGGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGcc                                   >  2:149433/1‑250 (MQ=255)                                                                                                                                                                                                        ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCCGGCGCTGGcacaccg                              >  2:146273/1‑248 (MQ=255)                                                                                                                                                                                                           aTCCTGCGGGCAAACCTGGGTACGCCCTATGCCCCCACACTGGAAGCGGTAAAACGCTAGCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGTTTGTTATGGTGGCCATTGCGGCGCGGCGTGATTTTGCCGCTTCGCTCGCCGCGTGTGGCGAATCTGTATGCCTCTCTCGGGGTGGAGGGTGTCTCGCCGCTGA‑GGTTGTTCAGCCGCCAGCAACAGCAGACGCTGGCTCAACGttt                           <  2:63126/251‑1 (MQ=255)                                                                                                                                                                                                              ctgctgGCAAACCTAGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTCAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTAcc                        >  2:187287/1‑251 (MQ=255)                                                                                                                                                                                                                 ctgGCAAACCTGGGGACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCCCAACGTTTACCggg                     >  2:56038/1‑250 (MQ=255)                                                                                                                                                                                                                      aaaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGc                <  1:378125/251‑1 (MQ=255)                                                                                                                                                                                                                       aaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACGTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGAt                                                                      <  1:435407/197‑1 (MQ=255)                                                                                                                                                                                                                       aaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACGTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGAt                                                                      >  2:435407/1‑197 (MQ=255)                                                                                                                                                                                                                       aaCCTGGGTACGCCCGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGcc               <  1:221074/251‑1 (MQ=255)                                                                                                                                                                                                                                   cccGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTggg   <  1:165786/251‑1 (MQ=255)                                                                                                                                                                                                                                     cGATGCCCCCACACTTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTCACGGAGATGCCCGTAGCGCTGGGa  >  2:106040/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              AGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCTGAAGCGGTAAAACGCTATCTG‑AAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTT‑TACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGA  >  NC_000913/497874‑498351

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