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	97.1%	253.5 / NA	70	D43N (GAC→AAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/37);  minor base G (0/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTT  >  NC_000913/389356‑389846                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      gtATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCa                                                                                                                                                                                                                                                   >  1:19614/1‑250 (MQ=255)        caGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTGCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTGCAACGGCTTTGTAGTTGTCAATttcttc                                                                                                                                                                                                                                           <  1:509960/251‑1 (MQ=255)           aggCAGGTGTCTGACATAACGATCATTTCTGTCTCGGTCTGCTTGCAGATGCGCGTCATACGCGCCCCCAGTCCATCTTCCCGCCAGGCATAGCTGCCGGTTTCAGCTGTATGGTGAGAGATGCCCAAAGTCATCGCGGAACGAATACCGGGGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCGCGCCTGGCATGGCTTCAACGGCTTTGTAGTGGTGAATttcttcttc                                                                                                                                                                                                                                        <  2:287993/249‑1 (MQ=255)            aGCAGGCGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttcttc                                                                                                                                                                                                                                        >  2:14831/1‑250 (MQ=255)              cAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaac                                                                                                                                                                                                                                     <  1:256220/251‑1 (MQ=255)                 gtgtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                                   >  2:245507/1‑250 (MQ=255)                  tgtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                                   >  2:461260/1‑249 (MQ=255)                   gtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa                                                                                                                                                                                                                                   >  1:297381/1‑248 (MQ=255)                   gtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGa                                                                                                                                                                                                                                <  1:266029/251‑1 (MQ=255)                     cTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTtaacaa                                                                                                                                                                                                                                   >  2:272447/1‑246 (MQ=255)                       gACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATTCGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATc                                                                                                                                                                                                                              <  2:430242/249‑1 (MQ=255)                         cATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCa                                                                                                                                                                                                                          <  2:426745/251‑1 (MQ=255)                              cGATCATTTCTGGCACTGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCGTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATTCCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAAGGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGAGCGGCAACAcc                                                                                                                                                                                                                     <  2:498705/251‑1 (MQ=255)                               gATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc                                                                                                                                                                                                                     >  1:584695/1‑250 (MQ=255)                                ataaTTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTATTTGTCAATTTCTTCTTCAACAAAGATCGGCacccac                                                                                                                                                                                                                     >  2:346864/4‑244 (MQ=255)                                aTCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc                                                                                                                                                                                                                     >  1:512189/1‑249 (MQ=255)                                  cATTCCTGGTACGGTCTGCTTGCAGAGGCGCGACATACGCGCAACCAGTCCATGTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGACCTAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTGTCTCTGGAATGCGCTTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCACCACCAgg                                                                                                                                                                                                                  <  2:27857/250‑1 (MQ=255)                                   aTTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTc                                                                                                                                                                                                                <  1:44433/251‑1 (MQ=255)                                            aCGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCt                                                                                                                                                                                                       <  1:599547/251‑1 (MQ=255)                                             cgttCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGAGCGGCAACACCAGGTCGTTAAGGCTa                                                                                                                                                                                                      <  2:339442/248‑1 (MQ=255)                                               gTCTGCTTTCAGATGCGCGACCTACGCGCCACCAGTCCATCATCCCGCCATGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATCCCGGCGTTGTCGATG‑CGTTCAATTTCGCGTGCCAGATGTGTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAg                                                                                                                                                                                                    <  2:629794/251‑1 (MQ=255)                                                     ttGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCCGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCATGGCATGGCGTCAACGGCGTTGTAGTTGTCAATTTCTTCTTCGACAAAGATCGTCAACACCAGGTCGTTAAGGCGAAGtgttgt                                                                                                                                                                                              >  2:83790/1‑251 (MQ=255)                                                       gCAGAGGCGCGACATACGCGCCACCAGTCCATCATCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGGCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtc                                                                                                                                                                                             <  2:301887/250‑1 (MQ=255)                                                         ggATGCGCGACATAAGCGCCACCATTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAATATCGGCAACACCAGGTCGTTAAGGTTAAGTGTTGTctct                                                                                                                                                                                          <  2:448076/250‑1 (MQ=255)                                                          gATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCtt                                                                                                                                                                                         >  1:435985/1‑251 (MQ=255)                                                             gcgcGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGACTACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTcc                                                                                                                                                                                       >  2:46603/1‑249 (MQ=255)                                                                  aCATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAt                                                                                                                                                                                  >  1:248782/1‑250 (MQ=255)                                                                   cATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAt                                                                                                                                                                                  >  2:289236/1‑249 (MQ=255)                                                                    aTACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgc                                                                                                                                                                               <  1:245507/251‑1 (MQ=255)                                                                      aCGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcg                                                                                                                                                                              <  1:176458/250‑1 (MQ=255)                                                                       cgcgCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgcg                                                                                                                                                                            <  2:469029/251‑1 (MQ=255)                                                                             accaGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCgcgcagcgc                                                                                                                                                                      <  1:280552/251‑1 (MQ=255)                                                                                     acctacccGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACTGAACGACTACCGGCTGTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTGAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGGGAtt                                                                                                                                                              <  2:75034/246‑1 (MQ=255)                                                                                         tCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATTGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAATATCGGCAACCCCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgc                                                                                                                                                           <  2:563585/250‑1 (MQ=255)                                                                                         tCCCGCCAGGCATCGCTGCCGGTTTCATCGGTAGGTTGAGAGATGCCAAAAGTCAGCACGGAACGAATACCGGCGTTGGCGATT‑CGTTCAATTTCGCGTGCCATATGGTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTGCAACAAAGATCGGCAACGCCAGGTCGGTAAGGCTAAGGGTTGTCTGTGCAAACATAGCGGGCAGCTCAGGAGAttt                                                                                                                                                             <  2:184775/248‑1 (MQ=255)                                                                                          cccGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcgc                                                                                                                                                         >  1:55846/1‑251 (MQ=255)                                                                                                 ggCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcgac                                                                                                                                                  <  2:254381/251‑1 (MQ=255)                                                                                                  gCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcgacg                                                                                                                                                 <  2:95999/251‑1 (MQ=255)                                                                                                    aTCGCTGCAGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGcgacg                                                                                                                                                 >  1:412612/1‑249 (MQ=255)                                                                                                           ccGGTTTCCTCGGAATGGAGAGCGATGCCAAAAGTCATCACGGAACGAACACCGACGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTATCACTTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCAAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGAGCAGGCGACGAGGGCGTTg                                                                                                                                        >  2:162478/1‑251 (MQ=255)                                                                                                             ggTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTgg                                                                                                                                       >  1:288575/1‑250 (MQ=255)                                                                                                             ggTTTCATCGGTATGGCGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATTTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTCTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTACGTGTTGACTCTACAAACATAGCGAGCAGCGCAGGAAATTTGCGCAGGCGACGAGGGCATTgg                                                                                                                                       >  2:604025/1‑250 (MQ=255)                                                                                                                     cGGTATGGTGAGAGATGCCAAAAGTCATCACGGCACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACAAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAACCATAGCGCGCAGCGCAGGAGGTTTGCGCAGGCTACGCGGGCGTTGGATTAAgcct                                                                                                                              >  2:425414/1‑248 (MQ=255)                                                                                                                         tgggTGAGAGATGACAAAAGACGTCAAGGATCGAATACCGGCGTTGGCGGTG‑CGTTCAATGTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGCCATGGCTTCAACGGCTTTGTAGTTGTAAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTACAAACATAGCGCGCAGCGCAGGAGATTTGCACAGGCGACGAGGGCGTTGGATTAAGTCTGTc                                                                                                                           <  2:596337/248‑1 (MQ=255)                                                                                                                         aTGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTc                                                                                                                           <  1:289236/250‑1 (MQ=255)                                                                                                                         aTGGTGAGAGATGCCAAAAGCCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTc                                                                                                                           >  1:568732/1‑250 (MQ=255)                                                                                                                                        aaaaGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGt                                                                                                           <  2:41296/251‑1 (MQ=255)                                                                                                                                             tcatcaCGGAACGAATACCGGCGTTGGCGATA‑CGTACAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGCGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTtgtg                                                                                                      >  1:467583/1‑251 (MQ=255)                                                                                                                                                  ccGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATc                                                                                                 <  1:418619/250‑1 (MQ=255)                                                                                                                                                     gAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGACTCg                                                                                                >  2:90736/1‑249 (MQ=255)                                                                                                                                                              ccGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGt                                                                                     <  1:78748/251‑1 (MQ=255)                                                                                                                                                                gtcGTTGGCGATGATGTT‑AATGTCGCGTGCCAGATGTTGCTCTGGAATGCGCATCACGACGGGCATGGCTTCAACGGCTTTGTAGTTGGCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGGCGTTAAGGCTGAGTGGTGTCTCGTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTGGGATTAAGTCTTTCATTGTCTGCCTGATGTCTGGGTAATCGGGGGGCATAGtat                                                                                   <  2:156444/249‑1 (MQ=255)                                                                                                                                                                  ggtttgCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATAc                                                                                 <  2:297381/246‑1 (MQ=255)                                                                                                                                                                        cGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAg                                                                           <  1:134975/251‑1 (MQ=255)                                                                                                                                                                         ggaggctttCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTGGTCAATTTCTGCTGCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACAGAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAg                                                                           <  2:568450/244‑1 (MQ=255)                                                                                                                                                                         gATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGc                                                                          >  2:170085/1‑251 (MQ=255)                                                                                                                                                                          aTG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGc                                                                          >  2:150508/1‑250 (MQ=255)                                                                                                                                                                            g‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAg                                                                        <  2:248782/250‑1 (MQ=255)                                                                                                                                                                            g‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAg                                                                        <  1:294749/250‑1 (MQ=255)                                                                                                                                                                                  cAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAgg                                                                  >  2:239527/1‑251 (MQ=255)                                                                                                                                                                                    aTTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTAAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACATGAAGCAGGGTAggg                                                                 >  2:312092/1‑250 (MQ=255)                                                                                                                                                                                      ttCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGAt                                                               <  1:22304/250‑1 (MQ=255)                                                                                                                                                                                           gTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTCTACCTGAAGCAGGGTAGGGATGtttt                                                          >  2:390602/1‑250 (MQ=255)                                                                                                                                                                                                          ctGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTTTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgc                                            <  1:293292/249‑1 (MQ=255)                                                                                                                                                                                                                                   ggCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAaga                  <  2:106143/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAag                   >  1:566120/1‑248 (MQ=255)                                                                                                                                                                                                                                    gCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAag                   >  1:429885/1‑248 (MQ=255)                                                                                                                                                                                                                                                ttGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAt     <  1:83790/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tGTAGTTGTCAATTTATTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCt   >  2:437333/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCtt  <  1:165463/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATG‑CGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTT  >  NC_000913/389356‑389846

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