Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	157,585	C→T	D50N (GAT→AAT)	folK ←	2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	157,585	0	C	T	98.7%	240.7 / NA	74	D50N (GAT→AAT)	folK	2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/31);  minor base G (0/1);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGT  >  NC_000913/157338‑157829                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     aTTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTGCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCGTTGCGGACGCGTCCTTGCTGCAATTGAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTg                                                                                                                                                                                                                                                   <  1:94232/251‑1 (MQ=255)            cGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCg                                                                                                                                                                                                                                        <  2:548595/251‑1 (MQ=255)               ccACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTg                                                                                                                                                                                                                                     <  2:202009/251‑1 (MQ=255)                     cATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGAATGATGTCGAGATCCAACGTGCGTGGTCCCCAGCGT‑TCAGCTTAGCGGACGCGCCCTTGCTGCAATTCAATACGCTGTGCGTGACTGACTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTCAGTAATCGGGTTGATTTTGCGGCCCCAGCGGCGGGGTGc                                                                                                                                                                                                                               >  1:581271/1‑251 (MQ=255)                             cACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCCCGGCTGCGTTTACGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTGAAAc                                                                                                                                                                                                                       >  2:506080/1‑251 (MQ=255)                              ccGATTCTTCAGATCGTAGTTCGGAACGGTTAGGCGTTCAGTAGTTATCACTTCATTACCAAACAGCATGATGGCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCATCGGTGGGTTGCGGTAAAACg                                                                                                                                                                                                                      <  2:500574/250‑1 (MQ=255)                                gATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTCGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACga                                                                                                                                                                                                                     >  1:518291/1‑250 (MQ=255)                                          atCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTTAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTaa                                                                                                                                                                                                          <  1:407119/251‑1 (MQ=255)                                          accggaGTGCGGAACGTTCAGGCTTTCAGTATTTATCACTTCAT‑ACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAg                                                                                                                                                                                                         <  1:83363/246‑1 (MQ=255)                                             gTAGTGCGGAACGGACAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGa                                                                                                                                                                                                        >  2:274901/1‑250 (MQ=255)                                                   cGGAACGGTCAGGCGTTCAGTATTTCTCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGCCCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCACGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCACCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGc                                                                                                                                                                                                 >  2:254649/1‑251 (MQ=255)                                                       aCGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTGTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGAGTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                                                                                                                             <  2:277331/251‑1 (MQ=255)                                                        cGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCTTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCa                                                                                                                                                                                            <  2:357365/251‑1 (MQ=255)                                                         ggTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCACCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                                                                                                                             >  1:288943/1‑249 (MQ=255)                                                            cAGGCGTTCAGTTTTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGCGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTACGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGACTGTGGCTTTCAggg                                                                                                                                                                                         >  2:287736/1‑250 (MQ=255)                                                            cAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAggg                                                                                                                                                                                         >  1:195109/1‑250 (MQ=255)                                                                cGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGGGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATc                                                                                                                                                                                    >  1:232405/1‑251 (MQ=255)                                                                    cAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCTT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGcc                                                                                                                                                                                 <  1:461718/250‑1 (MQ=255)                                                                    cAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGcc                                                                                                                                                                                 <  1:466428/250‑1 (MQ=255)                                                                      gTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTa                                                                                                                                                                               >  2:7468/1‑250 (MQ=255)                                                                          ttatacctccactACCAAACACCATGATGTCGAGATGCAACGTTCGTGATCCCCATCGT‑TCAGCTTTGCGGACGCGGCCTTGCTGCACTTCAATACGCTGTGAGTGATTGAGTGGATCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCTTTTAATTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGAGCGGTAAAACGAAGAAACGGCAAGATTGTGGATTTCAGGGATATCGCCTAATGCt                                                                                                                                                                          >  2:422876/13‑251 (MQ=255)                                                                            aTCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTACGAATGTGGCTTTCAGGGATATCGCCTAATGCttt                                                                                                                                                                        >  2:13700/1‑251 (MQ=255)                                                                                 ttCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAggg                                                                                                                                                                   >  1:422392/1‑251 (MQ=255)                                                                                  tCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAgg                                                                                                                                                                    <  1:7468/249‑1 (MQ=255)                                                                                  tCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGCGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgc                                                                                                                                                                  >  1:890/1‑251 (MQ=255)                                                                                   cATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGGGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgca                                                                                                                                                                 <  2:460457/251‑1 (MQ=255)                                                                                    aTTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgcag                                                                                                                                                                <  2:86928/251‑1 (MQ=255)                                                                                      tACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGTTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGgcagc                                                                                                                                                               >  1:463266/1‑250 (MQ=255)                                                                                            aCAGCAAGATGGCCAGAACCAGCGTACGTAGTCCCCAG‑GTGTCAGCTTTGTGGACGCGGCCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGATGCAAGAGAGTTTGCCAGCGCCACGGCTGCGTTTAAGTAATAGGGTTGAGGTTGCGGCCCCATCGGAGGGGTGCGGTAAAACCAAGAAACGTTATGCATGTGGCTTTCAGGAAAATCGCCTTATGCTTTCAGGGCAGCATTGAcc                                                                                                                                                        <  2:388360/251‑1 (MQ=255)                                                                                              aGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCACTTCACTACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCACGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTg                                                                                                                                                      >  2:38322/1‑251 (MQ=255)                                                                                              aGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCCAGAGAGGTTTCCAGCGCCCCGGCTGCGTTTAAGTACTCGGGTTGATTTTGCGGCCCCCGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGCTATCGCCTAATGCTTTCAGGGCAGCATTGACCTg                                                                                                                                                      >  1:571235/1‑251 (MQ=255)                                                                                               gCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTCAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGc                                                                                                                                                     >  1:480846/1‑251 (MQ=255)                                                                                                cATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCt                                                                                                                                                    >  2:584274/1‑251 (MQ=255)                                                                                                 atgatgTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTc                                                                                                                                                   >  1:189510/1‑251 (MQ=255)                                                                                                 atgatgTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGACAGCATTGACCTGCTc                                                                                                                                                   >  1:393972/1‑251 (MQ=255)                                                                                                       tCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCgg                                                                                                                                             >  1:41291/1‑251 (MQ=255)                                                                                                               cAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCa                                                                                                                                     <  2:422392/251‑1 (MQ=255)                                                                                                                   gtgcgtgGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGAtt                                                                                                                                 <  1:58312/251‑1 (MQ=255)                                                                                                                       gtgGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATtgctg                                                                                                                             >  1:537292/1‑251 (MQ=255)                                                                                                                         ggTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGcc                                                                                                                           <  2:327356/251‑1 (MQ=255)                                                                                                                                aGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGTGAGGCCAGATTGCTGCCTTTGGc                                                                                                                     >  1:419536/1‑250 (MQ=255)                                                                                                                                  cGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAat                                                                                                                  >  1:16867/1‑251 (MQ=255)                                                                                                                                       cAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATAc                                                                                                              <  2:476250/251‑1 (MQ=255)                                                                                                                                       cAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTGCAGTGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAAtata                                                                                                               <  1:70404/250‑1 (MQ=255)                                                                                                                                        aGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTACAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGCTATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg                                                                                                             >  1:390103/1‑251 (MQ=255)                                                                                                                                        aGCTTGGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGGGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCAGTGACCTGCTTCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACg                                                                                                             <  2:546625/251‑1 (MQ=255)                                                                                                                                                ggACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTc                                                                                                     >  1:412499/1‑251 (MQ=255)                                                                                                                                                 gACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTc                                                                                                     >  1:71283/1‑250 (MQ=255)                                                                                                                                                 gACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTc                                                                                                     >  2:575742/1‑250 (MQ=255)                                                                                                                                                     cgACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCAACGGCTGCGTTTACGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAACACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCg                                                                                                >  1:568252/1‑251 (MQ=255)                                                                                                                                                            gctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGCGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCt                                                                                         >  1:340643/1‑251 (MQ=255)                                                                                                                                                               gcAATTCAATACGCTGGGTGTGATTGAGTAGCTCTTCAGGGGCAAGAGAGGGTTCCAGCGCCACGGCTTCGTTTAAGGAATCGGGTTGATTTTGCGGCCACAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTCTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTc                                                                                        <  2:587498/249‑1 (MQ=255)                                                                                                                                                                caattcaattCGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCacg                                                                                     <  2:71283/251‑1 (MQ=255)                                                                                                                                                                caagtCAATACGCTGTGTGGGAGTGAGTAGCGTTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGGTGATTTTGCTGCCCCAGCGGGGGGGGGCGGTAAAACGAAGAAACGGTAAGAATGTGGCGTTCAGTGATAGCGCCCAATGCTTTCAGGGCAGCAGTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCGATGTCAATATACTCCACGGGCTTGCGGTACCCTCacg                                                                                     <  2:571235/247‑1 (MQ=255)                                                                                                                                                                cAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCacg                                                                                     >  1:205921/1‑251 (MQ=255)                                                                                                                                                                                 tgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGc                                                                    >  1:306440/1‑251 (MQ=255)                                                                                                                                                                                  gtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCg                                                                   >  2:346707/1‑251 (MQ=255)                                                                                                                                                                                       tGAGTAGCTCTTCAGTTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATGTCGCCTAATGCTTTCAGGGCAGCAGTGTCCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATGTTCACGAATGCCAGGCGGTACCCTCGCGACGTGTTGCGCTTTTGCGTGga                                                               <  2:213455/250‑1 (MQ=255)                                                                                                                                                                                          gTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGgacgac                                                           <  1:472085/251‑1 (MQ=255)                                                                                                                                                                                          gTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTGCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGGACCCTCACGACGTGGTGCGCGTTTGCGTGgacgac                                                           <  2:890/251‑1 (MQ=255)                                                                                                                                                                                                 ttCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTATGCGTGGACGACGAGTacg                                                    >  2:450375/1‑251 (MQ=255)                                                                                                                                                                                                  tCAGGAGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGCGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTacg                                                    >  1:353636/1‑250 (MQ=255)                                                                                                                                                                                                   cagttGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGTGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCGCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTacgac                                                  <  1:250429/247‑1 (MQ=255)                                                                                                                                                                                                         aaaGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGGGGGCGTTTGCGTGGACGAAGAGTACGACGAcgcg                                            >  2:166328/2‑251 (MQ=255)                                                                                                                                                                                                                   ttcaaGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAAGATACGCCACTGTCATGCGGTACCCTCACGACGTGGTTCGCGTTTGCTTGGACGACGAGTACGACTACGCGGTGACGGtt                                   <  2:109067/246‑1 (MQ=255)                                                                                                                                                                                                                      cAGCGCCACGGCTGCTTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGttcttc                               <  2:232405/251‑1 (MQ=255)                                                                                                                                                                                                                                gCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCg                     <  2:41291/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAg                <  1:346707/251‑1 (MQ=255)                                                                                                                                                                                                                                      tttAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGCGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGc               >  1:142707/1‑251 (MQ=255)                                                                                                                                                                                                                                        tAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGAGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCAt             >  1:180175/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTCGGGTTGATTTTGCGGCCGCAGCGGTGGGGTGCGGTCAAAAGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGACTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGGGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGGTCTTCATCCAGCTCGTTGAACATCCCttt       >  2:414030/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCttt       >  1:391079/1‑251 (MQ=255)                                                                                                                                                                                                                                                cGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGCCGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCttttt     >  2:110799/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGGTCTTCATCCAGCTCGGTGAGCATCCCTTTTTGGt  <  2:448594/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     ATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGT‑TCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGGGCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGT  >  NC_000913/157338‑157829

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