Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	989,959	A→C	17.6%	V133G (GTA→GGA)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,959	0	A	C	17.6%	128.5 / 4.1	75	V133G (GTA→GGA)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base C (2/12);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 5.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAA  >  NC_000913/989723‑990196                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 gggAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTAc                                                                                                                                                                                                                                                 <  1:992390/238‑1 (MQ=255)      aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACCAATTTGC‑TTTCCAGCGTGTCGGAGGGCTGGGCC‑AAC‑TC‑GGCCTGCGGTGAGCGCTAGCGCTGTACCActt                                                                                                                                                                                                                                            >  2:897698/1‑238 (MQ=255)      aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTAGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGGTTGCGC‑CAAGTC‑GGCCTGCGGGGAGCGAGAGCGATGTAcccattt                                                                                                                                                                                                                                           >  1:148373/1‑234 (MQ=255)                     tCGTATTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGA‑CGTTTAACTATGGTTTCTTTTACTACGAGAAAAAAGCGGCGACGGGTGCCAAAACCCATCAG‑ATGGAAGCGCGACATATCAAGCCCGGGGGTTAACGCCGAGAGGGCGACTAATTTGC‑TTCCCAGCGTGTC‑GAGGGTTGGCGC‑AACGCC‑GGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑cc                                                                                                                                                                                                                            <  2:270343/237‑1 (MQ=255)                             gcggtggCCCACAAACCAGGATTCCGGTTGCAGA‑CGTTTAACGATGTTTTCTTTTACTTCGCGAAAAAAAGGGCGACGGGTCCCAAAATCCACCAG‑ATGGAAGCCCGACATACAAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGC‑TTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGaa                                                                                                                                                                                                                  <  2:148373/233‑1 (MQ=255)                              cTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCCGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGGGATGGACCATTTTACCG‑ATAACCG‑CCAGCAAAAGa                                                                                                                                                                                                                   >  1:1209956/1‑237 (MQ=255)                               tGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GGGGGCTTGCGC‑CACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTG‑ATAACCC‑CCCGCACAAGAACt                                                                                                                                                                                                                >  2:428772/1‑239 (MQ=255)                                ggTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTGC‑TTTCCAGCGTGTC‑GGGGGGGTGCGC‑AAACTC‑GGCCTGCGGGGGGCGATAGCGCTGTCGCAATTTACTG‑ATAAACGCCCACCAAGGGAACt                                                                                                                                                                                                                >  2:44251/1‑239 (MQ=255)                                   gCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGGACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGACAAATCGACTAATTTTC‑TTTCCCGCCTGGT‑GAGGGCTTAAGC‑AAAGTC‑CGCCTGCGGTGAGCGAAAAGGGTTGACCCTTTTCCTG‑AGACACG‑CCAGCAAAGGGAATTcgc                                                                                                                                                                                                            >  2:649585/1‑237 (MQ=255)                                   gCCCACAAACCAGGATTCCTGTGGCAGA‑CGCTTAACGATGTTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAG‑ATGGAAGCGCGACATATCAAGCCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTcc                                                                                                                                                                                                             <  1:1212552/238‑1 (MQ=255)                                        caAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGC‑TTTTCAACGGGGC‑GAGGGCCTGCGC‑AACGCC‑GGCCCGCGGTGAGCGATAACGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCCGAgg                                                                                                                                                                                                       >  1:136552/1‑239 (MQ=255)                                          aaCCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCAAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTACCGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGAt                                                                                                                                                                                                     <  1:428772/239‑1 (MQ=255)                                           aCCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAAGTGACTAATTTGC‑TTTTCAACCGGTC‑GAGGGGCTGCGCAAACGTG‑GGCC‑GCGGGGGGTGGAAGAGGTGTACCATTTAACTG‑ATAACCC‑CCAGGAAAGAAACATCCCAGAGGaag                                                                                                                                                                                                    >  2:934169/1‑237 (MQ=255)                                                gATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTTCAACGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGGCCGCGGGGGGCGATAGCGGTTTACCATTTCACCG‑ATAAACG‑CCCGCAAAGGAAATTTCCCGAGGGTGACtt                                                                                                                                                                                                >  2:519363/1‑238 (MQ=255)                                                 aTTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                                               <  1:484931/238‑1 (MQ=255)                                                  tttccgTTGCAGGCCGCTTAAGGATGGTTTCTTCTACGTCGCGAAAAAACCGGCGACGGGTGCAAAAAGCCATAAG‑ATGGAAGGCCGAAATATCAAGACCGGGGTTTAACGCCGAAAAGTCGATTAATTTGT‑TTCCCAGCGTGTC‑GAGGGCTTGCGC‑ACCGTC‑GGCCTGCGGTGGGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                                               <  2:784158/233‑1 (MQ=255)                                                     cTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTTC‑TTTTCAGCGGGTC‑GAGGGCTTGCGC‑AACGGC‑GGCCTGCGGTGAGCGGTAACGATGTAACATTTTACTT‑ATAAACG‑CCCGCAAAGGAACTTCCCCGAGGGTGACTTTACGc                                                                                                                                                                                           >  1:532900/1‑238 (MQ=255)                                                     cTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAGTCGACTCATTTTC‑TTTCCCGCGGGTC‑GGGGGGTTGCGC‑AACGCC‑CGCCCGCGGGGGGGGGTAGCGGTGTACCATTTTCATG‑ATAAACC‑CCCGCAAAAGGACATCCCAGGGGGAGACtttccgc                                                                                                                                                                                           >  2:182688/1‑233 (MQ=255)                                                                            ttttCTTGTACTGCGGAAGAAAAAGGGCGAGGGGTGCCAAAATCCA‑CAGAAGGGAAGCGCGACATATCAAGACGGGCGTTTAACGCGGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGA‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAatat                                                                                                                                                                     <  1:967773/237‑1 (MQ=255)                                                                                 ttgtgcggcgaaagaAAAAGGGGGAGGGGTCCCAAAACCAAGAAA‑GGGGAGGCGCGAAATATCAAGACGGGCGGTAAACACCGAGAGGCCGAC‑AATTTGCTTTTCGAGCGTGTC‑GAGGGCTTGCAA‑AACGGC‑GGCCGGGGGTGAGCGATAGCGATGAACTATTCCACAG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCt                                                                                                                                                               <  2:1185405/227‑1 (MQ=255)                                                                                  tGTCCTTGGGGAAAAAAGCGGGGACGGGGCCAAAAACCCATCAG‑ATGGAAGCCGGACATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTGGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGA‑ACCGTC‑GGCCTGCGGGGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCtt                                                                                                                                                              <  2:961553/239‑1 (MQ=255)                                                                                         cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTTAAGCCGAGAAGTCGACTAATTTGC‑TTTTCAGCGGGTC‑GAGGGGTTGCGC‑AACGGC‑GGCCCGCGGGGAGCGATAGCGGTGTACCATTTCACCG‑ATAACCC‑CCAGCAAAGGAACTTCCCAGAGGGTGACTTTACGCCCCGGGCCGCCTAAAcagaaatacagatttccctta                                                                                                                                                       >  2:1117046/1‑219 (MQ=255)                                                                                         cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGGGTC‑GGGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCAATG‑ATAACCG‑CCCGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGc                                                                                                                                                                                    >  1:1235500/1‑210 (MQ=255)                                                                                         cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCCGCGGTGGGCGATAGCGATGTAACATTTCACTG‑CTAAACG‑CCAGCAAAAGAACTTCCCCCAGGATGACTTCACGCCACCGGCCGCTTAAACGCATATCCAACTTTCCCttt                                                                                                                                                       >  1:1001106/1‑238 (MQ=255)                                                                                         cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGGACAACAAATTGC‑TTTACAGAGGGTCGGAGGGC‑TGCGC‑GAAATC‑GGGCTTCAGGGGGCGATAGACATTTCACATTTCCTCG‑CAAAACG‑CCCGCCAAGGGACTTACCAGAAGAAGGGTTCCAGCCACCGGGCGCCTAAACCGATATTCCACTTGCCACTa                                                                                                                                                       >  2:210005/1‑239 (MQ=255)                                                                                                 aaccGGCGAGGGGGGCAAAAACCCATCAG‑ATGGAAGCGCGACATACCAAGACCGGCGGTAAACGCCGAGAAGTCGACTAATTTGC‑TTCCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGa                                                                                                                                               <  1:296234/236‑1 (MQ=255)                                                                                                    cgcgggcggtccccAAAACCCATCAG‑AGGGAACCGCAACATAGCAAGCCCGGCTGTTAACGCAGAAAGGTCGATTATTTTGT‑TTCCCAGCGGGTC‑GAGGCCTTGCGC‑AACGTC‑GGCCTGGGGTGAGCGATAGCGATGAACCATTTCACTG‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGacac                                                                                                                                            <  2:672391/227‑1 (MQ=255)                                                                                                               gcaaaaacacATCAG‑ATGGAAGCGCGACATATCAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGC‑TTCCCGGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCGGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  2:532900/230‑1 (MQ=255)                                                                                                               gacAAAACCCATAAA‑GTGGAAGGCCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTCT‑TTTCCAGCGTGTC‑GGGGGTTGGCGC‑ACCGCC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  2:849192/237‑1 (MQ=255)                                                                                                                ccAAAATCCATCAG‑AGAAAAGCGGAAAATATAAGGACCGGCTGTTAACGCGAATAAGTCGACTATTTGGC‑TTCCCGGCGTGTG‑TGGGGCTTGCGC‑ACCGTG‑CGCCGGCGGTGGTCGAGAGCGAGTTCCCATTTCACTG‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  2:943213/238‑1 (MQ=255)                                                                                                                ccAAAACCCATCAG‑ATGGAAGCGCGACATATCAAGGCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑tt                                                                                                                                <  1:519363/238‑1 (MQ=255)                                                                                                                          tCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTG‑ATAACCG‑CCCGCCAAGGAACTTCCCAAGAGATGACTTCACGCCAACGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCATTGGACACAGGGGCTTG‑TTCCGGGGTaaa                                                                                                                      >  2:572585/1‑238 (MQ=255)                                                                                                                                 ggAAGCGCGAAAGATCAAGACCGGGGGTTAACGCCGAGAAGCCGACTAATTTGC‑TTCCCAGCGTGTC‑GAGGTCTTGCCC‑AACGTC‑GGCCGGGGGTGAGCAATAGCGATGTACCTTTCCCCTG‑AAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGaa                                                                                                                <  2:1209956/238‑1 (MQ=255)                                                                                                                                         aaCATATCAAGCCCGGGGGTTAACGCCGGGAGGTCGACTATTTTGT‑TTCCCAGCGTGTC‑GGGGGTTTGCGC‑AACGGC‑GGCCTGCGGGGGGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                                       <  2:289057/238‑1 (MQ=255)                                                                                                                                          aCATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                                       <  1:44251/238‑1 (MQ=255)                                                                                                                                             tatCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTTC‑TTTCCACCGTGTC‑GAGGGCTTGCGC‑CACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCCTTTCACTG‑ATAACCG‑CCCGCAAACGCACTTCCCCGAGGATGACTTCCCCCCACCGGCCGCTTAAACGAATATCCCGCCTGCCCTTATTCTTGTACACGGTGGCTTGTTTCCGGGTTAAA‑CGGACAACGCGGAACCac                                                                                                   >  2:689178/1‑238 (MQ=255)                                                                                                                                                cAAGCCCGGGGTTTAACGTCGAAACGTCGGCTAGTTGGC‑TTCCCAGGGTGTG‑GAGGGTTTGCCC‑AAGGTC‑GGCGCGGGGGGACGGATAGCGTGGTCCCATTTCTCGA‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTa                                                                                                <  2:1226803/239‑1 (MQ=255)                                                                                                                                                       gggggTAAACGCCGGGAAGGCGAATAATTTGC‑TTTCCAGCTGGTC‑GAGGGCATGCCA‑ACCGTC‑GGCTGGCGGTGACGAATAGCGATGTACAATTGCCCTG‑AGAACCG‑CCAGCAAAGGAACTGCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                                         <  2:157560/236‑1 (MQ=255)                                                                                                                                                        gcggTTTAAGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGGGTC‑GGGGGCTTGCGC‑AACGTC‑GGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCGTTGGACACGGGGGATTG‑TTCCGGGGTAAAGCGGAACTCCGGGAAACAGTTaaaaaaaac                                                                                        >  2:686493/1‑232 (MQ=255)                                                                                                                                                        gcggTTAACGCCGAGAAGTCGACTAATTTTC‑TTTCCAGCGTGTC‑GGGGGCTTGCGC‑AACGGC‑GGCCTGCGGGGAGCGAGAGCGGTGTACCATTTTACTG‑ATAAACG‑CCCGCCAAGGAACTTCCCCGAGGGTGACTTTCCCCCACCGGCCGCTTAAACGAAAATACACCTTGCCCTTATCGTTTGACACGGTGACCTG‑TTTCGGGTTCAAGCGGAACTCCCGTAAACCGTTACGATAGt                                                                                         >  2:631283/1‑238 (MQ=255)                                                                                                                                                         cggTTACCGCCGAGAAGTCGACTAATTTGC‑TTCCCAGGGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGAGATTTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                        <  1:740722/238‑1 (MQ=255)                                                                                                                                                         cggTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCGGCGTTGC‑GTGTGCTTGCGA‑ACCGGC‑TGCGTGCGGGGATCGATAGCGATGTCCCATTTCCCTA‑AAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                        <  2:136552/238‑1 (MQ=255)                                                                                                                                                                     agaAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                                           >  2:379034/1‑239 (MQ=255)                                                                                                                                                                     agaAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                                           >  2:116645/1‑239 (MQ=255)                                                                                                                                                                       aaGTCGACTAATTGGC‑TTTCCAGCGTGTC‑GGGGGCTTGCCA‑ACCGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                                         <  1:901034/239‑1 (MQ=255)                                                                                                                                                                        agtcatagaaTTTGT‑TTCCCACTGTGCC‑GAGGGTTTGCGC‑CCCGTG‑GGCCTGCGGGGAGAGTTAGTGATGTACCATTGCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                                         <  1:1158402/230‑1 (MQ=255)                                                                                                                                                                         gTCGGCTAATTTGC‑TTTCGTTGGTGGG‑GAGGTCTTGCGC‑ATGGCC‑GCCCTGCGTGGAGGGAGAGCGATGTCCTTTCTCATTG‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                                        <  2:379445/238‑1 (MQ=255)                                                                                                                                                                          tCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag                                                                      >  2:1218794/1‑239 (MQ=255)                                                                                                                                                                             aCTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggc                                                                   >  1:204556/1‑239 (MQ=255)                                                                                                                                                                                    ttttctttccaggcttcggggggTTTGCGC‑ACCGTC‑GCCCTGCGGTGCGCGATATCAATGGACCATTCCACTG‑AAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaaa                                                              <  1:1023389/219‑1 (MQ=255)                                                                                                                                                                                     tgcttccgggCGTGTG‑GGGGGCTTGCG‑‑AACGGCGGGCGAGCGGTGGGGGATGGCGGTGTTCCATTTCGCTA‑ATACCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                                            <  2:20286/230‑1 (MQ=255)                                                                                                                                                                                         ttCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                                           >  1:707889/1‑236 (MQ=255)                                                                                                                                                                                            cAGTGTGTC‑GGGGGTTTCGGC‑AATCGC‑GGCCTGCGTGGGCCAATAGCTATGTACCATTTCACTT‑ATATCCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                                     <  1:865714/239‑1 (MQ=255)                                                                                                                                                                                                gtgtC‑GAGGGTTTGCGC‑AACGCC‑GGCCGGCGGGGAGCGATAGCGATGTACAATTTCACAG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                                  <  1:416752/238‑1 (MQ=255)                                                                                                                                                                                                gtggC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                                  <  1:1218794/238‑1 (MQ=255)                                                                                                                                                                                                 tgtC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                                                >  1:819315/1‑239 (MQ=255)                                                                                                                                                                                                             tGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  1:731040/1‑239 (MQ=255)                                                                                                                                                                                                              gcgc‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  1:1091819/1‑238 (MQ=255)                                                                                                                                                                                                              gcgc‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                     >  2:249316/1‑238 (MQ=255)                                                                                                                                                                                                                           gCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgct                          >  2:334793/1‑238 (MQ=255)                                                                                                                                                                                                                             gggcgtgaggggatagcgatGT‑CCCTTTCACCGAAAAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:548900/228‑1 (MQ=255)                                                                                                                                                                                                                             cggCGGTGAGCGATAGCGATGTCCCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:1152536/237‑1 (MQ=255)                                                                                                                                                                                                                             cgcgggggAGCGGAAGCGATTTAACTTTCCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:600998/232‑1 (MQ=255)                                                                                                                                                                                                                             cTGGGGTGAGCGATACCGTTGTGCCATTTCACTA‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:680435/239‑1 (MQ=255)                                                                                                                                                                                                                             cTGCGGTGGGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa                       <  1:182688/239‑1 (MQ=255)                                                                                                                                                                                                                             cTGCGGTGAGGGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctgc                        <  1:1117046/238‑1 (MQ=255)                                                                                                                                                                                                                                cGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc                     >  1:467246/1‑238 (MQ=255)                                                                                                                                                                                                                                cGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc                     >  2:678394/1‑238 (MQ=255)                                                                                                                                                                                                                                     agcgatagcgatGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaa                                                                               >  1:490444/1‑175 (MQ=255)                                                                                                                                                                                                                                          tagcgatGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGaa              >  1:988364/1‑235 (MQ=255)                                                                                                                                                                                                                                             cgatGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc       >  2:1020312/1‑239 (MQ=255)                                                                                                                                                                                                                                              gttGTACCATTTCCCTG‑ACAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc       <  2:386943/236‑1 (MQ=255)                                                                                                                                                                                                                                               atGTACCTTTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                                   <  1:747319/161‑1 (MQ=255)                                                                                                                                                                                                                                               atGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                                                                                   >  2:747319/1‑161 (MQ=255)                                                                                                                                                                                                                                                  tACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGaa  >  1:290137/1‑239 (MQ=255)                                                                                                                                                                                                                                                   aCCATTTCACTG‑AAAACCG‑CCATCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGaa  <  2:669835/238‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 GGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGA‑CGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGC‑TTTCCAGCGTGTC‑GAGGGCTTGCGC‑AACGTC‑GGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTG‑ATAACCG‑CCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTG‑TTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAA  >  NC_000913/989723‑990196

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