Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,100,070	G→A	C144Y (TGT→TAT)	ycdY →	chaperone protein YcdY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,100,070	0	G	A	99.0%	299.0 / ‑4.2	97	C144Y (TGT→TAT)	ycdY	chaperone protein YcdY
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (45/51);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTA  >  NC_000913/1099825‑1100314                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     cGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGTATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCtggtatggtg                                                                                                                                                                                                                                                 >  1:201828/1‑251 (MQ=255)  gATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGc                                                                                                                                                                                                                                                >  1:662445/1‑251 (MQ=255)  gATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGc                                                                                                                                                                                                                                                >  1:777063/1‑251 (MQ=255)   attaaaaggtgttttttATCGGTGTTGAATGTGCTGGGCCGCCCTATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCg                                                                                                                                                                                                                                               <  2:860930/237‑1 (MQ=255)    ttACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGCGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGGATCTGTTACCCTGGTATGGTGCGt                                                                                                                                                                                                                                              >  1:557843/1‑251 (MQ=255)      aCAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGtt                                                                                                                                                                                                                                             <  1:298739/250‑1 (MQ=255)      aCAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTc                                                                                                                                                                                                                                            >  1:422786/1‑251 (MQ=255)                gttTATCGGCGATGAATGTGCTGTGCCGCCATATCGGAGCGCCTGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAAGGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGcca                                                                                                                                                                                                                                    >  2:543502/1‑247 (MQ=255)                      gggCGATGATTGTGCTGGGCCGCTATATCGTTGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACAGCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAg                                                                                                                                                                                                                             <  2:621789/249‑1 (MQ=255)                         cGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGccc                                                                                                                                                                                                                         >  2:279360/1‑251 (MQ=255)                          gATGAATGGGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCa                                                                                                                                                                                                                        <  2:372387/251‑1 (MQ=255)                            tGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCACCGGAAGATGAGAGCGAAGCACTGGAACCACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAATGTGGAGGcccttg                                                                                                                                                                                                                      >  2:221581/1‑248 (MQ=255)                                  ggcggtgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTGCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCa                                                                                                                                                                                                                <  1:221581/247‑1 (MQ=255)                                   gctgtgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCAc                                                                                                                                                                                                               <  1:234665/251‑1 (MQ=255)                                                     ggcgctggggTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑TGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTGCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCAt                                                                                                                                                                                             <  2:712968/244‑1 (MQ=255)                                                          aTGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCAc                                                                                                                                                                                        >  1:227673/1‑251 (MQ=255)                                                            gggTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCAcc                                                                                                                                                                                       <  1:230306/250‑1 (MQ=255)                                                              gTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCt                                                                                                                                                                                    <  1:216895/251‑1 (MQ=255)                                                                tGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGACTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGCGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAACGTGGAGGCCCATGAAACCACGCCTTTATGCAGCACCATGGCACCCCTa                                                                                                                                                                                   >  2:844437/1‑250 (MQ=255)                                                                   gggCGCGACGGAAGCGGAAGCGCTCTCTTTTCTTTCTGAGCGAGGGATGCCCTT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTTTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAAccc                                                                                                                                                                               <  2:185111/251‑1 (MQ=255)                                                                      cgcgACGGAAGCGGAAGTGCGCTCTTTTCTTTCTGAGCGAGGGAGGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGATGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCcgcg                                                                                                                                                                            <  1:381731/251‑1 (MQ=255)                                                                         gACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCCTT‑AGCGGAAACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCCACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGGATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGGGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCAACGCTAACCCGCGATg                                                                                                                                                                         >  2:780949/1‑251 (MQ=255)                                                                            ggaagcggaagTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGcc                                                                                                                                                                       >  1:74747/1‑250 (MQ=255)                                                                              aagcggagtTGCCCGCTTTTCTTTCTGCGCGAGGGATGC‑ATTGAGCGGATCCGCCAGCCGATCACATCGGCACATTGCGGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACGGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                    <  2:510720/251‑1 (MQ=255)                                                                              aagcggaagTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCAtt                                                                                                                                                                    <  1:153366/251‑1 (MQ=255)                                                                                    aagTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                               >  1:875261/1‑250 (MQ=255)                                                                                    aagTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                               >  1:10600/1‑250 (MQ=255)                                                                                      gTGCGCGTTTGTCTTTTGGAGCGATGGATGCCATT‑GGCGGTTGCGCCATGCGATCACATCGGCACAGTGCTGCTCGCAGCTTCCTGGCTGTGAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGGGAGTAGCTGTTACCCGGGTATGGTTCGTTCCTTGGCAAAGTGGAGGCCTATGCAACCACGCCTTTCTGTCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAAt                                                                                                                                                            <  2:249484/251‑1 (MQ=255)                                                                                           cgcTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGTCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTggg                                                                                                                                                       >  1:407436/1‑251 (MQ=255)                                                                                               tgtcttccttggggAGGGATGCTATT‑AGCGGTTACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACg                                                                                                                                                    <  2:865016/237‑1 (MQ=255)                                                                                                 gtttactgcgggagggagtccATT‑AGTGGATACGCCAGCCGATCACTGCGGCAAATTGCTGCTCGCAGCTTTCGTGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACGGGAAACACTGTTCAGTGAGGATCTGTTACCCTGGTATGGGGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCGTTCTGGCGCACCATGGCACCGCTAACCCGCGATTCCATTAGTGCAATGTGGGActag                                                                                                                                                  <  2:639687/231‑4 (MQ=255)                                                                                                    ttctGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTg                                                                                                                                               <  1:384067/250‑1 (MQ=255)                                                                                                     tctGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTgg                                                                                                                                              <  1:543502/250‑1 (MQ=255)                                                                                                         aGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGAATCTGTTACCCTGGTATGGTACGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCTCACCCAGGCACCGCTAACCCGCGATGCCATTATCGCAATGTGGGACGAGCTGgaaga                                                                                                                                          >  2:552665/1‑250 (MQ=255)                                                                                                            gagtgATGCCATT‑AGCGGATACGCCAGCCGCTCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGAATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCCCCGCTAACCCGCGCTGCCATTAGTGCAATGAgggccggactgcaagaagct                                                                                                                                      >  2:279248/5‑234 (MQ=255)                                                                                                            gAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGAt                                                                                                                                      <  1:133730/251‑1 (MQ=255)                                                                                                             aGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGcc                                                                                                                                                                                                                          <  1:797267/166‑1 (MQ=255)                                                                                                             aGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGGGGAGGcc                                                                                                                                                                                                                          >  2:797267/1‑166 (MQ=255)                                                                                                                 aTGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTg                                                                                                                                  >  1:10893/1‑250 (MQ=255)                                                                                                                   gCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCCTTAGt                                                                                                                                                                 >  2:109756/1‑217 (MQ=255)                                                                                                                   gCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGt                                                                                                                                                                 <  1:109756/217‑1 (MQ=255)                                                                                                                          aGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAg                                                                                                                         <  1:349413/251‑1 (MQ=255)                                                                                                                          aGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAg                                                                                                                                                  <  2:563012/226‑1 (MQ=255)                                                                                                                          aGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAg                                                                                                                                                  >  1:563012/1‑226 (MQ=255)                                                                                                                                 aCGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGGGCTGGAAGAAGATTCTGAAGAATAAGTGTGATc                                                                                                                  >  2:16557/1‑251 (MQ=255)                                                                                                                                        ccGATCACATCGGCATATAGCTGCTCGCAGCTGCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACGTGAAACACCGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCa                                                                                                           <  2:86173/251‑1 (MQ=255)                                                                                                                                        ccGATCACATCGGCACAGTGCTGCTCGCAGCTTCCGGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGCGCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCGCCGCTAACCCGTGATGCCAGTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTTTGATCTACGTCa                                                                                                           <  2:662445/251‑1 (MQ=255)                                                                                                                                         cGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTAGCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCAc                                                                                                          <  2:355691/251‑1 (MQ=255)                                                                                                                                          gATCACATCGGCACATTGCTGCTCGCAGCTTCCTGTCTGGAAGATCAGTCAACGGCAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGCGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCAATGCAACCACGACTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGCAGAATAAGTGTGATCTACGTCACt                                                                                                         >  1:811348/1‑251 (MQ=255)                                                                                                                                           aTCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTc                                                                                                        >  1:177120/1‑251 (MQ=255)                                                                                                                                             cacaTCGGCTCATTGCTGCTCGCAGCTTCCTGGCTGGAAGAGCAGGCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCAATAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCa                                                                                                       <  2:6160/250‑1 (MQ=255)                                                                                                                                             cacaTCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCa                                                                                                       <  1:795700/250‑1 (MQ=255)                                                                                                                                                aTCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAAc                                                                                                   >  1:404169/1‑251 (MQ=255)                                                                                                                                                aTCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAAc                                                                                                   >  1:240085/1‑251 (MQ=255)                                                                                                                                                  cGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACt                                                                                                  <  1:801070/250‑1 (MQ=255)                                                                                                                                                    gCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGACGAATAAGTGTGATCTACGTCACTCATAACTGc                                                                                                >  2:21010/1‑250 (MQ=255)                                                                                                                                                       caTTGCTGCTCGCAGCTGCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTGCGTCACTCATAACTGCAACg                                                                                            <  2:393173/251‑1 (MQ=255)                                                                                                                                                              gctCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGTAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTTCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCTATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAtt                                                                                     <  1:86195/251‑1 (MQ=255)                                                                                                                                                                   cAGCTTCCTGGCTGGAAGATCAGTCAACGGCAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGAATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTTGAGGCCCATGCAAACACGCATTTCTGGCGTACCATGTCACCGCTAACCCGCGATGCCATTATACCCATGTGGGGCGAGCCGGAAGAAGATTCAGAAGAATACGTGTGACCTACGACACTCATACATGCAACAGATAATttgtct                                                                                >  2:163492/1‑249 (MQ=255)                                                                                                                                                                           tggctggAAGCTCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATACTTTGTTGTTGCATa                                                                         >  1:548688/1‑250 (MQ=255)                                                                                                                                                                            ggctggAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATaaa                                                                       <  1:554284/251‑1 (MQ=255)                                                                                                                                                                               tggAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATaaaagg                                                                    >  2:732622/1‑249 (MQ=255)                                                                                                                                                                                gggagatacgtaaaCGGAAGATGGGTGCGTAGTACTGGAAACACTGTACAGTGTGTGTCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAGGTGGCGGCCAATGCAACCACGCCTTGCCGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGGAGAATAAGTGCGATCTACGCCACTCATAACTGCAACGGATACTTTGTTGTTGCATAAAAtgt                                                                   >  2:489831/13‑251 (MQ=255)                                                                                                                                                                                 gAAGGTCAGTCAACGGAAGATTAGAGCGAAGCACGGGAAACTCTGTTCAGTGAGTAACTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCGTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTCCGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtgtg                                                                  <  2:252754/251‑1 (MQ=255)                                                                                                                                                                                   agagtaGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCTCCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtgtgtg                                                                <  2:279950/246‑1 (MQ=255)                                                                                                                                                                                   aGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGAGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCAAGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtgtgtg                                                                >  2:863578/1‑251 (MQ=255)                                                                                                                                                                                    gATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGAGGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGGCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGc                                                               <  2:683171/251‑1 (MQ=255)                                                                                                                                                                                      gcAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTc                                                             <  2:207413/250‑1 (MQ=255)                                                                                                                                                                                       cAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTc                                                             <  1:206385/250‑1 (MQ=255)                                                                                                                                                                                       cAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCg                                                            >  1:741446/1‑251 (MQ=255)                                                                                                                                                                                            aaCGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGAtctc                                                       >  1:483705/1‑251 (MQ=255)                                                                                                                                                                                               ggAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCtcatt                                                    >  1:504491/1‑251 (MQ=255)                                                                                                                                                                                                 aaGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCtcattca                                                  <  1:655083/251‑1 (MQ=255)                                                                                                                                                                                                 aaGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCtcattca                                                  <  1:824808/251‑1 (MQ=255)                                                                                                                                                                                                  aGATGACAGCGAAGCACGGGAAACACGGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAGTTGTTGTTGCATAAAATGTGTGCTCGATCtcattcat                                                 <  2:59319/251‑1 (MQ=255)                                                                                                                                                                                                    aTGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATgg                                               >  2:258205/1‑251 (MQ=255)                                                                                                                                                                                                       agagCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcg                                            <  1:599468/251‑1 (MQ=255)                                                                                                                                                                                                       agagCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcg                                            >  1:527786/1‑251 (MQ=255)                                                                                                                                                                                                             aaGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGAATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGAGGGACGAGCTGGGAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGtttt                                      >  2:363123/1‑251 (MQ=255)                                                                                                                                                                                                              aGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTc                                     <  1:21010/251‑1 (MQ=255)                                                                                                                                                                                                                cACTGGGAACCCTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGATCCTTGGCAAAGTGGAGGCCCATACAACCACGCCTTTCTGGCGAACCATGGCACCGCTAACACGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCACCGGATAATTTGTTGTTGCGTAAAATGTGTGCTCGATCTAATTCATGGCCGCGTTTTCTg                                   >  2:113774/1‑251 (MQ=255)                                                                                                                                                                                                                   tGGAAACACTGTACAGTGAGTATCGGTTACCCTGGGATGGAGCGTTCCTTGGCCAAGTGGAGGCCCATGCAACCACGCCTTTCTGTCGCACCAGGGCACCGCTAACCCGCGCTGCCCTTTGTGCAATGTGGGACGGGCTGGCGGGAGATTCTGAAGAATAAGTGTGATCTACGTCGCTCATACCTGCACCGGCTACTTTACTGTTGCTAAAAATGTGTGCTCGATCTCCTTCATGTACGCGTTTTCTGCtt                                >  2:299327/1‑250 (MQ=255)                                                                                                                                                                                                                        acacTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATg                           >  1:710257/1‑251 (MQ=255)                                                                                                                                                                                                                         cacTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                <  1:270871/165‑1 (MQ=255)                                                                                                                                                                                                                         cacTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTa                                                                                                                >  2:270871/1‑165 (MQ=255)                                                                                                                                                                                                                         cacTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgc                          <  1:589403/251‑1 (MQ=255)                                                                                                                                                                                                                           cTGTTCAGTGAGTTTCTGTTAACCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTACCCCGCGATGCCATTAGTGCAATGAGGGACGAGCTGGAAGAAGATTCTGAAGCATACGTGTGATCTACGTCACTCATAACTGCAACGGATAAttt                                                                                    >  2:282938/1‑191 (MQ=255)                                                                                                                                                                                                                           cTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCGTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAttt                                                                                    <  1:282938/191‑1 (MQ=255)                                                                                                                                                                                                                               tCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGc                     >  1:499937/1‑250 (MQ=255)                                                                                                                                                                                                                                 aGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATg                  <  1:279360/251‑1 (MQ=255)                                                                                                                                                                                                                                   tGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGa                 >  1:470764/1‑250 (MQ=255)                                                                                                                                                                                                                                    gggTATCTGTTATCCGGGTATTGTGCGTTCCGTTGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAGGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGaa                <  2:198655/248‑1 (MQ=255)                                                                                                                                                                                                                                           tGTTACCCTGGTATGGTGCGTTCCTAGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtc        >  1:624915/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACt          <  1:780949/248‑1 (MQ=255)                                                                                                                                                                                                                                            gTTACACTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACt          <  2:158747/248‑1 (MQ=255)                                                                                                                                                                                                                                              tACCCTGGTATGGTGCGTTCCTTGTCAAAGTGGAGGCCCATGCAACCACGCCTGTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctc      <  2:134957/250‑1 (MQ=255)                                                                                                                                                                                                                                              tACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTGTGCTATCATGCGCGGCATGAACATACTtctc      <  1:238766/250‑1 (MQ=255)                                                                                                                                                                                                                                                 ccTGGTATGGTGCGTTCCTTGGACAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTa  <  2:201828/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATT‑AGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTA  >  NC_000913/1099825‑1100314

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: ‑