Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	389,601	C→T	D43N (GAC→AAC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	389,601	0	C	T	100.0%	235.0 / NA	64	D43N (GAC→AAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/42);  total (22/42)

GTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAT  >  NC_000913/389356‑389843                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gtATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCaa                                                                                                                                                                                                                                               <  1:109287/251‑1 (MQ=255)     tCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAAttt                                                                                                                                                                                                                                            <  1:98244/250‑1 (MQ=255)        caGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttc                                                                                                                                                                                                                                        <  1:236692/251‑1 (MQ=255)             gCAGGTGTCGGACAGAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCGTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACTGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGAGGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGTCTTTGTAGTTGTCAATTTCGTCTTcaa                                                                                                                                                                                                                                   <  2:22009/251‑1 (MQ=255)              cAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACTCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaac                                                                                                                                                                                                                                  >  2:477298/1‑251 (MQ=255)                       gACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATc                                                                                                                                                                                                                           <  1:435278/249‑1 (MQ=255)                               gATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc                                                                                                                                                                                                                  >  1:29286/1‑250 (MQ=255)                                  catttctggcacctgttgtttGCGGATGCGGGTCATTCGCGGCACCATGCCATTTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGTCGATGCGTTTAATTTCGCGTTCCAGATGTTTCTCTGGAATGCGCATCTCGCCTGGCATGGCTGCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGGTCGGCAACACCAGGt                                                                                                                                                                                                              <  2:449290/232‑1 (MQ=255)                                         ggCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGTAACACCAGGTCGtt                                                                                                                                                                                                          <  1:251648/248‑1 (MQ=255)                                         ggCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTa                                                                                                                                                                                                         >  1:385084/1‑249 (MQ=255)                                         ggCACGGTCGGCTTGCAGATGCGCGACATCCGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAg                                                                                                                                                                                                       <  2:127926/251‑1 (MQ=255)                                                 cTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTTTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtg                                                                                                                                                                                               >  2:417186/1‑251 (MQ=255)                                                   gCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGACCGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCACTTTCTTCTTCAACAAAGATCGGCAACACCAGGGCGTTAAGGCTAAGagtt                                                                                                                                                                                             >  2:302241/1‑251 (MQ=255)                                                      tGCAGATGCGCGACATGCGCGCCCCCAGTCNATCTTCCCTCCAGTCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGGCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGtgttgt                                                                                                                                                                                           <  2:248369/250‑1 (MQ=255)                                                       gCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct                                                                                                                                                                                         <  1:116437/251‑1 (MQ=255)                                                       gCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtc                                                                                                                                                                                          <  2:29286/250‑1 (MQ=255)                                                          gATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGTCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCtt                                                                                                                                                                                      >  1:256039/1‑251 (MQ=255)                                                                 gACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACCAAGATCGGCAACACCAGGTCGTTAAGGCTAGGTGTTGTCTCTTCAAACAt                                                                                                                                                                               >  2:107129/1‑251 (MQ=255)                                                                 gACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCAGCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAt                                                                                                                                                                               <  1:401951/251‑1 (MQ=255)                                                                        gcgcCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgcgc                                                                                                                                                                        <  1:177970/251‑1 (MQ=255)                                                                                aGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAgg                                                                                                                                                                <  1:431897/251‑1 (MQ=255)                                                                                 gTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAgg                                                                                                                                                                >  1:451784/1‑250 (MQ=255)                                                                                          cccGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcgc                                                                                                                                                      <  2:65844/251‑1 (MQ=255)                                                                                            cGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcgc                                                                                                                                                      >  1:294616/1‑249 (MQ=255)                                                                                              ccAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGc                                                                                                                                                  <  2:368750/251‑1 (MQ=255)                                                                                                aGGCATCGCTGCAGGTTTCATCGGTATGGTGCGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGCATGCGAATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCACTTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGGCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcg                                                                                                                                                 >  2:38863/1‑250 (MQ=255)                                                                                                  gCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcgacg                                                                                                                                              >  1:138590/1‑251 (MQ=255)                                                                                                          gccgtgttcatcgtgctGGTGAGAGATTACAAAAGTCATCACTGAACGAATACCGGCGTTGGCGATGCGTTCATTGTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTGCTTCTGCAACAAAGATCGGCAACGCCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGAAGGCGACGAGGGCg                                                                                                                                        <  2:342962/233‑1 (MQ=255)                                                                                                              gTTTCTTCGGTATGGTGAGAGATGGCAAAAGTCATCACGGAACGAATACCGGCGTGGGGTATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCTCGTCTGGCATGGCTTCAACGGCTTTGTAGTTGGCAATTTCTTCTGCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGAt                                                                                                                                  <  2:463307/251‑1 (MQ=255)                                                                                                                             tGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGt                                                                                                                   >  2:103270/1‑251 (MQ=255)                                                                                                                                 agaTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGc                                                                                                               <  2:451784/251‑1 (MQ=255)                                                                                                                                   aTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCt                                                                                                             <  2:395626/251‑1 (MQ=255)                                                                                                                                     gCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGa                                                                                                           <  1:360230/251‑1 (MQ=255)                                                                                                                                        aaaaGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTTGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGt                                                                                                        <  1:20644/251‑1 (MQ=255)                                                                                                                                           aGTCAGCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTtg                                                                                                     <  1:480350/251‑1 (MQ=255)                                                                                                                                            gTCAACACGGAACGAATACCGGCGTTGGCGATGCGGTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACGCCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGTCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTtgt                                                                                                    <  2:138590/251‑1 (MQ=255)                                                                                                                                              catcaCGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTtgtg                                                                                                   <  1:3879/250‑1 (MQ=255)                                                                                                                                                   cGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCCGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCg                                                                                             >  1:136701/1‑251 (MQ=255)                                                                                                                                                             aCCGGCGTTGGCGATGCGTTCAATTTCGCGNGCCAGATGTTTCTCTGGAATGCGCATCTCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGGCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATa                                                                                    <  2:231690/250‑1 (MQ=255)                                                                                                                                                                   gTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATAcc                                                                             <  1:205170/251‑1 (MQ=255)                                                                                                                                                                         gATGCGTTCAATTTCGCGTGCCAGATGCTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGc                                                                       >  1:400732/1‑251 (MQ=255)                                                                                                                                                                         gATGCGTTCAATTTCGCCTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCGTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGc                                                                       <  1:284340/251‑1 (MQ=255)                                                                                                                                                                           tGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCCTGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGACTGCCTGATGTTTGTGGAATCGAGGGGCATCGTATACCTGCAGc                                                                       >  1:155438/1‑249 (MQ=255)                                                                                                                                                                              gtttaaTTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGTCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTGCAACAAAGAGCGGCAACACCAGGTCGTTAAGGCTTAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCGGTCATGGTCTGCCGGATGTTTGTGGAATCAAGGGGCATAGTATACCTGAAGCAGGGt                                                                  <  2:197451/247‑1 (MQ=255)                                                                                                                                                                                            gCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                                                    >  2:18103/1‑251 (MQ=255)                                                                                                                                                                                               aGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGTCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                                                 <  1:465764/251‑1 (MQ=255)                                                                                                                                                                                                gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGGTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTTATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAg                                                <  1:466479/251‑1 (MQ=255)                                                                                                                                                                                                 aTGTTTCTCTGGAATGCGCATCACGCATGGCATGGATTCAACGGCTTTGTAGTTGTCAATTTCTGCTTCAACAAAGATCGGCCCCACCAGGTCGTTAAGGCTATGTGGTGTCTCTTCAAACATAGCGCGCAGCGCCGGAGATTTGCGCAGGGGACGAGGGCGTTGGATTAAGGCTGTCATGGTCTGCCTGATGTTTGTGGAATCGGGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAgg                                               >  2:190587/1‑250 (MQ=255)                                                                                                                                                                                                  tGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCATACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGAATAGTATCACTGAAGCAGGGTAGGGATGTTTTACGAAAGtt                                              >  2:205051/1‑251 (MQ=255)                                                                                                                                                                                                   gTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGtt                                              <  1:485638/250‑1 (MQ=255)                                                                                                                                                                                                     ttCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTc                                           >  2:137516/1‑251 (MQ=255)                                                                                                                                                                                                      tctcTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATTCCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTc                                           <  1:281498/250‑1 (MQ=255)                                                                                                                                                                                                       ctctGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAAGGAATACTTGAAGCAGGGAAGGGAGGTTTTACGAAAGTTGTcgc                                         >  1:188958/1‑251 (MQ=255)                                                                                                                                                                                                         ctGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcg                                          <  1:302241/248‑1 (MQ=255)                                                                                                                                                                                                           ggAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgcg                                        <  1:107129/248‑1 (MQ=255)                                                                                                                                                                                                           ggAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATg                                     <  1:18103/251‑1 (MQ=255)                                                                                                                                                                                                           ggAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATg                                     >  1:394307/1‑251 (MQ=255)                                                                                                                                                                                                           ggAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAAGTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGATAAGTGTTGTCTCTGCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGAt                                      <  2:103281/250‑1 (MQ=255)                                                                                                                                                                                                            gAATGCGCATCCCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGGAGGGATGTTTTACGAAAGTTGTCGCGATGt                                    >  2:264859/1‑251 (MQ=255)                                                                                                                                                                                                                             tgCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAAt                     <  1:190587/248‑1 (MQ=255)                                                                                                                                                                                                                               cATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGGGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGaa                 <  2:204546/251‑1 (MQ=255)                                                                                                                                                                                                                                      tCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAgg           <  1:137516/250‑1 (MQ=255)                                                                                                                                                                                                                                            gCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGAc    <  2:385084/251‑1 (MQ=255)                                                                                                                                                                                                                                               ttGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTTGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAt  <  1:417186/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAT  >  NC_000913/389356‑389843

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