Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,985,672	G→A	N112N (AAC→AAT)	aslA ←	putative sulfatase AslA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,985,672	0	G	A	98.7%	234.0 / NA	78	N112N (AAC→AAT)	aslA	putative sulfatase AslA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (35/42);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTG  >  NC_000913/3985431‑3985906                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           cccccATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCAcc                                                                                                                                                                                                                                   <  1:333300/251‑1 (MQ=255)       tatGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCAcg                                                                                                                                                                                                                             >  1:20192/1‑251 (MQ=255)        atGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCAcgc                                                                                                                                                                                                                            >  2:266196/1‑251 (MQ=255)           ccATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGAGCACCAACTGCGGCAGCGTGGTTCACCCTTGCAGCCCGCCCGGTTGGCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGACAACGGCGGCGATATCTGTGGA‑AGTGTTACCCACCGCCAcgccgc                                                                                                                                                                                                                         >  1:355473/1‑251 (MQ=255)             aTTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCa                                                                                                                                                                                                                       <  1:569073/251‑1 (MQ=255)                ttCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCCACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGCATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACCCCGCCACCg                                                                                                                                                                                                                    >  1:63372/1‑251 (MQ=255)                  ccGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGTTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGt                                                                                                                                                                                                                   <  1:605084/250‑1 (MQ=255)                   cGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTa                                                                                                                                                                                                                 <  1:160878/251‑1 (MQ=255)                       ggCCTGAGGGACGTAGCCCTGACCGTGCAGCAACTGAGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGAGGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGCGAAGAGCTTGGTTGAGAATACGCCGAAGATAAAATCAGCCCCTGGCTGGCAACGGCGTTGATATCTGGTGT‑AGGATTACCCACCGCTACGCCCCCACTGTTAAAAc                                                                                                                                                                                                             >  2:614906/1‑251 (MQ=255)                                 aCGTAGCCCTGATGGGGCAGCAGGTGCGGCAGCGTGGTTAACCCTTGAAGCCCGCCGGGTTGCTCGTACATTGGCGGCATCAGAATGCCGTGGGGGATGGAGTATTGTACGGTGAGAAGGGTGGCGCGGGTTGGGGAAGAGGTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCGGGCAACGGGGTTGATATCTGGTGT‑CGTATTACCCCCCGCCACGCCGCCACCGTTAAAATCGACGtccat                                                                                                                                                                                                   <  2:15723/251‑1 (MQ=255)                                 aCGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGtccat                                                                                                                                                                                                   >  1:379982/1‑251 (MQ=255)                                    tAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGtccancc                                                                                                                                                                                                 >  2:587366/1‑250 (MQ=255)                                        ccTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATAAGCCGAAGTTAAAATCAGCCCCTGGCGGGCAACGGCTTCGATATCTGGTGT‑AGGATTACCCCCCGCCACGCAGCCACCGTTAAAACCGACGTCGATCCAGccc                                                                                                                                                                                            >  2:347968/1‑251 (MQ=255)                                            aTCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACAt                                                                                                                                                                                        <  1:116464/251‑1 (MQ=255)                                                      aaCTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCa                                                                                                                                                                               >  1:213376/1‑250 (MQ=255)                                                      aaCTGCGACAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCGTCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCAACGTTCAAACCGACGTCCATCAAGCCCACATCGTCCAGCa                                                                                                                                                                               >  2:75363/1‑250 (MQ=255)                                                          gcggcggcgTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGa                                                                                                                                                                            <  1:331730/249‑1 (MQ=255)                                                          gCGGCAGCGTGGTTAACCCCTGCAGCCCGCCCGGTTGCCGTTGCATTGGCGGCATCAGAAGGCGGTGGTGGATGGATTATTTTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTGGGTTAAGAATACGCCGAGGTTAAAATCAGCCCCTGGCTGGCAACGTCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGaaa                                                                                                                                                                          <  2:146521/251‑1 (MQ=255)                                                          gCGGCAGCGGGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGACTAAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATCCGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCAACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGaa                                                                                                                                                                           <  1:188796/250‑1 (MQ=255)                                                            ggCAGCGTGGTTAACCCTTGCAGCCCGACCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGATGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGAGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATTCAGCCCACATCGTCCAGCAAGAAaac                                                                                                                                                                        >  1:473750/1‑251 (MQ=255)                                                             gCAGCGGGGTTAACCCTTGCAGCCCGCCGGGGTGCCCGTAAATGGGCGGCATCAGAAGGCGGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGTTTGAGAATACGCCGAAGTTAAAATCAGCCCCGGGCTGGCAACGGCGTCGATATCTGGGGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCAGCCAGCCCACATCGTCCAGCAAGAAaaca                                                                                                                                                                       <  2:205354/251‑1 (MQ=255)                                                              cAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGCTGGAGTATAGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAAAGGCGCCGATATCTGGTGT‑AGGATTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAaaca                                                                                                                                                                       >  1:47879/1‑250 (MQ=255)                                                              cAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAaaca                                                                                                                                                                       >  1:276284/1‑250 (MQ=255)                                                               ggggtgtttAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCGGGGGTGGATGGAGTATTGTCCGGTGAAAATGGGGGCGCGGGTTGGGGAAGAGCTTGGTGGAGAATACGCCGAAGTGAAAAGCAGCCCCGGTCTGGCAACGGCGTCGATATCGGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCTCACATCGTCCAGCAAGAAaacaac                                                                                                                                                                     <  2:643246/244‑1 (MQ=255)                                                                     gTTAACCCTTGCAGCCCGCCCGGTTGCCCGGACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGGGGCGCGGGTGGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCACCCCCTGCCTGGCAACGGCGTCGATATCTGTTGT‑AGTATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACAtt                                                                                                                                                               <  2:97368/251‑1 (MQ=255)                                                                       tAACCCTTGCAGCCCGCCCGGTTGCCCGTACATGGGCGGCATAAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTc                                                                                                                                                              <  1:152945/250‑1 (MQ=255)                                                                        aaCCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCGCATTCg                                                                                                                                                             <  1:182851/250‑1 (MQ=255)                                                                        aaCCCTTGCAGCCCGCCCGGTTGCCCGAACATTGGCGGCATCAGAATGCCGTTGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACACCGAAGTTAAAATCAGCCCCGGGCTGGCAACGGCGTCGATATCTGGAGT‑AGGATTACCCACCGCCACGCCGCCACCGTTCAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAAACACATTc                                                                                                                                                              >  1:356088/1‑249 (MQ=255)                                                                         aCCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGt                                                                                                                                                           >  1:42973/1‑251 (MQ=255)                                                                          cctggcatccctcccgggTTGCC‑GTTCATGGGCGGCATCAGAGGGACGTGGTGGATGGAGTTTTGTCCGGTGAAAATGGTGGCGCGGGTTGGGGAGGAGCTTTGTTAAGAATAGACCGAATTTAAAACCACCCCCTGGCTGGCAGCGGCGTCTATTTCTG‑TGTTAGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCTTCCAGCAAGAAAACAACCACCTTCTGttt                                                                                                                                                         <  2:330843/235‑1 (MQ=255)                                                                           ccTTGAAGCCGGCGCGTTTGCCCGTACATTGGCGGCATCAGAAGGCCGTGGTGGATGGAGTATTGGCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGGTAATACGCCGAATGTAAAGTCAGCCCCGGGCTGGCTACGGCGTTGATATCTGGTGT‑AGCATTACCCACCTCCACGCCGCCACCGTTAAAGCCTACGTCCATCCATCCCACATCGTCCAGCAAGAAAACAACCACATTCGGttt                                                                                                                                                         <  1:601963/251‑1 (MQ=255)                                                                            cTTGCAGCCCGCCAGGTTGCCCGTACCTTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCCCCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGnttt                                                                                                                                                        >  2:593335/1‑250 (MQ=255)                                                                              tGCAGCCCGCCCGGTTGCCTTTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCAGGTGGGAATGGTGGGGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTAAAAATCAGCCCCTGGCGGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATGCGGTTTCt                                                                                                                                                       <  2:63372/250‑1 (MQ=255)                                                                              tGCAGCCCGCCCGGTTGCCTGTACATTGGGGGCATCAGAATGCGGGGGTGGATGGAGTAGTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAACCAGCCCCTGGCTGGCAACGCCGCCGATATCTGGTGT‑AGAATCACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCt                                                                                                                                                       <  1:254781/250‑1 (MQ=255)                                                                               gCAGCCCGCCCGGTTGCCCGTACTTTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTTGCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTGCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTa                                                                                                                                                     <  1:320285/251‑1 (MQ=255)                                                                               ccAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCTTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGTTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTa                                                                                                                                                     <  1:65932/250‑1 (MQ=255)                                                                                  gcccgcccGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCg                                                                                                                                                  <  1:348118/251‑1 (MQ=255)                                                                                     cgcccgGTGGCCCGTACATGGGCGGCATCAGAATGCCGGGGTGGAGGGAGTATTGTCCGGTGAGAAGGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGtt                                                                                                                                               <  2:20192/251‑1 (MQ=255)                                                                                        ccgGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGGGGGTTGGGGAAGAGCTTGTTTGATAATACGCCGAAGTTAAAATCAGCCCCTGCCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGtttttt                                                                                                                                            <  1:516434/251‑1 (MQ=255)                                                                                        ccgGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGCGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGttttt                                                                                                                                            >  1:447025/1‑251 (MQ=255)                                                                                          gtttGCCCGTACATTGGGGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTGGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGttttttt                                                                                                                                          <  1:362656/249‑1 (MQ=255)                                                                                             tGCCCGTACATGGGCGGCATCAGAATGCCGGGGGGGATGGAGTATTGTCCGGTGAGAAGGCTGGCGCGGGTTGGGGAAGAGCTTTGTTGAGAATACGCCGAAGTTAAAATCAGCCCCGGGCTGGCACCGGCGTCGATATCTGGTGT‑AGTATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGttttttgct                                                                                                                                        <  2:191763/250‑4 (MQ=255)                                                                                                      aTTGGCGGCATCAGAATGCCGGGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGc                                                                                                                              <  1:39065/251‑1 (MQ=255)                                                                                                                 cAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAACCAGCCCCTGGCTGGCACCGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTctgc                                                                                                                    <  1:453070/250‑1 (MQ=255)                                                                                                                  aGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTctgctg                                                                                                                  >  1:545584/1‑251 (MQ=255)                                                                                                                    aaTGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGCAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCCCCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTggg                                                                                                                >  2:661041/1‑251 (MQ=255)                                                                                                                        cggtggtgGATGGGGTATTGTCGGGTGAGAATGGAGGGGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAATTTAAAAACAGCCCCGGGCTGGAAACGGCGTCGATATCTGGTTT‑AGGATTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGGCCACATCGTCCAGCAAGAAAACAACCACATTCGGTGTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTc                                                                                                            <  2:364778/249‑1 (MQ=255)                                                                                                                          gtggtgGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCtt                                                                                                          <  1:266196/251‑1 (MQ=255)                                                                                                                          gtggtgGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCtt                                                                                                          >  2:183920/1‑251 (MQ=255)                                                                                                                            gttgtttggtggatatgccGGTGAGAATGGTGGCTCGGGTTGGTGGAGAGCTTGGTTGAGAATACCCCGATTTTAAAATCGCCCCCCTGCTGGCCCCGGCTTCGATACCGGGTGT‑CGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATGCGGTTTCTTACCGGTTTTTTTCTCAACTTCTGCCGGCTTCTCCTGGGTTTCttt                                                                                                         <  2:383984/233‑1 (MQ=255)                                                                                                                                gatggaGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCt                                                                                                    >  1:449572/1‑251 (MQ=255)                                                                                                                                   ggaGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTgcg                                                                                                 >  2:174729/1‑251 (MQ=255)                                                                                                                                     aGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGcc                                                                                               <  1:183920/251‑1 (MQ=255)                                                                                                                                      gTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCg                                                                                              <  1:105135/251‑1 (MQ=255)                                                                                                                                            tCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGACGTTAAAATCAGCCCCTGGCTGGCAACGGCGTAGATATCTGGTGA‑AGGATTAACCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGATTCTGCTGGGTTTCTTTATACTGCGCCGGAtgct                                                                                        >  2:621673/1‑251 (MQ=255)                                                                                                                                                   agaATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGAGATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGACCACATCGTCCAGCAAGAAAACAACAACATTCGGTTTCTTAGCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGCTGCTGCATTa                                                                                  >  2:410274/1‑250 (MQ=255)                                                                                                                                                        ggtggCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCTTTACTGGcc                                                                            >  2:524310/1‑250 (MQ=255)                                                                                                                                                          gggCGCGGGTGGGGGAATGGCTTGGTTGAGAATACGCCGAAGTTAAAATCGTCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGTTGCATTACTGGcatc                                                                          <  2:655741/250‑1 (MQ=255)                                                                                                                                                             cgcgGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCAtt                                                                                   <  2:175610/239‑1 (MQ=255)                                                                                                                                                             cgcgGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCAtt                                                                                   >  1:175610/1‑239 (MQ=255)                                                                                                                                                                  gTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTc                                                                  >  1:454165/1‑251 (MQ=255)                                                                                                                                                                     ggggAGGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGTTATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTCCTGGCATCATATTGTCGGc                                                               <  2:493197/251‑1 (MQ=255)                                                                                                                                                                         aaGAGCTAGGATGAGAATACGACGAAGTTAAAATCAGCCCCTGGCAGGCAACGGAGTCGATATCTGGTGT‑AGGATTACACACAGCAACGCCGCCACCGTTAAAACCGACGTCCATCCAGACCACATCGTCCAGCCAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTAAAGTTCTGCAAGCTGCTGCTGGGGTTCTTTATCCTGCGCCGGATGCTGCATTACCTGCATATTATTGTCGGCNAt                                                            >  2:608854/1‑250 (MQ=255)                                                                                                                                                                          tgagattGGATGAGAAGACGCCGAAGTTAAAATCAGCCCCAGGCTGGAAACGGCGTCAATATCTGATGT‑AGGATTACCCACCGCCACGCCGACACAGTTAAATCCGACGTACATCTAGCACACAACGTCAAGCAAGAAACTAACCACATTCGGTTTCTTACCGGTTTTTTTCTCCAGTTCCGTCAGCTTCTGCTGGGATTCTTTATCCTGCGCAGGATGATGCATTACTGGCATCATATGGTCGGAAAtcg                                                          >  2:249756/6‑249 (MQ=255)                                                                                                                                                                                tGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCg                                                    >  2:390077/1‑251 (MQ=255)                                                                                                                                                                                     gagaATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCTCATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTGCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTAGCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCgg                                                <  2:65360/250‑1 (MQ=255)                                                                                                                                                                                       gaATATGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGCCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACGGACGTCCATCCAGCCCACATCGTCCAGCAAGCAAACAACCACATTCGGTTTCTTACCGGTTTATTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGTCTCATATTGTCGGCAATAGCGGTCGCTGGttt                                             >  2:332718/1‑251 (MQ=255)                                                                                                                                                                                                gAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATa                                    <  1:587366/251‑1 (MQ=255)                                                                                                                                                                                                  aGTTAAAACCAGCCCCGGGCTGGCCACGGCGTCGATATCTGGTGT‑AGGATTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTATCCAGATACt                                  <  2:149748/251‑1 (MQ=255)                                                                                                                                                                                                          tCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTggg                          >  2:500519/1‑251 (MQ=255)                                                                                                                                                                                                          tCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTggg                          >  2:543792/1‑251 (MQ=255)                                                                                                                                                                                                           cAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTAAAGTTCTGCCAGCTTCTGCTGGGTTTATTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGt                         >  1:31382/1‑251 (MQ=255)                                                                                                                                                                                                              ccccTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGAt                      <  1:174729/251‑1 (MQ=255)                                                                                                                                                                                                                ccgggctggcAACGGCGTCGATATCTGGTGT‑AGGATGACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCg                    <  2:566717/248‑1 (MQ=255)                                                                                                                                                                                                                       ggcACCGGCGTCGATATCTGTTTT‑AGGATTACCCACCTCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATcc               <  1:75363/249‑1 (MQ=255)                                                                                                                                                                                                                             ggCGTCGATCTCTGGTGT‑AGGATCACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAACACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGc       >  2:221411/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGATATCTGGTGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAACACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCCGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCACAGCCTTtg  >  1:24256/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGATATCTGGGGT‑AGGATTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGACCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTtg  <  1:570242/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           CCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGT‑AGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTG  >  NC_000913/3985431‑3985906

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