Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	989,901	G→T	16.9%	S152R (AGC→AGA)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,901	0	G	T	16.9%	172.3 / 7.2	94	S152R (AGC→AGA)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/46);  new base T (12/5);  total (41/53)
Fisher's exact test for biased strand distribution p-value = 2.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGC  >  NC_000913/989666‑990138                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ggCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAC‑GCCGGGAAG‑TCGGC‑GAATTTGCtt                                                                                                                                                                                                                                                   >  2:173858/1‑239 (MQ=255)       tttGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑Gc                                                                                                                                                                                                                                            <  2:1354006/239‑1 (MQ=255)                gtgCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGCGGGTTAAC‑‑CCGCGGAA‑GCGACATAATTTG‑TTTCCA‑GCGTGTCGaag                                                                                                                                                                                                                                   >  1:1311075/1‑237 (MQ=255)                  gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGc                                                                                                                                                                                                                                 <  1:659136/239‑1 (MQ=255)                  gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGc                                                                                                                                                                                                                                 <  1:738982/239‑1 (MQ=255)                  gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTACTGTTGCAGACGCTTAACGATGGTTTCTTGTGCTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGc                                                                                                                                                                                                                                 <  2:1076533/239‑1 (MQ=255)                   ccTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAg                                                                                                                                                                                                                                    <  1:61450/235‑1 (MQ=255)                   ccTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAA‑GCCGCGAAG‑TGGAC‑TAAATTTTTTTTCA‑GCGGGTCGAggggt                                                                                                                                                                                                                                >  1:587505/1‑237 (MQ=255)                        aaCCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAA‑GACGAGAAG‑GCGAC‑TCAATTGCTTTTCA‑GCGTGTTGGGGGCCTgcg                                                                                                                                                                                                                            >  1:731928/1‑238 (MQ=255)                             ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTTAC‑GCCGTAAAG‑GAGGC‑GAATTTGCTTTCCT‑GCCGGTGTAGGGCGTgccggcag                                                                                                                                                                                                                       >  2:998206/1‑232 (MQ=255)                             ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAA‑GCCCAGAAG‑TCGAC‑TAATTTTCTTTTCA‑CCGTGTCGAGGGGTTGCG‑CCAC‑g                                                                                                                                                                                                                     >  2:1175243/1‑238 (MQ=255)                             ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTTAA‑GCCGCGAAG‑GCGGC‑TAATATGCTTTTCC‑CCCTGTCGAGGGACTGCG‑CGActg                                                                                                                                                                                                                     >  2:61450/1‑237 (MQ=255)                                     ccTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAA‑ACCGCGAAG‑ACGGC‑ACAATTTCTTTTCA‑GCGGGGTGTGGGCGTGCG‑GAAC‑GCCgggccgc                                                                                                                                                                                                            >  2:557080/1‑234 (MQ=255)                                      cTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCg                                                                                                                                                                                                           <  2:150905/239‑1 (MQ=255)                                      cTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAAGCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCg                                                                                                                                                                                                           <  2:374055/239‑1 (MQ=255)                                      cTGTGTTCCCATCAGCGTGAGGGAAAGCCGACGCGCCAGTTAGTAGTGGCTGGGTCCTTCAAACCAGGATTCCTGTTGCAGACGGTTAACGATGGTTTCTTGTACTGCCCGAGACAGACGGCT‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCTACATATCAAGACCGGC‑GGTTAAC‑TCTGAGAAG‑TCGAC‑TAATTTTTTTTCCT‑GCGTGTCGAGGGCCTGCG‑CATA‑GTCGGCTTGCg                                                                                                                                                                                                           <  2:1037453/239‑1 (MQ=255)                                        gtgtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAC‑GCCGAGAAG‑ACGAC‑TAATTTGCTTTTCA‑GCGTGTCGGGGGGTTGCG‑CAAC‑GTCCGCCTGCGGt                                                                                                                                                                                                         >  1:139925/1‑239 (MQ=255)                                         tgtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAC‑GCCGAGAAGATGGGC‑AAATTT‑CTTTTCA‑GCGTGGCGAGGGCCTGCG‑CAAC‑CCCGGCCCGCgggg                                                                                                                                                                                                        >  1:1009214/1‑237 (MQ=255)                                          gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTg                                                                                                                                                                                                        <  2:366970/238‑1 (MQ=255)                                              cctccGGCGTGAGGGAAAGCCGCCGCGCCAGATCGTAGTTTCTGGTTCCCACAAACCAGGATTCCTGTAGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGagcga                                                                                                                                                                                                   <  2:392466/235‑1 (MQ=255)                                              ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGagcga                                                                                                                                                                                                   <  2:848897/239‑1 (MQ=255)                                              ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGagcga                                                                                                                                                                                                   <  1:106280/239‑1 (MQ=255)                                              ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGagcga                                                                                                                                                                                                   <  2:1349062/239‑1 (MQ=255)                                                aTCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTTAA‑GCCGAGAAG‑TCGGC‑TAATTTGCTTTTCC‑GCGTGTCGGGGGGTTTGG‑CCAC‑GTCGGCCCGCGGGGaggggta                                                                                                                                                                                                 >  1:386035/1‑234 (MQ=255)                                                aTCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGGGAAG‑ACGAC‑CAATTTGCTTTTCA‑GCGCGTCGAGGGGTTTCG‑CAAC‑GTCGGCCGCCGGTGggcgat‑a                                                                                                                                                                                                >  1:884048/1‑239 (MQ=255)                                                    gCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTAAC‑GCCCCGGAG‑GCGAC‑TAATTTGGTTTCCACGCG‑GTCGAGGGCTTGCG‑CAAC‑GTCGGGCTGCgggggggtgaagggg                                                                                                                                                                                             >  1:214580/1‑226 (MQ=255)                                                     cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGCTGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGCGAAG‑TCGGA‑CAATTTTTTTTTCA‑GCGTGTGCGGGGCTTGCG‑CAAC‑ATCCGCCTGCGGTGggcggt‑a‑tcgat                                                                                                                                                                                          >  1:531185/1‑239 (MQ=255)                                                     cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGagcgat‑a‑gcga                                                                                                                                                                                           <  1:1268279/238‑1 (MQ=255)                                                      gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATg                                                                                                                                                                                         <  2:1184633/239‑1 (MQ=255)                                                      gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATg                                                                                                                                                                                         <  1:660447/239‑1 (MQ=255)                                                       tGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCCCGAAG‑TTGAC‑CAATTTTCTTTTCA‑GCGTGTCGCGGGCTTGCG‑CAAC‑CCCGGCGCGCGGTGggggat‑a‑gcgtagt                                                                                                                                                                                        >  1:261511/1‑235 (MQ=255)                                                            gAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCa                                                                                                                                                                                    <  1:988885/238‑1 (MQ=255)                                                                                  aGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGCGAAG‑TCGAC‑TAATTTGCTTTTCA‑GCGGGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGGGAGCGAT‑A‑GCGATGTACCATTTCTCTGCTAAACGCCCGCaaa                                                                                                                                                             >  2:1103907/1‑239 (MQ=255)                                                                                  aGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTgcg                                                                                                                                                                                                                            <  1:1231273/180‑1 (MQ=255)                                                                                         ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTTAC‑CACGAGAAG‑GACAA‑TAATTTGCTTTCCA‑GCGTGCCGCAGGCTTGCG‑CCAC‑GCCGCCCCGCGGTGTGCGGT‑A‑GCGGTGTCACATTTTACTGATAACCGCCAGCaacgacact                                                                                                                                                       >  2:7315/1‑231 (MQ=255)                                                                                            gCCCACAACCCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAAAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTcc                                                                                                                                                    <  1:499147/238‑1 (MQ=255)                                                                                            gCCAACAAACCAGGATCCCGGTGGAAGCCGCTTAACGATGTTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGAAATATCAAGCCCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTCCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTcc                                                                                                                                                    <  2:587505/238‑1 (MQ=255)                                                                                             cccACAAACCAGGTTTCCTGTTGAGGACGCTTAAGGATGTTTTCTTGTACTTCGCGAGAAAAACGGCG‑AGGGGGGCAAAAATCCA‑TCAGAGGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCGGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑AAAC‑GTCGGCCGGGGGTGAGCGAT‑A‑GCGATGTACCATTTCACGGATAACCGCCAGCAAAGGAACTTCCCa                                                                                                                                                  <  2:238203/239‑1 (MQ=255)                                                                                                            tCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTTCTTTCCA‑GCGTGTCGCGGGCTTGCG‑CAAC‑GTCGGCCCGCGGGGAGCGAT‑A‑GCGATGTACCCTTTTCCTGTCAACCGCCAGGAAAGGAAATTCCCCGAGGAGGACTTCAc                                                                                                                                    >  2:908378/1‑238 (MQ=255)                                                                                                            tCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAA‑GCCGCGAAG‑GCGAA‑CAATTTTCTTTTCA‑GGGTGTCGAGGGCTTTCG‑CCAC‑GTCGGGCTGCGGGGGGCGGT‑A‑GCGATGTACCCTTTCACTGATAAACGCCAGCAAAGGAAATTCCCCGAAGAAGACtttcc                                                                                                                                    >  2:1191028/1‑235 (MQ=255)                                                                                                                  tGCAGCCGCTTAACGATGGTTTCTTGTCCTTCGCGAAAAAAACGGCG‑AGGGGTGCCAAAATCCA‑TAAGAGGGAGGCGCACCATAACAAGCCGGGC‑GGTTACC‑GCGGAGAAG‑TCGAT‑TAATTTGCTTCCCG‑GCGTGTCGAGGGTTTGCG‑CAAC‑GTCGGCCTGGGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACg                                                                                                                              <  2:26679/238‑1 (MQ=255)                                                                                                                  tGCAGACGCTTAACGATGTTTTCTTGTACTGCGAAAAAAAAACGGCG‑CCGGGTGCCAAAACCGA‑TCAGATGGAAGCGGGACATATCAAGCCCGGC‑GGTTAAC‑GCCGAGACG‑TCGAC‑TATTTTGCTTCCGA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGGGCGAT‑A‑GCGATGTCCCATTTCACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACgg                                                                                                                             <  2:82009/239‑1 (MQ=255)                                                                                                                   gCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGGTTAC‑GCCGAGCAA‑GCGAC‑TAATTTGCTTTCCA‑GCGGGGCGGGGGGCTGCG‑GAAC‑GGCGGCCTGCGGTGAGCGATAA‑GCAGTG‑ACCCTTTTACTGATAACCGCCAGCAAAGGGACTTCCCCGCGGATggatttcccccaccg                                                                                                                             >  2:195995/1‑224 (MQ=255)                                                                                                                   gCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATACCAAGACCGGC‑GGTTAAA‑GCCCAGAAG‑TCGAC‑CAATTTGTTTTCCA‑GCGGGTTGGGGGGTTGCGACAAC‑‑TCGGCCTCCGGGGGGCGAT‑A‑GCAGTGTGCCATTTTCCTGGTAAACGCCAGCAAAGGAACCTCCCAGAGGACGGCTTCCcaccactg                                                                                                                             >  2:363273/1‑236 (MQ=255)                                                                                                                         gataacggttGGTTTTTTGATCTTCG‑GAAAAAAACCGTGCACGGGTGCAAAAATCCA‑TAAGATGAAGTCGCAAAATAGCAAGACGGGC‑‑GTTAACGGCCGAGAGG‑CCGAC‑TATTTGGCTTTCCA‑GCGTGTCAGGGGCTTGCG‑AATC‑GTCGGCCTGCGGTGGGCATT‑G‑GCAATGTCCCTTTCCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑Ctt                                                                                                                     <  2:1125025/230‑1 (MQ=255)                                                                                                                             aaCGATGGTTTTTTGTACTTGGAGAAAAAAACGGCG‑AGGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGCCCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTGGAGGGCTTGCG‑AAAC‑GTCGGCCTGGGGTGAGCGAT‑A‑GCGATGTCCCTTTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAAc                                                                                                                 <  1:557080/239‑1 (MQ=255)                                                                                                                                   ggTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGGGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAA‑GGCGGCCTGCGGTGAGCGGT‑A‑GCGATGTACCATTTTACTGATAACCGCCAGCAAAGGAACTTCCCCGAGGATGACTTCACGCCCCGGGCCG‑CTTAAACGAata                                                                                                            >  1:641121/1‑238 (MQ=255)                                                                                                                                   ggTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAA‑ACGAC‑TAAATTGCCTTTCA‑ACCGGTCGAGGGGTTGCG‑CACC‑GTCGGCCCGCGGGGAGGGCT‑A‑GAGGTGGACCCTTTTTCTGATAACCGGCCGAAAAGGGACCTCCCAGCGGGAGAATTTACGCCACGGGGCG‑CTTAAAAGAata                                                                                                            >  2:510409/1‑238 (MQ=255)                                                                                                                                       tCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTTCTTTTCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGGTGTACCATTTCACCTATAACCCCCAGCAAAGGAACTTTCCCGAGGGTGACCTCACGCCACGGGCCG‑CTTTAACCAAAATccc                                                                                                        >  2:455202/1‑237 (MQ=255)                                                                                                                                       tCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCCAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCCGGTCGAGGGCTTGCG‑CAAC‑CTCGGCCCGCGGTGGGCGCA‑A‑ACGATGTACCCTTTTCATGTTAACCCCCAGCCAAAGAACTTCCCCCAGGATTACTTCCCCCCACGGGCCC‑CTTAAACgccaacccc                                                                                                        >  2:625725/1‑230 (MQ=255)                                                                                                                                        cTTGTACTGCGCGAAAAAAACGGCG‑AGGGGTGCAAAAATCCA‑TCAGAGGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTCCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTGGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGc                                                                                                      <  1:1103907/239‑1 (MQ=255)                                                                                                                                               tgcgagagaAAAACGGCG‑A‑GGGGGCACAAACCCAGACAGGTGAACGCGGGAAATATCAAGCCGGGG‑GGTAACC‑GCCGAGAGG‑CCGTC‑TAATTTTCTTTCCA‑GGGTGTCGAGGGTTTGCG‑AACC‑GTCGGCCGGGGGTGGGCGAT‑G‑GCGATGTACCTTTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCAt                                                                                               <  2:733590/234‑1 (MQ=255)                                                                                                                                               tccgggaaaaaaaCGGCG‑CCGGGTGCCAAAATCTA‑TCAGATGGAAGCGCGACATATCAAGACGGGC‑GGTAAGC‑GCCGAAATG‑TCGAC‑AAATTCTCTTGCCA‑GCGTGTCGAGGGCTTCCG‑CAAC‑GGCGGCCGGCGGTGAGCGAT‑A‑GCGATGTCCCATTCCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCAt                                                                                               <  2:47946/231‑1 (MQ=255)                                                                                                                                               tccgcgAAAAAACCGGCG‑CGGGGGGCCAAAGCCCT‑ACCGAGGGAACCGCGAAATAACAAGCCCGGC‑GGGTACC‑CCCTAGATT‑CCGAC‑AAATTTCTTTTCCA‑GCGTGGCGAGGGTTGGCC‑CACC‑GTGCGCCTGCGGTGGGTGAT‑A‑GCGATGTCCTTTTCCACGGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCAt                                                                                               <  2:252479/237‑1 (MQ=255)                                                                                                                                                 cgcgGGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCCAGAAG‑GTGGC‑TAATTTGCTTTCCA‑GCGGGTGGAGGGCTTGCG‑CCAC‑GGCGGCCCTCGGGGGGCGGT‑ATGGCTTGGACC‑TTTTACTGATAAACGCCCGCCAAAGAACGTCACCAAGGGTGGCATTCCGCCACCGGCCGCCTTaaaaaaaaaacacgcttcccata                                                                                              >  2:754975/1‑219 (MQ=255)                                                                                                                                                 cgcgAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGGGGGGTTGCG‑CAAC‑GTCGGCCCGCGGGGAGCGCG‑A‑GCGATGGACCATTTTCCTGATAACCCCCCGCAAAGGGACTTTCCCGAGGGTGACTTCACGCCCCCGGCCG‑CTTAAAACAATAACCacccttccctct                                                                                             >  2:1294113/1‑228 (MQ=255)                                                                                                                                                 cgcgAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑GCGAC‑TAATTTGCTTTCCC‑GCGTGTCGGGGGCTTGCG‑CAAC‑GTCGGCCTGCGGGGAGCGAT‑A‑GCGATGTACCATTTTACTGATAAACGCCAGCAAAGGAACTTCCCCGAGGATGACTTCACGCCCCGGGCCG‑CTTAAACGAATATCCAACTTGCCCTTa                                                                                             >  1:1160199/1‑239 (MQ=255)                                                                                                                                                 cgcgAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAA‑GCCGAGAAG‑TCGAC‑TAATTTGGTTTTCA‑GCGTGTCCGGGGCGTGCG‑CCAC‑GTCGGGCTGCGGTTAGCGAG‑A‑GCGATGTACCATTTTAATGATAACACCCAGCCAAGGAACTTCCCAGAGGATGGCGTTACGCCACGGGCCG‑CTTAAACGAAAAACCAACCTGCCagtt                                                                                             >  2:291933/1‑236 (MQ=255)                                                                                                                                                 cgcgAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGCTAAC‑GCCGCGAAA‑TCGGC‑TAAATTTCTTTCCC‑ACGTGTCGAGGGTTTGCG‑CAAC‑GCCCGCCTTCGGTGAGCGAC‑A‑GCGCTGGACCATTTGCCTGATAACCGCCAGCATCAGAACCTCCCAGAGGGTGACTTCCTCCCCCAGGCCG‑CCTTCACGAATACTGACCCTGCCATTa                                                                                             >  2:1232069/1‑239 (MQ=255)                                                                                                                                                      aaaaacCGGCG‑ACGGGTGCCAAAATCCA‑TAAGATGGAGGCGCGAATTACAAAGGCGGGC‑GGTTAAC‑GCCGAAAAG‑TCGAC‑TAATTTGTTTTCCA‑GCGTGCGGAGGGCGTGCG‑AAAC‑GCCGGTGGGCGGTGAGCAAT‑A‑GCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGt                                                                                         <  2:641121/237‑1 (MQ=255)                                                                                                                                                       aaaacCGGCG‑GCGGGGGCAAAAACCCA‑GCAAGTGGAGGCGCGAAAAATCAAGCCGGGC‑GTTTACC‑GCCGAAAGG‑TCGAC‑TAATTTGCTTTCCA‑GGGTGTGGAGGCTTTGCG‑AAAC‑GTCGGCCGGGGGGGAGCGAT‑G‑GCGTGGTACCATTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTg                                                                                       <  2:71015/239‑1 (MQ=255)                                                                                                                                                        aaaaCGGCG‑ACGGGTGCCAAAACCCA‑TCAGAGGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTgg                                                                                      <  1:455202/239‑1 (MQ=255)                                                                                                                                                             ggcg‑gcGGGTGCCAAAATCCA‑TCAGAAGGAAGCGCGACATAACAAGACCGGG‑GGTTAAC‑GCGGAGAAG‑CCGAC‑TAATTTGCTTCCCA‑GCGGGGCGAGGGCTGGCG‑CACC‑GTCGGCCGGCGGTGAGCGAT‑A‑GCGATGTACCATTCCACGAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                                                                 <  2:1190240/239‑1 (MQ=255)                                                                                                                                                              gcg‑gcgGTTGCAAAAACCAC‑TAAGAGGGAGCCGCGACATAACAAGCCGGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTTCTTTCCA‑GCGTGCCGAGGTTTTGCG‑CACC‑GTGGGCCTGCGGTGAGCGTG‑A‑GCGATGTCCTATTTCATTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATAGCCAGCTTGCCATTATCGTTGGACACg                                                                                 <  2:1315994/238‑1 (MQ=255)                                                                                                                                                                 gagggggGCAAAAACCCA‑TCAGAGGGAAGGGAAACATACCAGGCCCGGG‑TGTTACC‑GCCGAAAAG‑CCGAT‑TAATTTGCTTTCCG‑GGGTGTCGAGGGTTTGCG‑AACG‑GTGGGCCTGCGGTGGGAGTT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGt                                                                               <  2:731928/235‑1 (MQ=255)                                                                                                                                                                         ccAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGCGTTGAGGGGTTGCG‑CCAC‑GTCGGCCCGCGGTGAGCGAT‑A‑ACGATGTACCATTTTACTGATAACAGCCAGCAAAGGCACTTCCCACAGGATGACTTCACCCCCCGGGCCG‑CTTAAACCAATATCCAGCCTGCCCTTTTCGTTGGACAACGGGAC‑Ttttt                                                                      >  2:767963/1‑238 (MQ=255)                                                                                                                                                                          aaaaaCCCA‑TGAGATGGAAGCGGGACATATCAAGACCGGC‑GTTTACC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTCCCA‑GCGTGTCGAGGGCTTGCG‑AAAC‑GTCGGCCTGGGGTGAGCGAT‑A‑GCGATGTACCTTTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTc                                                                     <  2:386035/237‑1 (MQ=255)                                                                                                                                                                               tCCA‑TCAGATGGAAGCGCGACAAAGCACGCCCGGT‑GTTTACC‑GCCGAGAAG‑TCGAC‑TATTTTGCTTTCCG‑GCGTGCCGAGGGCTTGCG‑CACC‑GTCGGCCTGCGGTGAGCTAT‑A‑GCGATGTCCCATTCCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGtt                                                               <  2:150269/239‑1 (MQ=255)                                                                                                                                                                                ccA‑TCAGATGGAAGCGCAACTTATCAGGCCCGGC‑GGTTACC‑GCCGAGAAG‑TCGCT‑TATTTTGCTTCCCG‑GGGTGTGGAGGGCTGGCG‑CAAC‑GCCGGCCGGCGGGGAGCGAT‑G‑GCGAGGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTa                                                              <  2:1160199/239‑1 (MQ=255)                                                                                                                                                                                 cA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCCAGAAG‑TCGAA‑TAATTTGCTTTTCA‑GCGTGGCGAGGGGTTGCG‑CCAC‑GTCGGCCCGCGGGGAGGGGT‑A‑GCGATGTACCATTTCCATGATAACCCCCCGCACAGGGAATTTCCCGAGGATGGCTTTACGCCACGGGCCC‑CCTAAACGAAAATCCAGCTTGCCCTTCTCGTTGGACACGGCGCC‑TTTTTCCGGGGTac                                                             >  2:1287093/1‑238 (MQ=255)                                                                                                                                                                                          ggAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGGTTGCG‑CAAC‑GGCGGCCTGCGGTGAGCGGT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGGTGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCCTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAACGGa                                                       >  1:1022130/1‑237 (MQ=255)                                                                                                                                                                                           aacgcggcaCATAGCAAGCCCGGC‑GTTAACC‑GCCGAGAAG‑GCGAC‑TAATTTGCTTCCCG‑GGGGGTCGGGGGCTTGCG‑CAAC‑GTCGGCCGGGGGTGAGCGAT‑A‑GCGATGTACCATTCCCCTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGaa                                                      <  2:884048/229‑1 (MQ=255)                                                                                                                                                                                            aaGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAA‑TCGAC‑TAATTTGCTTTTCC‑GCGGGTCGAGGGCCTGCG‑CAAC‑GGCGGCCTGCGGGGGGCGGG‑A‑ACGATGTACCCTTTTACCGATAACCGCCCGCCAAGGAAATTCCCCCAGGATGACTTCACGCCACGGGCCG‑CGTAAACGAAAATCCCCATTGCCCTTATCTGTGGACACAGGGGC‑ATTGTCCGGGGTTAAGCGGGAcac                                                   >  2:289260/1‑237 (MQ=255)                                                                                                                                                                                                  gaaagctcacggcgggc‑ggTTAGC‑GCGGAGAAG‑TCGCC‑TAATTTGCTTTCCG‑GCGTGTCGAGGGCTGGCG‑AAAC‑GCCGTCCTGCGTGGAGGAAT‑G‑GCGAAGAACCAGTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                             <  2:531185/225‑1 (MQ=255)                                                                                                                                                                                                  gACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCG‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGTTGAGCGAT‑A‑GCGATGTACCATTTCACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                             <  1:1175243/239‑1 (MQ=255)                                                                                                                                                                                                   gcacgatcagggcgggtggTTAAC‑GGCGCGAAC‑GCGAC‑TAATCTTCTTTCCA‑GCGTGGCGCGGGCTTCCG‑AAAC‑GCCGGGCTGCGGTGAGCGAG‑A‑GCTATGTACCTTCTCTCTGATACCCGCCGGCAAAGGAACTTCCCAGAGGATGGCTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                             <  2:306707/222‑1 (MQ=255)                                                                                                                                                                                                   caatatCAGGCCCGGG‑GGTTAAC‑GCCGAGAAG‑TCGAC‑AAATTTGTTTCCCA‑GCGTGTGGGGGGTTGGCG‑CAAC‑GTCGCCCGGGGGTGAGCGAG‑A‑GCGATGTACCATTTCACTGAAACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                             <  2:139925/236‑1 (MQ=255)                                                                                                                                                                                                   acaaaacAAGACCGGG‑GGTTAAC‑GCCGAGAAG‑TAGCT‑TAATTTGTTTCCCA‑GCGTGGCGAGGGCTTGCG‑AAAC‑GCGGGCTGGGGGTGGGCGAT‑A‑GCGATGTCCCATTCCACTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                             <  2:314892/232‑1 (MQ=255)                                                                                                                                                                                                     agaatagcgccggggggTTAAC‑GCGGAGATG‑TAAAC‑GATTTCGCTTTCCG‑GCGTGGGGGGTGCTTGCG‑CAGC‑GTCGGCCGGCGGTGGCGGAT‑A‑GTGATGTACAATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                             <  2:897168/222‑1 (MQ=255)                                                                                                                                                                                                       atCAAGACCAGC‑GCTTAAG‑CCCAAAAAC‑TCGTG‑CAATTTGTTCTCCT‑GCGTGACGAGGGCTGGCA‑CATC‑GACAGCTTGCGGTAAGTCAT‑A‑AGGTTGTACTAATCCCCGAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCAAGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGCGTTAAAGCAGAACTCGCGTAACCAGt                                        <  2:1220753/239‑1 (MQ=255)                                                                                                                                                                                                         aaaGCCCGGC‑GGTTAAC‑GCGGAGAAG‑TGGAC‑TAATTTGCTTTGCA‑GCGTGTGGAGGGTTTGCG‑CAAC‑GTGGCCCTGCGGTGAGGAAT‑A‑GCGTTGTACAATTTCATGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTa                                      <  1:980556/238‑1 (MQ=255)                                                                                                                                                                                                               cggc‑ggTTTAC‑GCCGAGAAG‑TCGAC‑TAATTTTCTTTTCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GGCGGCCTGCGGGGAGCGAT‑A‑GCGATGTACCATTTTACTGATAACCGCCCGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGGGAC‑TTGTTCCGGGTTAAAACGGGACTCGCGCAACCCGTTAAGATAg                                >  2:120841/1‑239 (MQ=255)                                                                                                                                                                                                                ggc‑gtTTACC‑GCCGAAAAG‑TCGAC‑TAATTTGCTTTCCG‑GCGTGTCGGGGGCTTGCG‑AAGC‑GGGGCCCTGCGGTGAGCGAT‑A‑GCGATTTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                               <  1:448960/239‑1 (MQ=255)                                                                                                                                                                                                                ggc‑ggTTAAC‑GCCGAGACG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCGGGGGGGAGCGAT‑A‑GCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                               <  1:1287093/239‑1 (MQ=255)                                                                                                                                                                                                                 gggggTTAAC‑GCGGAGAGG‑TCGAC‑AAATTTGCTTTCCA‑GCGTGTGGAGGGCTTGCT‑CAAC‑GTCGGCCGGCGGTGGGGTAT‑G‑GCGATGTACATTTTCACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGGTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                               <  1:414600/236‑1 (MQ=255)                                                                                                                                                                                                                 gc‑ggTTACC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGTTTGCG‑CAAC‑GCCGGCCTGGGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                               <  1:173858/238‑1 (MQ=255)                                                                                                                                                                                                                 gc‑ggTTAAC‑GCCGGGAAG‑TCGAC‑TAATTTTCTTTCCA‑GCGGGTCGAGGGCTTGCG‑CAAC‑CTCGGCCTGCGGTGAGCGAT‑A‑GCGATGGACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCCGAGGGTGACTTCACGCCCCGGGCCG‑CTTAAACGAATATTCCGCTTGCCATTATCGTTTGACACCGGGAC‑TTGTTCCGGGGTAAAACGGGAATCGCGGCACCAGTTAAGATAATc                              >  2:1256634/1‑239 (MQ=255)                                                                                                                                                                                                                 gc‑ggTTAAA‑GCCGCGAAG‑TCGCC‑TAAATTTCTTTCCA‑GCGGGTCGAGGGGTTGCG‑CAAC‑GCCGGCCCGCGGTGAGCGAT‑A‑GCGATGTACCATTTCAATGATAACCGCCAGCAAAGGAAATTCCCCCAGGACGCCTTCACGCCCCGGGCCG‑CTTAAACGAACATCCCACTTGCCATTATCCTTGGGCCCGGTGACTTTTTTCTGGGTAAAAGC‑GCACCCGCGTAACCCGTTTAGATa                                 >  2:1074776/1‑236 (MQ=255)                                                                                                                                                                                                                        tatggCCGAGAAG‑CCGAT‑TATTTGCCTTTCCA‑GCGGGTCG‑GGGGTTGGC‑CAACAGTCGGCTGGCGGGGAGGCAG‑A‑GCGATGTCCCATTTCACTTATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                         <  2:373708/235‑1 (MQ=255)                                                                                                                                                                                                                                agaAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                 >  2:335453/1‑239 (MQ=255)                                                                                                                                                                                                                                      tCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag            >  1:406275/1‑239 (MQ=255)                                                                                                                                                                                                                                              aTTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaa     >  1:305231/1‑239 (MQ=255)                                                                                                                                                                                                                                              aTTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaa     >  1:197312/1‑239 (MQ=255)                                                                                                                                                                                                                                                 tGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CACC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc  <  1:1276399/239‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     GGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCG‑ACGGGTGCCAAAATCCA‑TCAGATGGAAGCGCGACATATCAAGACCGGC‑GGTTAAC‑GCCGAGAAG‑TCGAC‑TAATTTGCTTTCCA‑GCGTGTCGAGGGCTTGCG‑CAAC‑GTCGGCCTGCGGTGAGCGAT‑A‑GCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCG‑CTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAC‑TTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGC  >  NC_000913/989666‑990138

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