Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	783,293	A→G	16.1%	E43G (GAG→GGG)	pnuC →	nicotinamide riboside transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	783,293	0	A	G	16.1%	179.6 / 10.3	81	E43G (GAG→GGG)	pnuC	nicotinamide riboside transporter
Reads supporting (aligned to +/- strand):  ref base A (29/38);  new base G (12/2);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 6.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGATTGGTCT  >  NC_000913/783055‑783530                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      atgaAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATTGATCTAAGCGGTAAGGACG‑CTCGCCGGGTTGCTGTG‑TATTGGCCCTGCCGGTC‑TTga                                                                                                                                                                                                                                                       >  2:1020550/1‑239 (MQ=255)           gCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACACATAGGGGCAGGCGGTTATTGACTCTCATGGGTCGGAGCGGTAGGCCCG‑CTCGCCGGGGTGGTGGG‑TATTTGCTTTTCCAGTC‑TGGGGA‑AAAtaaac                                                                                                                                                                                                                                            >  1:1090926/1‑235 (MQ=255)                 gagGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAacta                                                                                                                                                                                                                                       <  2:953008/238‑1 (MQ=255)                      ggctgggtccgtttaatgttttGCTTGTTTTTGCGGCTGCACTACGTGGATAACGAATAATAAGGCGTAACGTTACGGTGTGGGGGAAAGATGGATGTTTTTAGTGTGCAGAATAGCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAc‑ttct                                                                                                                                                                                                                                 <  1:16816/217‑1 (MQ=255)                                     aTCGGAAGCTGGATTTTTCGGCTTCACTACGTGGATAACGGATAATAAGGCGGAACGTTACGGTTTGGGGGAAAGATGGATTTTTTTAGTTGGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGt                                                                                                                                                                                                                <  2:570379/239‑1 (MQ=255)                                                gATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTTTGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCCGGGTGCTGTG‑TATTGGCCTTGCCCGGC‑TTGAGA‑AGATCAGCAACTAC‑TTCTTTTGCCTGTT‑TCACGTCACCTTGtttt                                                                                                                                                                                                     >  2:781666/1‑238 (MQ=255)                                                   tttGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCa                                                                                                                                                                                                  <  2:1310588/239‑1 (MQ=255)                                                             cacTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATtttc‑tt                                                                                                                                                                                       <  2:191370/239‑1 (MQ=255)                                                                   gTGGATAACGAATAATAAGGCATAACGTTACGTTTTGGGGGAAAGATGGATTTTTTTAGTGTGAAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGa                                                                                                                                                                                  <  2:1030384/238‑1 (MQ=255)                                                                          aCGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTTTGATCTCTCATGGGTCGAAGCGGGAGGCACG‑GTCGCCGGGGTGCTGTG‑GTTTTGCCTTGCCAGCC‑TGGGGAGAGA‑CAGAAAATAC‑TTTTTTTGCCTGAT‑CAAAGTACCCTTGTTTGGCATTTTTTTT‑TTTTAGATTC‑a                                                                                                                                                                             >  2:909377/1‑235 (MQ=255)                                                                              ataataAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATg                                                                                                                                                                     <  2:195349/239‑1 (MQ=255)                                                                                aataaGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGTCGGTTATGATATCTCCAGGATCGAAGCGGTAAGCACC‑ATCCGCGGGTTTCTGGG‑TGTTTGCCTTGCCAGCT‑TCGGGA‑AAAACCGCCACTACTTTCTTGGCCCGGTT‑AAACTTCCCCTTG‑TTGGACTTAATTTC‑TTTTAAATTC‑aagcggtagga                                                                                                                                                                   >  1:540980/1‑229 (MQ=255)                                                                                    aGGCGTAACGTTACGATTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATTCCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGcc                                                                                                                                                               <  2:231258/239‑1 (MQ=255)                                                                                    aGGCGTAACGTTACGATTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGcc                                                                                                                                                               <  2:804707/239‑1 (MQ=255)                                                                                       ggTAACGTTTCTGTTTGGGGTAAAGTTGGATTTTTTTAGTGGGCAGAATTTCCGGGTACATATACCAATAGGGGCAGGCGGTTTTGATCTCTCATGGAACGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTg                                                                                                                                                             <  1:788774/237‑1 (MQ=255)                                                                                       ggTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCt                                                                                                                                                              <  2:1247042/236‑1 (MQ=255)                                                                                         tAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGGGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGTCATTATTTTC‑TTTCAGATTC‑CACTGTATGCCAGCCTCC‑Ta                                                                                                                                                         >  2:685725/1‑239 (MQ=255)                                                                                               tACGCATTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATTGATCGAAGCGGTAGGCGCG‑ATCGGCGGGGTTCTTTG‑TATTTGGCTTGCACGTC‑TTGGGA‑CGATGCGCCACAAC‑TTCTTTTGCCTGATCCAACG‑CACCCTGGTTGGGCTTTTTTTCTTTTCAGTTTC‑ACGTGTTTGCC‑GCCTgtctttaaag                                                                                                                                                    >  1:376699/1‑230 (MQ=255)                                                                                                        ggggAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGTTCTTTTCTTGATCGAAACGGGAGGGACG‑ATCGCCGGGGTGC‑GTGTTTTTTGCCTTGGCCGTC‑TGGGGA‑AGATCACGGACTACTTTCTTTGG‑CTGGT‑CATCCTGAACCTGGTTTGCCTTATTTTT‑TTTCCATTTC‑AGCCGTTTTCCAGCCCGC‑TTTTTAAG‑AGG‑tttttt                                                                                                                                        >  1:1098783/1‑238 (MQ=255)                                                                                                             aaGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGGTGCTGTG‑TATTGGCCTTGCCAGGC‑TGGGGA‑AGATCAGCAAATAC‑TTCTTTGGCCTGAT‑CAACGTCCCCTTGTTTTGCATTTTTTTT‑TTTTCGATTT‑ACCTGTATGCCCGCCTGC‑TTTTTCAA‑GGG‑TTTTTTTTTg                                                                                                                                    >  2:852917/1‑238 (MQ=255)                                                                                                                 tGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCcg                                                                                                                                 <  2:267113/237‑1 (MQ=255)                                                                                                                 tGGATTTTTGTAGTGTGCAGAATAGCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCcgcg                                                                                                                               <  2:651047/239‑1 (MQ=255)                                                                                                                  ggATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATGGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGa                                                                                                                              <  2:540980/239‑1 (MQ=255)                                                                                                                   gATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTTTGATATCTCATTGGTCCAAGCGGTCAGGACG‑ATCGCCCGGGTGGTTGG‑TATTGGGGTTTTCAGGT‑TGGGGA‑AGATCAGCAAACAC‑CTTTTTGGCCTGGT‑TCAAGTGCCCCTGTTTGGGCTTATTTTT‑TTTTATATTC‑ACCTTTAATCCCACCTGT‑TTTTACAGAGTGTTTTTTTTTTGCcgca                                                                                                                               >  1:122175/1‑238 (MQ=255)                                                                                                                   gATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGAAGGCGGTTATAATCCCTCAATGATCGCAACGGTAAGAAAGCTTCGCC‑GGTCGTTGTG‑TATTTGCCTTCCCAGCCTTGGAGA‑GCAT‑AGCAACCAC‑TTCTTTTGCCTGAT‑CCACGCCCACTTGTTTTGTCTTATTTTT‑CTTTAGACTT‑ACCTGTATGCCAGCCTGC‑TTTTAACG‑CGG‑Tttttttttttcccgaa                                                                                                                             >  2:1229185/1‑231 (MQ=255)                                                                                                                            agggtgCAGAATATCCTGGTACATATTCCAATGGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACgg                                                                                                                    <  2:376699/236‑1 (MQ=255)                                                                                                                                               tACATATACCAATAGGGGCAGGCGGTTATTATCTCTTCTTGATCGAAACGGTAGGCACG‑GTCGCCGGGTTTCTGTG‑TATTGGCCTTGCCCGGC‑TGGGGA‑AGATCCGCAACAAC‑CTTTTTTGCCTTTT‑CTACGCCCCCTTGTTTGGCATTATTTTT‑TTTCACTTTC‑ACCTGTTTTCCCCACCGG‑TTTTTCAG‑GGT‑TTTTTTTTTTTCGCGAATTTTTTTGGTTGGTGTTCGTGGtgcgg                                                                                                  >  1:656121/1‑234 (MQ=255)                                                                                                                                                 cATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGGGA‑AGATCAACCACTAC‑TTCTTTGGCCTGAT‑CAACGTCCCCTTGTTTGGCATTATTTTC‑TTTCCGATTC‑AGCTGTATGCCAGCCTGC‑CATTACAG‑GTG‑TTTTTTTTTGCCGCGAATATTTACGGTTGGGATGCGTGGTCGCGGCa                                                                                               >  2:549935/1‑239 (MQ=255)                                                                                                                                                   tataCCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGGAGGCACG‑AACGCCGGGGTGCTGTG‑TATTGGCCTTGCCCGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTGACCTTGTTTGGCATTTTTTTC‑TTTTCGATTC‑AGCTGTATGCCAGCCTGC‑CATTACAG‑GCG‑TTTTTTTTTTCCGCGCATATTTACGGTTGGTATGCGTgggtcggcacaa                                                                                             >  2:10405/1‑229 (MQ=255)                                                                                                                                                    ataCCAATAGGGGCAGGCGGTTATGTACTCTCATGGGTCGGAACGGCAGGGACG‑CTCCCCGGGGTGCTGTG‑TATTTGCCTTGCCAGTC‑TGGAGA‑AAATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAAAGTTCCCTTGGTTGGGCTTATTTTC‑TTTTAGTTTC‑AGCTGTATGCCAGCCTGC‑TTTTACAG‑GGG‑GTTTTTTTTTTCTCGCATTTTTACGGGTGGTGTGctggggtggggtaaac                                                                                            >  2:675707/1‑224 (MQ=255)                                                                                                                                                    ataCCAATAGGGGCAGGCGGTTATGATCTTTCCCGGGTCGAAGCGGTAGGCACG‑ATCGCCGGGGTTCTGTG‑TTTTTGCCTTGCCCATC‑TGGAGA‑AGATCCGCCACTAC‑CTCTTTTGTCCGAT‑CAACCGCCACTTCTGTGGGCTTATTTTC‑TTTCTTCTTC‑TGCTGTCTGGCCCCCCGC‑GCTTACAC‑GTG‑TTTTTCTTTTTCGCGCATATTTttggttgggcttgatgttcaggacaaa                                                                                             >  1:390527/1‑219 (MQ=255)                                                                                                                                                    ataCCAATAGGGGCAGGCGGTTATGATCTCACCAGGATTGAAGCTGTAGGCGCG‑AGCGACGGGGTGCTGTG‑TAATTGGCTTTCTAGTA‑TGGAGA‑AGGTCCGCAACCAC‑ATTTTTGTGCCGGT‑CAACGTGCCCTTTGTTGGTCTTATTTTC‑TTTTGGATTCTAGTCGTA‑GCCAGCACGCACATTAAGG‑GGT‑TTTTTTTTTGCCGCAAATATTTAAGATTggattcggggtggagaaaaac                                                                                             >  1:557035/1‑220 (MQ=255)                                                                                                                                                       ccAAGAGGGGCAGGCGGTTAAGATCTTCCCGGGAACGAAGCTGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGGATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCGTTGCCTCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAg                                                                                         <  2:180622/239‑1 (MQ=255)                                                                                                                                                        gggtataaggggcaggtgttGATCTCTCATGGCTCAGAGCGGGAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCGTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGt                                                                                        <  1:93753/220‑1 (MQ=255)                                                                                                                                                        cAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGGAGGCACG‑ATCGCCGGGGTGCTGTG‑GATTTGGCTTTCCAGTC‑TGGAGA‑AAATCAGCAACTAC‑TTTTTTTGCCTGGT‑CAACGTCACCTTGTTTGGGATTATTTTC‑TTTTAGATTC‑AGCTGTATGCCAGCCTGC‑TCTTACAG‑GGG‑TTTTTCTTTTCCGCGAATATTTACGGTTGGTATTGGTGGTTCGGGCAAACCaac                                                                                        >  2:621242/1‑237 (MQ=255)                                                                                                                                                         aaTAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGGAGGCACG‑ATCGCCGGGTTGCTGGG‑TATTGGCCTTGCCAGTC‑TGGGGA‑AAAACAGCAACTAA‑TTCTTTTGCCTGAT‑CAAAGTCACCCTGTTTGGCATTTTTTTC‑TTTCAGATTT‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GGG‑TTTTTTTTTTCCGCGAATATTTTCGGGTGGTATGGGTGGGTCCGGCAAAACAGTc                                                                                       >  2:889776/1‑239 (MQ=255)                                                                                                                                                            aGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑GATTGGCCTTGCCCGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTTGCCTTAT‑CCACGTCACCTTGTTTGGCCTTATTTTT‑TTTTAGATTT‑AGCTGTATTCCAGCCTCC‑TATTACAG‑GGG‑GTTTTTTTTTCCGCGGATATTTACGGTTGGTATTCCGTGTGGGGGCAAAACAGTCAGa                                                                                    >  2:821255/1‑239 (MQ=255)                                                                                                                                                                gCAGGCGGTTATGATCTCTCCTGGAACGAAGCGGGAGGCAAG‑ATAGGCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGGGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTT‑TTTTAGATTC‑TGCTGTATGCCAGCCTGC‑TATTATCG‑GTG‑TTTTTTTTTGCCGCGAATATTTACAGGTGGGATGCGGGGTGGCGGAAAACACGACAGAAAGa                                                                                >  1:597230/1‑239 (MQ=255)                                                                                                                                                                  aGGCGGTTATGATCTCTCATGGGTCGAAGCGGGAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AAATCAGCCACTAC‑TTCTTTGGCCTTGT‑CAACGTCACCTTGTTTTGCATTATTTTC‑TTTTCGGATC‑AACTGTATGGCCGCCTGC‑TATTTCAA‑GAG‑TTTTTCTTTGCACCGCATATTTTCGGTTGGGATGCGTGGTGGCGCCACACCAGTCCGAAcacgg                                                                              >  2:1078245/1‑235 (MQ=255)                                                                                                                                                                   gggtgtttgtatctttgctGGATTAAAGCGGTTGGCAAG‑ATCGCCGGGTTTCTGTG‑TATTGGCCTTGCCCGTC‑GGGAGA‑AGGTCAGCAACTTC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCAATATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAgg                                                                              <  1:593131/220‑1 (MQ=255)                                                                                                                                                                   ggCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAgg                                                                              <  2:1090926/238‑1 (MQ=255)                                                                                                                                                                     tgttttttttttctttttgATGGAGGCTGTAGGCACG‑GCCGCGGGGTTGCTGTG‑TATTGTCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCCGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCg                                                                            <  1:631095/220‑1 (MQ=255)                                                                                                                                                                                tctcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGgag                                                                                                                                                                                                                                                      >  2:188537/1‑65 (MQ=255)                                                                                                                                                                                tctcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGgag                                                                                                                                                                                                                                                      <  1:188537/65‑1 (MQ=255)                                                                                                                                                                                 ctcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTc                                                               <  2:318147/239‑1 (MQ=255)                                                                                                                                                                                  tcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCcgcg                                                                                                                               <  2:675060/174‑1 (MQ=255)                                                                                                                                                                                  tcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCcgcg                                                                                                                               >  1:675060/1‑174 (MQ=255)                                                                                                                                                                                  tcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGCCAAACCAGTCAGAACGAGGCGGAGTTGAAAATTc                                                               >  2:559957/1‑238 (MQ=255)                                                                                                                                                                                  tcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAgg                                                                              >  2:93753/1‑223 (MQ=255)                                                                                                                                                                                  tcATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACAAGTCAGAACGAGGCGGGGTTGAAAATTc                                                               >  2:1181429/1‑238 (MQ=255)                                                                                                                                                                                        atggAAGGGGTGGGCACG‑ATCGCCGGGGTTGTGGG‑TTGTGGCCTGGCCAGTC‑TGGATA‑AGGCCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTGCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGt                                                         <  1:559957/235‑1 (MQ=255)                                                                                                                                                                                         tCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGC‑AGATCAGCAACTCT‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAAGATTTACGGTTGGTATGCGTGGTCGCGCCAAACCAGTCAGAACGAGGCGGGGTTGAAAATTCGCTGGtt                                                        >  2:76689/1‑238 (MQ=255)                                                                                                                                                                                         tCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCTTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTAGCTGGtt                                                        >  2:1242757/1‑238 (MQ=255)                                                                                                                                                                                            caGCGGTAGGCACG‑TTCGGCGGTTTTCTGTT‑TATTGGCCTGGCCAGTC‑TGGAAA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCa                                                    <  1:781666/238‑1 (MQ=255)                                                                                                                                                                                            aaGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGCCAAACCAGTCAGAACAGGTCGGAGTTGGAAATTCGCCGGTTGCCa                                                    >  2:1079613/1‑239 (MQ=255)                                                                                                                                                                                             gtgggTAGGCCCG‑TTCCCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTTC‑TTCTTTGGCCTGAT‑CAACGTCCCCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCa                                                    <  1:909377/235‑1 (MQ=255)                                                                                                                                                                                             aGCGGTAGGCACG‑ATCGCCGGGTTGCTGGT‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCa                                                    <  2:727986/238‑1 (MQ=255)                                                                                                                                                                                                   aGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGa                                             <  2:597230/239‑1 (MQ=255)                                                                                                                                                                                                    ggtccgcttcgtctggtttttgtt‑TATTGGGCTTGCCGGTC‑TAGGGA‑AGTTCAGCTACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTA‑TTTCAGATTC‑AGCTTTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGGCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGa                                             <  1:621242/219‑1 (MQ=255)                                                                                                                                                                                                    ggCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTTTTTGCCGCGAATATTTACGGTTGGTTGTCGTGGTCGCGACAAACCAGTCAAAACGAGGGGGGTTTGAAAATTCGTGGGGTGCCATTGCCaaa                                            >  1:753488/1‑239 (MQ=255)                                                                                                                                                                                                       aCG‑ATCGCCGGGTTGCTGGG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAgg                                          <  2:656121/238‑1 (MQ=255)                                                                                                                                                                                                        caggtCGCCGGGTTGTGGTG‑TATGGTCCTTTCCAGTC‑TGGAGA‑AGATCAGGAACTCC‑TTCTTTGGCCTGAT‑CAACGGCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCCGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAgg                                          <  2:1006132/234‑1 (MQ=255)                                                                                                                                                                                                                       cTGTT‑TATTGGCCTTGCCAGTC‑TGGATA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTTACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGc                          <  1:852917/239‑1 (MQ=255)                                                                                                                                                                                                                       cTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGCGGCGGGGTTGAAAATTTCCTGGGTGCCATTACCGAAGGAACTCAGCTGGTTgga                          >  1:1067771/1‑238 (MQ=255)                                                                                                                                                                                                                                tGGCCTTGCCAGTT‑TGGGAA‑TGTTTCGCAATTTC‑TTTTTTGGCCTGAT‑CAACGTCCCCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCg                  <  1:685725/239‑1 (MQ=255)                                                                                                                                                                                                                                 ttaagggacaggttaggAGA‑ATCTCAGAAACTGC‑TTCTTTGGCCGGCC‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTTTATGCGCGGCGGC‑TATTACAG‑GTG‑TTTTTCATTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCg                  <  1:1279660/224‑1 (MQ=255)                                                                                                                                                                                                                                  gCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCCGTCAGAACAAGGCGGAGGTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGGTGGCGGttttggtt                >  1:387730/1‑234 (MQ=255)                                                                                                                                                                                                                                    cTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCgttgt              <  2:449235/239‑1 (MQ=255)                                                                                                                                                                                                                                     ttGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCgttgtt             >  2:926204/1‑239 (MQ=255)                                                                                                                                                                                                                                     ttGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCgttgtt             >  2:813177/1‑239 (MQ=255)                                                                                                                                                                                                                                     ttGCCAGCC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGGCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GAG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTAGCCGTGGTCGGGACAAACCAGTCAGAACGAGGCGGGGTGGAAAATTCGATGGGTTCCATTGCCGCAGGAACCTGGCGGGGTGGCGGtgttcggttgt             >  1:910951/1‑229 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCTACTAC‑GTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGTCGCGAATATTTACGGTTGGTATTCGTGGTCGCGCCAATCCATCTCGAACGAGTCGGAGTTGCACGTTTGCTGGTTGCCATCGCCGAAGGCACACAGCTGGTTGGCGGTGTTGGTTTTTTc           >  1:1214120/1‑239 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCAACTTC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACt                                      <  1:1241561/212‑1 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGCGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTc           >  2:165145/1‑239 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACt                                      >  2:1241561/1‑212 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTc           >  2:1135024/1‑239 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTc           >  2:163153/1‑239 (MQ=255)                                                                                                                                                                                                                                       gCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGCCGGAGTTGGAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTGGTTTc           >  2:87110/1‑239 (MQ=255)                                                                                                                                                                                                                                         cAGTT‑TGGAGA‑AAATCAGCAACTTC‑TTCTTTGGCCTGAT‑CAACGGCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGa         <  2:705891/239‑1 (MQ=255)                                                                                                                                                                                                                                         cAGTC‑TTGATA‑AGTCCTGCAATTTT‑TTCTTTGTCCTGCT‑CAACTGACCCTTTTTTGGAATTATTTTC‑TTTGATACTC‑AGTTGTATGCCGACTTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCTGTAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGGCAGAACGAGGCGGATTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTGAGCTGGTTGGCGGTTTGCGTTGTTTCGa         <  1:421455/239‑1 (MQ=255)                                                                                                                                                                                                                                          cgTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGAt        <  2:1067771/238‑1 (MQ=255)                                                                                                                                                                                                                                          cgTC‑TAGAAA‑TGATCAACAATTTA‑TTTTTTTGCCTGTA‑CCAGGTTACCTTGTTTTGCTTTATTTTC‑TTTCAGAGTC‑AGCTGTATGCCAGCCGGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGGCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGAt        <  1:614901/238‑1 (MQ=255)                                                                                                                                                                                                                                          ccttttGGAGA‑ATATCAGCAATTCA‑TTCTTTGGCCTTAT‑CAACGTCACCTTGTTTTGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCCGCCTGC‑TATTACAG‑GTT‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGa         <  1:1242757/234‑1 (MQ=255)                                                                                                                                                                                                                                            tC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACAAGTAAGACCGAGGCGGAGGTGAAAATTCGCTGGTTGCCATTGCCAAAGGCACTAAGATGGTTGGCGGTTTACGTTGTTTCGATTg      >  1:31496/1‑239 (MQ=255)                                                                                                                                                                                                                                                  aga‑agaTCAGCAACTAC‑TTCTTTGGCCTGAT‑TAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGATTGGTCt  <  2:1165707/238‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      ATGAAAGGCTGCGAGAGAGGGCGCTGGTGCCGCTCAATCGTATGCTGGATTTTGCGGCTACACTACGTGGATAACGAATAATAAGGCGTAACGTTACGCTTTGGGGGAAAGATGGATTTTTTTAGTGTGCAGAATATCCTGGTACATATACCAATAGGGGCAGGCGGTTATGATCTCTCATGGATCGAAGCGGTAGGCACG‑ATCGCCGGGTTGCTGTG‑TATTGGCCTTGCCAGTC‑TGGAGA‑AGATCAGCAACTAC‑TTCTTTGGCCTGAT‑CAACGTCACCTTGTTTGGCATTATTTTC‑TTTCAGATTC‑AGCTGTATGCCAGCCTGC‑TATTACAG‑GTG‑TTTTTCTTTGCCGCGAATATTTACGGTTGGTATGCGTGGTCGCGACAAACCAGTCAGAACGAGGCGGAGTTGAAAATTCGCTGGTTGCCATTGCCGAAGGCACTCAGCTGGTTGGCGGTTTGCGTTGTTTCGATTGGTCT  >  NC_000913/783055‑783530

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