Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,969,982	C→T	A335T (GCG→ACG)	tap ←	methyl‑accepting chemotaxis protein Tap

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,969,982	0	C	T	98.8%	281.2 / NA	84	A335T (GCG→ACG)	tap	methyl‑accepting chemotaxis protein Tap
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (45/37);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTT  >  NC_000913/1969750‑1970228                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          tctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCAGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGc                                                                                                                                                                                                                                         >  1:335890/1‑251 (MQ=255)  ctTTTGCTGCCTGCGCGCTACGGCTGGCATGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGc                                                                                                                                                                                                                                         >  1:153690/1‑250 (MQ=255)    tttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAtt                                                                                                                                                                                                                                      >  1:280527/1‑251 (MQ=255)     ttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAttt                                                                                                                                                                                                                                     >  2:119583/1‑251 (MQ=255)     ttGCGGCCCGCGCGCTACGGCTGGCAAG‑TTGCCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAtt                                                                                                                                                                                                                                      <  1:189894/250‑1 (MQ=255)      tgcggcctgcGCTCGACGGCTGGCAAGATGGCTC‑ACTCCACCTGCCCCTACCGCAAAACCACGCCCCTGCTCGCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAttt                                                                                                                                                                                                                                     <  1:521265/246‑1 (MQ=255)          gcttgcGCGCTACGGCGGCCAAGATGGCCA‑ACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCAtttttt                                                                                                                                                                                                                                  <  1:312678/246‑1 (MQ=255)           cctgcGCGCTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGcc                                                                                                                                                                                                                               >  2:138706/1‑251 (MQ=255)                 cgcTACGAATGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCAACTGCCGCCTTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCCATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTAATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCATTTcc                                                                                                                                                                                                                         >  2:57646/1‑251 (MQ=255)                   cTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGCACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCg                                                                                                                                                                                                                        >  2:447286/1‑250 (MQ=255)                     ccgctgggaaataTTGCGC‑TCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCGCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGGCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACg                                                                                                                                                                                                                     <  1:46873/239‑1 (MQ=255)                        gCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGTTCCCCTGCATGGGGTGGGTCAGGGAACGAACCTGACCACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCGGACGCCt                                                                                                                                                                                                                  >  1:203407/1‑251 (MQ=255)                            gCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCg                                                                                                                                                                                                              >  1:312161/1‑251 (MQ=255)                               aGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGc                                                                                                                                                                                                           >  2:216666/1‑251 (MQ=255)                                    gcgc‑CCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTAt                                                                                                                                                                                                      <  1:518303/251‑1 (MQ=255)                                      gc‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGACTATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGAGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCg                                                                                                                                                                                                    >  1:270914/1‑251 (MQ=255)                                      ccccttcccccTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATc                                                                                                                                                                                                     <  2:583474/242‑1 (MQ=255)                                       cccctccACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCg                                                                                                                                                                                                    <  1:558117/245‑1 (MQ=255)                                           ttCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGt                                                                                                                                                                                                >  2:439626/1‑251 (MQ=255)                                                  gCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGccc                                                                                                                                                                                         >  1:402798/1‑251 (MQ=255)                                                        cccgaaaaaCCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGt                                                                                                                                                                                   <  2:147327/246‑1 (MQ=255)                                                            cAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCg                                                                                                                                                                               >  1:87604/1‑251 (MQ=255)                                                                       cccTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTc                                                                                                                                                                    <  2:139627/251‑1 (MQ=255)                                                                           cctctCCGGCCCTAGCCGCTTCCTCTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATa                                                                                                                                                                <  2:159979/250‑1 (MQ=255)                                                                                ccGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGGGGCGGTTAATTGCTCTATACt                                                                                                                                                              <  1:130571/248‑1 (MQ=255)                                                                                    cgcgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATAATGGCGGcc                                                                                                                                                       >  2:434771/1‑251 (MQ=255)                                                                                     gcgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCg                                                                                                                                                      <  1:216666/251‑1 (MQ=255)                                                                                      cgAGCCGCTTCCACTGCAGCATTCAGGGCCAGACTATTGGTCTGGAAAGCAATTCCGCCGATAACGCTGCTAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCgg                                                                                                                                                     >  2:39899/1‑250 (MQ=255)                                                                                       gAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGtt                                                                                                                                                    >  1:14122/1‑251 (MQ=255)                                                                                            gCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCGCCATACTGGCGGCCGTTTGT‑Gc                                                                                                                                              <  2:402798/251‑1 (MQ=255)                                                                                            gCTTCCACTGCCGCATTCAGGGCAAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTTGCGGCCGTTTTT‑Gc                                                                                                                                              >  2:417514/1‑251 (MQ=255)                                                                                               tCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGCTAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGGTAATTGCTCCCTACTGGGGGCCGtgtgt‑ggagg                                                                                                                                           >  2:470386/1‑247 (MQ=255)                                                                                               tCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTGAATTGCTCCATACTGGCGGCCGTTTGT‑Gc                                                                                                                                              <  1:470386/248‑1 (MQ=255)                                                                                                      ccGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCTTCGATAACGCTGAGAAGGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCAGGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTTCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCGACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGG‑GCCAGCGATgca                                                                                                                                    <  2:398201/251‑2 (MQ=255)                                                                                                            acAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAGAACTCTTATAATGGCGCCAATGTTGTGCGAGCTGGTGGCGATCTCCTGCATTGGGTGAGTCATGGTACTGTCCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTGCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTCCCGTGGCGGTTAATGGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGTctg                                                                                                                               <  2:308391/249‑1 (MQ=255)                                                                                                              aGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATATCGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCctgctg                                                                                                                            <  2:445009/251‑1 (MQ=255)                                                                                                               gggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGGTGCCGCCTGCTGt                                                                                                                           >  2:104603/1‑251 (MQ=255)                                                                                                                ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTTCGTTATCGGCGTTTTGcc                                                                                                                                                                                          <  1:508803/188‑1 (MQ=255)                                                                                                                ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGcc                                                                                                                                                                                          >  2:508803/1‑188 (MQ=255)                                                                                                                 gCCAGAATAATGGTATGGAAAGCAATGCCTTCTATAACTCTTATAATGGCGCGAATTTTGTGCGAGCTGGTAGCGAGCTGCTGTATGGTGTATGTCAGGGTACTGAGCGGACAACAGCCTGTCTGCGCCGTTGTCGCATTATTTTGTGCCAGTTCCAACGCCTGTCGTTCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTTCATACTGGCTGCCGTTTGT‑GCCAGCGATGCTGCCTTCTGttt                                                                                                                         <  2:11122/251‑2 (MQ=255)                                                                                                                   cAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGTTTGT‑GCCAGCGGTGCCGCCTGCTGTTCgg                                                                                                                       >  2:457582/1‑251 (MQ=255)                                                                                                                             gTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGGTCGGTACGACTTgg                                                                                                             >  1:388163/1‑250 (MQ=255)                                                                                                                               cTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACtttgg                                                                                                            >  1:72585/1‑247 (MQ=255)                                                                                                                               cTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCCTTTGT‑GCCAGCGCTGCCGCCTGCTGTTCTGTACGACTTgggg                                                                                                           >  2:146448/1‑248 (MQ=255)                                                                                                                                    aaGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTGTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGGTGTCGCGGCATTGTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTGTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTAc                                                                                                                    <  1:43355/237‑1 (MQ=255)                                                                                                                                    aaGCAATTCCGGCGATAACGCTGCTAATGTCGCCAATTTTTTGCGAGCTGGTGGCGCTCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGACGCCTGCTGTTAGGTagcggctcgaagggcg                                                                                                      >  2:43355/1‑236 (MQ=255)                                                                                                                                      gCAATTCCGGCGATAACGCTGATAATGGCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGtt                                                                                                    <  1:259771/251‑1 (MQ=255)                                                                                                                                       cAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACCCCGCCCGTCTGCGCCGTTGGCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTa                                                                                                   <  2:209682/251‑1 (MQ=255)                                                                                                                                       cAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGtt                                                                                                    >  1:358321/1‑250 (MQ=255)                                                                                                                                        aaTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGGTTGT‑GCCAGCGATGCCGCCTGCTCTTCGGTACGACTTgagaagtcgtttt                                                                                                  >  2:307279/1‑242 (MQ=255)                                                                                                                                         aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGGTGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTAtt                                                                                                 >  1:586723/1‑251 (MQ=255)                                                                                                                                         aTTCCGTCGATAACGCTGATAATGTCGCCAACTGTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAATACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGAGTGATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGGTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTAtt                                                                                                 <  2:153690/251‑1 (MQ=255)                                                                                                                                           tCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGc                                                                                               >  1:56284/1‑251 (MQ=255)                                                                                                                                             cGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGAGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACCGGCCGCCGATTGT‑GCCAGGGATGCCGCCTGCTGTTCGGGAAGACTTGGGAGATCGTTATTCCCt                                                                                             >  1:255160/1‑251 (MQ=255)                                                                                                                                              gTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTg                                                                                            <  1:69066/251‑1 (MQ=255)                                                                                                                                              gTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GACAGCGCTGCCGCCCGCTGTTCGGGACGACTTGAGAGATCGTTATTGACTg                                                                                            >  1:368207/1‑251 (MQ=255)                                                                                                                                                cGATAACGCTGATTATGGCGTTAATTTTTGGCGGGCGGGTGGCGAGCTCCTGCATGGTGTGAGTCATGGTACTTACCTGAACACCTCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGGTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGTCGGTTAATTGCTCCATACGGGCGGCCGTTTGT‑GCTATCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTTCCTGCg                                                                                          <  2:353917/251‑1 (MQ=255)                                                                                                                                                 gATAACGCAGATAAGGTCGCCAATTTTTTGCGAGCTGGTGGCGATCACCTGCATGGTGTGAGTCATGGTTCTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTGTTGCCAGTTCCGACGACTGTCGTGCGTTATCGGCGTTTTGCCCTCCCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCTTTATTGCCTGcc                                                                                          <  2:329905/250‑2 (MQ=255)                                                                                                                                                      cGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTAAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAAtc                                                                                    <  1:457627/251‑1 (MQ=255)                                                                                                                                                        cTGATAATGTCGCCAATTTTTTACGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTAGCCCTACCGTGGCGGTTAATTGCACCATAATGGCTGACGTTTGTCGCGA‑CGGTGCCGCCTGCATTTCGGGACGCCTTGAGAGCTAGTTATTGCCTGCGGCAAtctc                                                                                  >  1:409657/1‑251 (MQ=255)                                                                                                                                                              aTGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGATAGATCGTTATTGCCTGCGACAAGCTCcgccat                                                                            <  2:26690/251‑4 (MQ=255)                                                                                                                                                              aTGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAt                                                                            <  2:419202/251‑1 (MQ=255)                                                                                                                                                                   gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCa                                                                       >  1:481849/1‑250 (MQ=255)                                                                                                                                                                             tGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTACATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCAtctc                                                           >  1:352080/1‑251 (MQ=255)                                                                                                                                                                             tGCGAACTGGTGGCGATCTCCAGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGATCCATACTGGCGGCCGTTTGT‑GCCAGCGAGGCCCCCTGCTGTTCGGTACGACTTGAGGGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTTCAtctc                                                           >  2:267125/1‑251 (MQ=255)                                                                                                                                                                                gAGCTGGTGGCGATCTCATGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGGGGACCTTTTT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGAACTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑AGGGCATCTCCTg                                                        >  1:424731/1‑251 (MQ=255)                                                                                                                                                                                  gCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCTCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTTTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGc                                                      <  2:554591/251‑1 (MQ=255)                                                                                                                                                                                  gCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCCA‑TGTGCCTCTCCTGGc                                                      >  1:466056/1‑251 (MQ=255)                                                                                                                                                                                               tcCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCCTACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGGGAGGTCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGACTCTCCCGGCTTCCCTTACGcac                                         >  2:172457/1‑251 (MQ=255)                                                                                                                                                                                                ccTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTa                                              <  1:434771/245‑1 (MQ=255)                                                                                                                                                                                                ccTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACg                                            <  1:504993/247‑1 (MQ=255)                                                                                                                                                                                                 cTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGGGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACAt                                       >  1:458770/1‑251 (MQ=255)                                                                                                                                                                                                   gCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATTGTTATTGCCTGGGGCCATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGcacagct                                     >  2:376999/1‑248 (MQ=255)                                                                                                                                                                                                   gCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGcac                                         <  1:376999/247‑1 (MQ=255)                                                                                                                                                                                                     aTGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACt                                   <  1:387660/251‑1 (MQ=255)                                                                                                                                                                                                              gTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTTCTGTTCGGTACGAATTGCGAGATCGTTTTTGTCTGCGCCAATCTCCCCGATCACCA‑‑TGGGCATCTCCTGGCTTCCCTTTCGCACATCACTTACCGAccc                          >  2:152260/1‑251 (MQ=255)                                                                                                                                                                                                              gTCATGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTccc                          <  2:56284/251‑1 (MQ=255)                                                                                                                                                                                                                 aTGGTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGGCCCACg                       <  1:439626/251‑1 (MQ=255)                                                                                                                                                                                                                  tGGTACTGACCTGAACACCGCCCGTCTGCGCCGCTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTCCCGTGGCGGTTAATTGCTCCCTAACGGCCGCCCTTTGT‑GCCAGCGACGCCCCCTCCCGTTCGGTACCCCTTTCGAGCTCGTTATTGCCTGCGACCATCTCCCCGAT‑ACCAA‑TATGCCTCCCCTGGCTTCCCTTACCCACATCCCATACCGTCCCCCGc                      >  2:482504/1‑251 (MQ=255)                                                                                                                                                                                                                   ggTACTGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTCTCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTTT‑GCCAGCGATGCCGCCTGCTGTTCGGGACGACTTGAGAGGTCGGTATTGCCTGCGGCAATATCCGCGCT‑AACAACTG‑GCCTCTCCTGGCTTCCCTTACCCACATCACTTACCCTCCCACGcc                     >  2:92854/1‑250 (MQ=255)                                                                                                                                                                                                                        tGACCTGAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTATTTGCTCCATACTGGCGGCCGTTTGT‑GCCCGCGATGCCGCCTGATGTTTGGTACGACTTGTGTGGTCTTTTTTGCCTGCGGCAATCTCCGCGGT‑ACCAA‑TGTGCCTCTCCTGGCTTCCCTTACGCACATCACTTTCCGTCCCACCCAAAGcc                >  2:284222/1‑251 (MQ=255)                                                                                                                                                                                                                         gACCTGAACACCGCCCGGCCGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCTCATCACTTACCGTCCCACGCAAAGCct               <  2:352080/251‑1 (MQ=255)                                                                                                                                                                                                                              gAACACCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTTGCTTCCCTTACGCACATCCCTTTCCGTCCCACGCAAAGCCtgctgc          >  1:503244/1‑251 (MQ=255)                                                                                                                                                                                                                                 cacCGCCCGTCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACTCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACAGCACTTACCGTCCCACGCAAAGCCTGCTGCATg       <  1:490597/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccgcccgtCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGt     <  1:576248/251‑1 (MQ=255)                                                                                                                                                                                                                                      cccgtcTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCtt  <  2:481849/251‑1 (MQ=255)                                                                                                                                                                                                                                      cccgtcTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCtt  <  1:447286/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          TCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGC‑ACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGT‑GCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGATCGTTATTGCCTGCGACAATCTCCGCGAT‑ACCAA‑TGTGCATCTCCTGGCTTCCCTTACGCACATCACTTACCGTCCCACGCAAAGCCTGCTGCATGGTCTT  >  NC_000913/1969750‑1970228

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