Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,587,112	C→T	A461V (GCA→GTA)	narQ →	sensor histidine kinase NarQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,587,112	0	C	T	86.2%	243.8 / 29.0	88	A461V (GCA→GTA)	narQ	sensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base T (49/26);  total (55/33)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTCGCTCTGCGGA  >  NC_000913/2586817‑2587406                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                        gCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTg                                                                                                                                                                                                                                                                                                     <  1:650653/300‑1 (MQ=255)          gCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCa                                                                                                                                                                                                                                                                                             <  2:602668/299‑1 (MQ=255)                gAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑ATACTTTCCGCCTGACGCTGCAGCAGGTGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCA‑CCAAACTGACCCTCGAC‑TACCGTCTGCC‑AACCCTGGAACTGGACGCGCAAATGCa                                                                                                                                                                                                                                                                                       >  2:420644/1‑299 (MQ=255)                  aTTGCACGACTCGCTTGCTAAGGTACTTTCTGACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCACACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCCGGTACTg                                                                                                                                                                                                                                                                                                     <  2:127396/283‑1 (MQ=255)                   ttGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAA‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGt                                                                                                                                                                                                                                                                                    >  1:526884/1‑299 (MQ=255)                                        gTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAAGTGCATTTGTTGCAAATTATTcgc                                                                                                                                                                                                                                                              >  2:622174/1‑300 (MQ=255)                                                        gTATCCAGTTGACTTTACTGAAGCGTTCGATACCGTAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGATGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAAt                                                                                                                                                                                                                                               <  2:209115/299‑1 (MQ=255)                                                          aTCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGc                                                                                                                                                                                                                                             >  1:424626/1‑299 (MQ=255)                                                               gTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAACTGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGAACCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGaa                                                                                                                                                                                                                                       >  2:611096/1‑300 (MQ=255)                                                                      ttACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGTCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑acgcc                                                                                                                                                                                                                               <  2:403060/300‑1 (MQ=255)                                                                      ttACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑acgcc                                                                                                                                                                                                                               >  2:649459/1‑300 (MQ=255)                                                                      ttACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑acgcc                                                                                                                                                                                                                               >  1:312446/1‑300 (MQ=255)                                                                      ttACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCATGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑acgcc                                                                                                                                                                                                                               >  2:609790/1‑300 (MQ=255)                                                                       tacggAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑acgcca                                                                                                                                                                                                                              <  1:216071/296‑1 (MQ=255)                                                                        acttaaGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑acgccaa                                                                                                                                                                                                                             <  2:120346/296‑1 (MQ=255)                                                                                gTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCATGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCAAGc                                                                                                                                                                                                                      >  2:502231/1‑299 (MQ=255)                                                                                        ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGNGGGAAATGCTGGATACGTTACAAACTCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACTCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGA‑ACGCAAACGAAAGCGAAATTGcc                                                                                                                                                                                                             >  2:548194/1‑300 (MQ=255)                                                                                           gAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCATGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTc                                                                                                                                                                                                          >  2:720846/1‑300 (MQ=255)                                                                                                          gCACAAAGTATCATGGCCGATTTTTCCCCGATGCGTTGAATGATGGCTTATCGGCAGTTACGCGAGCTGTTGAGCTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTTAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGGCCAAACTGACCCTCGACCTGCCGTCTGCCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCCACGCCAACGCCAGCGAAATCGCCGTCAGTTGc                                                                                                                                                                                                    <  2:621636/299‑1 (MQ=255)                                                                                                          gCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCCTTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGCATGCGATGAAGC‑ACGCCACCGCCAGCGAACCCGCCATCAGTTGCGTCACcgc                                                                                                                                                                                            >  1:388753/1‑299 (MQ=255)                                                                                                          gCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTCCTGGATGCGCCAATGCAGGTGCCTTTTTTGCAAATTCTTCGCGAAGCGGTGCTGACTGCGCTGAAGC‑CCACCCACGCCAACGAACTCACCGTCCGTTTCGTccccac                                                                                                                                                                                            >  1:388744/1‑294 (MQ=255)                                                                                                          gCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTNCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACcg                                                                                                                                                                                             <  1:150815/298‑1 (MQ=255)                                                                                                          gCACAAAGTATCATGGCCGATTTTTCCA‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATATGTTACAAAATCAAACCAGCG‑CCAAACCGACCCTCGAC‑TGCAGTCTGCC‑AACCATGGTACTGGATGCGCAAATGCAGGt                                                                                                                                                                                                                                                                                    >  2:385543/1‑212 (MQ=255)                                                                                                          gCACAAAGTATCATGGCCGATTTTTCCA‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTACGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAATTGCGTCACcgc                                                                                                                                                                                            >  2:264061/1‑299 (MQ=255)                                                                                                              aaaGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCgg                                                                                                                                                                                       <  1:515136/300‑1 (MQ=255)                                                                                                               aaGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGTTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGa                                                                                                                                                                                      >  1:116449/1‑300 (MQ=255)                                                                                                               aaGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGa                                                                                                                                                                                      >  1:116403/1‑300 (MQ=255)                                                                                                                                        aGGCGTTGAATGATG‑CTTATCTGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGTCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGTTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATc                                                                                                                                                                <  1:4167/299‑1 (MQ=255)                                                                                                                                        aGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGTTTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAGATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATc                                                                                                                                                                <  2:20740/299‑1 (MQ=255)                                                                                                                                        aGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGNGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATc                                                                                                                                                                <  2:587608/299‑1 (MQ=255)                                                                                                                                        aGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATc                                                                                                                                                                <  1:228249/299‑1 (MQ=255)                                                                                                                                        aGGCGTTGAATGATG‑CNTATCGGCAGTTACGCGAGNTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCANTCACACGGTTTATATc                                                                                                                                                                <  1:625026/299‑1 (MQ=255)                                                                                                                                         ggCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATcc                                                                                                                                                               >  1:747505/1‑299 (MQ=255)                                                                                                                                          gCGTTGAATGATG‑CTTNTCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTNACCGCGCCGGACGGCAATNACACGGTTTATATCCGt                                                                                                                                                             >  1:624223/1‑300 (MQ=255)                                                                                                                                          gCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGt                                                                                                                                                             >  1:564748/1‑300 (MQ=255)                                                                                                                                          gCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGt                                                                                                                                                             >  1:52372/1‑300 (MQ=255)                                                                                                                                            gTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAACTCGACGTCAGTTGCGTCACCGCGCAGGACGGCAATCACACGGTTTATATCCGTGa                                                                                                                                                           >  2:745774/1‑300 (MQ=255)                                                                                                                                                 atgatg‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTc                                                                                                                                                                                                 <  1:610641/257‑1 (MQ=255)                                                                                                                                                 atgatg‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCATGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTc                                                                                                                                                                                                 >  2:610641/1‑257 (MQ=255)                                                                                                                                                        cTTATCGGCAGTTACGNGAGCTGTTGA‑CTCCTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACNGccaccaacagggtta                                                                                                                                                                    >  1:629475/1‑266 (MQ=255)                                                                                                                                                        cTTATCGGCAGTTACGNGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGGGCCGGACNGCAACCCCACGGttt                                                                                                                                                                     >  1:629161/1‑279 (MQ=255)                                                                                                                                                        cTTATCGGCAGTTACGNGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACNGCAATCACAAGGTTtatcnccatgaccaaagcatt                                                                                                                                                 >  1:629468/1‑281 (MQ=255)                                                                                                                                                        cTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGAtt                                                                                                                                                 >  1:219386/1‑299 (MQ=255)                                                                                                                                                            tCGGCAGTTACGCGAACTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCTTGATAACGGGATTGGTAt                                                                                                                                            <  2:57483/300‑1 (MQ=255)                                                                                                                                                               gCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGNG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACNCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATc                                                                                                                                           >  2:574253/1‑298 (MQ=255)                                                                                                                                                               gCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGAGGTGCTGAATGCGATGAAGC‑ACGCCAACGCAANCGAAATCGCCGTCAGTTGCGGCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATc                                                                                                                                           >  1:140180/1‑298 (MQ=255)                                                                                                                                                                           gCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                                                                                                             >  1:723208/1‑300 (MQ=255)                                                                                                                                                                           gCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGNCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                                                                                                             >  1:152059/1‑300 (MQ=255)                                                                                                                                                                            cTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                                                                                                            >  2:462323/1‑300 (MQ=255)                                                                                                                                                                              gTTGA‑ATACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCACTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTAACCCTCGAC‑TGCCGTCTATC‑AACCCCGGTACTGGATGCGCAAACGCAGGTGTATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCAAACGCCAGCGAAACCGCCGTCATTTGCGTCACCgcgc                                                                                                                                                                                          >  2:108622/1‑236 (MQ=255)                                                                                                                                                                                tGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGTG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTAAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGGTTGGTATCGGTGAACCGAAAGaa                                                                                                                            >  2:168897/1‑296 (MQ=255)                                                                                                                                                                                 tggctactGTCCGCCTGACGCGGCAGCGGGCGGAGCTCCCCTCCGCAGTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACGGACCCTCGAT‑TGCCGTCTGTC‑AACCCTGGTACTGGATGCGCGAATGCGGATGCAGGTGTTGCAAATTATTCGCGAAGCGGTGCTGAAGGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACggtggtcgcgcggtttgtatc                                                                                                                                                                <  2:388753/257‑20 (MQ=255)                                                                                                                                                                                     tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGAATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAgg                                                                                                                    <  1:605936/300‑1 (MQ=255)                                                                                                                                                                                       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTTCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCAGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTc                                                                                                                  >  1:270559/1‑300 (MQ=255)                                                                                                                                                                                       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCAGAAGGTc                                                                                                                  >  1:364261/1‑300 (MQ=255)                                                                                                                                                                                       ctTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCAACGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTc                                                                                                                  >  1:364197/1‑300 (MQ=255)                                                                                                                                                                                             gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTAt                                                                                                             <  1:682805/299‑1 (MQ=255)                                                                                                                                                                                                 gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCt                                                                                                        >  1:467037/1‑300 (MQ=255)                                                                                                                                                                                                  aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTTATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCt                                                                                                        >  1:73608/1‑299 (MQ=255)                                                                                                                                                                                                  aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGAAGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCt                                                                                                        >  2:492350/1‑299 (MQ=255)                                                                                                                                                                                                  aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGNTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCt                                                                                                        >  2:549223/1‑299 (MQ=255)                                                                                                                                                                                                  aCGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCATGGTACTGGATGCGCAAATGCAGGTGCATTTGNTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGtat                                                                                                        >  2:549213/1‑297 (MQ=255)                                                                                                                                                                                                      tGCAGCAGGCGGATCTCCCCTCCGCATTTAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAat                                                                                                    <  2:747505/299‑1 (MQ=255)                                                                                                                                                                                                                         cTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCAGTCAGTTGCGTCACCGCGCCGGACGGCAATCACAAGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGGGCGGa                                                                                >  2:411327/1‑300 (MQ=255)                                                                                                                                                                                                                          tCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGNTGCATTTGTTGCAAATTATTCGCGAAGCGGTGNTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATc                                                                                                                                                                >  2:569144/1‑219 (MQ=255)                                                                                                                                                                                                                          tCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATc                                                                                                                                                                <  1:569144/219‑1 (MQ=255)                                                                                                                                                                                                                             gCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCNGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAACCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAc                                                                              <  1:552906/298‑1 (MQ=255)                                                                                                                                                                                                                                   aGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGt                                                                      <  1:437439/300‑1 (MQ=255)                                                                                                                                                                                                                                     ggAAGTGCTGGATACGTTACAAAAGCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTgg                                                                    <  2:153231/300‑1 (MQ=255)                                                                                                                                                                                                                                     ggAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTgg                                                                    <  1:574253/300‑1 (MQ=255)                                                                                                                                                                                                                                              ggATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGAc                                                           >  1:677613/1‑300 (MQ=255)                                                                                                                                                                                                                                              ggATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCAGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCAGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGAc                                                           >  2:471683/1‑300 (MQ=255)                                                                                                                                                                                                                                                         aaaaTCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCTTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAAc                                                <  1:609435/300‑1 (MQ=255)                                                                                                                                                                                                                                                           aaTCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCt                                              >  1:699973/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 aCCAGCG‑CCAAACTTACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCAAACGCGCGGAACGGCTAGGTTGGACGCTGACTTTTTCACAACCTTCcggc                                        <  2:626735/300‑1 (MQ=255)                                                                                                                                                                                                                                                                     gcg‑cCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCANTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAAACTTCcggcggc                                     >  2:547907/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                     tCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCTAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGt                      <  2:649780/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gccttCTGCC‑AACCCTGGCACTGTATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCTAAGCGGTGCTGAATTCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTTGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCt               <  2:186306/296‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTTGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCt               <  1:49456/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCt               <  1:303506/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gCCGTCTGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCCGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCt               <  2:715414/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGCATTAGCt               >  2:417754/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                  tGCC‑AACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTCGct         >  2:286404/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                  tGCC‑AACCATGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTc            >  2:285240/1‑297 (MQ=255)                                                                                                                                                                                                                                                                                                         cccTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCANTTGCGTCACCGCGCCGGACGGCAATCACACGGNTTATATCCGTGATAACGGGATTGGTATCGGTGACCCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTCGCTCTGc     >  2:565712/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                         cccTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCANTTGCGTCACCGCGCCGGACGGCAATCACACGGNTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTCGCTCTGc     >  2:565692/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                         cccTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCAAGTTAGTGAGTATTAGCTTTCGCTCTGc     >  2:148220/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                           cTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGANCGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTCGCTCTGCGGa  >  1:148121/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                        GCGACCATCGCCCGCGAATTGCACGACTCGCTGGCTCAGGTACTTTCTTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCC‑‑AGGCGTTGAATGATG‑CTTATCGGCAGTTACGCGAGCTGTTGA‑CTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCG‑CCAAACTGACCCTCGAC‑TGCCGTCTGCC‑AACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGC‑ACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTTCACAACCTTCCGGCGGCGGCACGTTAGTGAGTATTAGCTTTCGCTCTGCGGA  >  NC_000913/2586817‑2587406

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑