Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,542,110	A→G	22.6%	S161P (TCC→CCC)	narZ ←	nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,542,110	0	A	G	22.6%	140.6 / 7.1	82	S161P (TCC→CCC)	narZ	nitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (36/24);  new base G (2/18);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 8.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGG  >  NC_000913/1541872‑1542339                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         tCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGa                                                                                                                                                                                                                                          <  1:1024390/239‑1 (MQ=255) tCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCCGGGTTAGTTTCT‑GCCCGTTGGa                                                                                                                                                                                                                                          >  1:27296/1‑239 (MQ=255)             gTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGGTCCT‑GCCCGTTTGA‑GC‑GGATCAA‑ccc                                                                                                                                                                                                                           >  2:454895/1‑239 (MQ=255)                  cTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGGTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑TCCCGTTGGA‑GC‑GGATAAA‑Cccgcccc                                                                                                                                                                                                                      >  2:963041/1‑237 (MQ=255)                    ccccccAGGCCGTCGGGGAGGCGGGCGGCAGGC‑CGCAATACCAGACGTAGAAACTTAAACAGG‑CGCCGCCAGGCCGCGAC‑A‑AATAACGCGTGCGGGCGGCGTAAGAACCCGCCGACCTCGCGGGGATCGGAGAGAAACC‑GG‑CAACGCGCCCCGGGCCGTAGTTTTGGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCAAGTGGGA‑GC‑GGATAAA‑CCCGCCgcgc                                                                                                                                                                                                                    <  1:1305505/237‑1 (MQ=255)                       cccAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTCCCGGCGGCGTAAGAACCCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCa                                                                                                                                                                                                                 <  2:689983/239‑1 (MQ=255)                         acggTTATCGGCGAGGCGGGCGG‑AGGTCCGCATAACCAGTCGGAGAAACTAAAAAGGG‑GCCCGCCAAGCAGCGAC‑A‑AATAACGCTTCCCGGCGGCTTAAGAACCCATCGACATCGCCGGGATCGGCGAAAACCC‑GG‑AACCGCGACCCGGCCCGTGGTTTTTGAGGGCCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGG‑GC‑GGAAAAA‑CCCGCCGCGCCCACg                                                                                                                                                                                                               <  1:976320/237‑1 (MQ=255)                          aGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGGTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGc                                                                                                                                                                                                              >  2:773113/1‑239 (MQ=255)                            gtacggggggAGGCGGGCGGGAGGC‑CGCAATCCCGGTCGTAGAAACTTAACCAGG‑TCCCCCCAGGAGGCGAC‑A‑GATACCGCGTCCCGCGGGCGTAAAAAACCAGCGACAGCGCGGGGAGCGGCGAGAAACC‑GG‑CACCGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACg                                                                                                                                                                                                               <  1:1713495/227‑1 (MQ=255)                            gcatgcggcggGGCGGGGGGAAGGT‑CGAAATCCCAGTCGGAGAACCTTAAACAGG‑TCCCGCCAGGCAGCGAC‑A‑GAGACCGCGTTCCGGCGGCGTAAGAACCCATCGACATCGCCGGGATCGGCGAGAACCC‑GG‑CAACGCGATCGGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCG‑GCCAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCa                                                                                                                                                                                                             <  1:1690445/233‑1 (MQ=255)                             tcgccggcggGGGGGGCGGCAGGT‑CCAAATCCCAGGCGTAGAAACTTAAACGGG‑TGCCCCAAACCGGCGAC‑A‑GAAACCGCGTCCGGGCGGCGAAAGAACCCAGCGACACCGCGGGGACCGGCGAGAACCC‑GG‑CACCGCGATCCGGCCGGAAGTTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTGGGG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCAc                                                                                                                                                                                                            <  1:957226/234‑1 (MQ=255)                             tcacggggggagggggcggGAGGT‑CGAAATACCAGTCATAAAAAATGAGCCAGG‑GGCCGCCAAGCGGCGGG‑A‑GATAAGGCGTTACGGGTGCGTAAGAAACCATCGACATCGCCGGGATGGGCGAGAACCC‑GG‑CAACGCGATCCGGGCCGTAGGTTTGGATGGGCCAGGCGTTAGCGGC‑GGCAATCAGGTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCAc                                                                                                                                                                                                            <  1:699902/225‑1 (MQ=255)                              cATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAAAGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACtt                                                                                                                                                                                                          >  2:910520/1‑239 (MQ=255)                                tCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGGTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACttgt                                                                                                                                                                                                        >  2:1036023/1‑239 (MQ=255)                                 cGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCCCttgtt                                                                                                                                                                                                       >  2:976320/1‑239 (MQ=255)                                 cGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGGTTTGATGGTCCAGACGTTAGCGGC‑GGCCATCCGCCGGTTTTAGTCCT‑TGCAGTTTGG‑GC‑GGGTAAA‑CCCCCCCCGCCCACCCCCctgtt                                                                                                                                                                                                       >  1:1143216/1‑239 (MQ=255)                                       ggcgggcgGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGGTCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAg                                                                                                                                                                                                  >  2:641836/1‑238 (MQ=255)                                                ggcgggacgcaaaaccgtcggtaaaaCTAAAACAGGTTGCCGCCAGGC‑GGGACAA‑GCTAACGCGTTGCGGAGGCATAAAAAACCATCAACATCGGCGGGATGGGCGAAAAACC‑GG‑AACCGCGATCCAGGTCGTAGGTTTTGTTGGTCCTGACGTTAGCGCC‑GGCAATAAGCTGGTTTATTTCGT‑GCCAGGGGGA‑GG‑GGATAAA‑CGCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGc                                                                                                                                                                                           <  1:1208353/215‑1 (MQ=255)                                                 aGGT‑CGCAATACCAGTCGTAAAAACTAAAACAGG‑TCCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAACCCATCGACATCGCCGGGATCGGCGAGAACCC‑GG‑CAACGCGTTCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACt                                                                                                                                                                                        <  2:1332910/238‑1 (MQ=255)                                                   gT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑g                                                                                                                                                                                                                                        <  1:728465/190‑1 (MQ=255)                                                   gT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑g                                                                                                                                                                                                                                        >  2:728465/1‑190 (MQ=255)                                                         aatccaggccgagaaaaCTAAAAGAGG‑TCCCGCCAGCCAG‑GAC‑AGGAAAGCCGGTCCCGGCGGGGAAAGAAACCACCGACATCGCGGGGTGGGGGGAAAAACC‑GG‑CACCGCGATCGGGGCGGTAGTTTTGGATGGCCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGa                                                                                                                                                                                <  1:869707/225‑1 (MQ=255)                                                         aaaccgagtgggaaaactatAAAAAGG‑GGCCGCCAACCAGGGAC‑A‑CATACCGTGTCGCGGGGGGGAAAAAAACCATCGCCATGGCGGGGAGGGGAGAAAAACG‑GG‑GACCGCGATCCGGGCCGTAGGTTTTGAGGGTCCGGACGTTAGCGGC‑GGCAATCAGCTGTTTTATTTCCT‑GCCAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGa                                                                                                                                                                                <  1:1420685/220‑1 (MQ=255)                                                           tACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCACTCCGCTGGTTTAGTTCCT‑GCCAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCCCTTGCTTGTAGCCCAGGCAATTTTGCCGATc                                                                                                                                                                              >  1:1103639/1‑239 (MQ=255)                                                                gTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCa                                                                                                                                                                         >  2:640418/1‑239 (MQ=255)                                                                gTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCa                                                                                                                                                                         >  2:1100658/1‑239 (MQ=255)                                                                   gTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGATTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTCCTGATCGTTCATaa                                                                                                                                                                      >  1:1715532/1‑239 (MQ=255)                                                                     aGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATc                                                                                                                                                                    >  2:1122435/1‑239 (MQ=255)                                                                      caaccaaaaaCGGG‑CGCCGCCAGCCGGCGAC‑A‑GAAAACGCGCTCCGGGGGCGTAAGAACCCCCCTCCACCGCGGGGACCGGCGAAAACCC‑GG‑CAACGCAATCGGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCGGGTTTAGTCCCT‑GCCAGTTGGC‑GG‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCg                                                                                                                                                                   <  1:1243768/232‑1 (MQ=255)                                                                        aaCTAAACCCGG‑TCCCGCCAGGCAGGAAC‑A‑GATAACGCGTGCCGGCGGCGTAAGAACCCACCGACATCGCCGGGGCCGGCGAGAACCC‑GG‑CAACGCGATCCGGGCCGTAGTTTTTGATGGCCCAGCCGTTAGCGGC‑GGCAATCGGCTGGTTTAGTTCCT‑GCCATTTGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCg                                                                                                                                                                   <  1:1600823/237‑1 (MQ=255)                                                                           ttAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACtt                                                                                                                                                                                                          >  2:161314/1‑195 (MQ=255)                                                                           gtaaccggg‑TGCCGCCAAGCGGCGCC‑A‑AAAAACGCGTTCCGGGGGCTAAAGAAACCATCGACAGCGCGGGGATGGGCGAGAAACG‑GG‑CACCGCGTCCGGGGCGGTAGGTTTTGAGGGTCCAGCCGTGGGCGGC‑GGCAGTCGGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAC‑CCCGCGGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGcc                                                                                                                                                              <  1:910520/232‑1 (MQ=255)                                                                              aaCAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTTCGGATCGGTCATAATCCATGCCCa                                                                                                                                                            >  1:121721/1‑238 (MQ=255)                                                                                     tGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGt                                                                                                                                                     >  2:780492/1‑239 (MQ=255)                                                                                     tGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGt                                                                                                                                                     >  1:555000/1‑239 (MQ=255)                                                                                     ggccccaaaGCAGCGAA‑A‑AAAACCGCTTCCCGGCGGGGAAAAAACCCACCACCACCGCGGGGACGGGCGAAAAACC‑GC‑CAACGCGACCCGGGCCGTAGTTTTTGAGGGCCAGGACGTTAGCGGC‑GGCAATCAGCTGTTTTAGTCCCT‑GCCAGTTGGG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGc                                                                                                                                                         <  1:1122435/228‑1 (MQ=255)                                                                                      gccgccAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGGTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGGTCGTTCATAATCTATGCCCACGCCCGt                                                                                                                                                     >  1:1151829/1‑238 (MQ=255)                                                                                              gcagcGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCCGTACCGGATc                                                                                                                                             >  1:1721946/1‑238 (MQ=255)                                                                                               cagcGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGC‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGGTTGTAGCCCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGCTc                                                                                                                                             >  1:920960/1‑237 (MQ=255)                                                                                                   gAC‑A‑GATAACGCGTTCCGGCGGCGAAAAAAACCTGCGGCCTGGGCGGGACGGGGGAGAACCC‑GG‑GCACGGGACCGGGGCCGTAGGTTTGGATGGTCCAGCCGTTAGCGGC‑GGCAACCACCTGTTTTAGTTCCT‑GCCAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCtgt                                                                                                                                        <  1:1570700/238‑1 (MQ=255)                                                                                                    aC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTa                                                                                                                                                    >  2:1712518/1‑225 (MQ=255)                                                                                                              cgcgTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCCCTTGTTTGTAGCTCAGGCCCTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTGGCTTGAg                                                                                                                              >  2:39534/1‑239 (MQ=255)                                                                                                                 gTTCCGGCGGCGTAAGGAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCCATCAGCAGGTTTAGTTCCT‑GACAGTTGGA‑AC‑GGATAAA‑CCCGCCGCCCCCACGCACTTGTTTGTAGCACCGGCACTTTTTCGGTTCGATCTTAATCGCTGCCCACGCCAATACCGGATCGCGGTGTGCCTTGAgg                                                                                                                             >  1:170036/1‑237 (MQ=255)                                                                                                                  ttCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCt                                                                                                                          >  1:586422/1‑239 (MQ=255)                                                                                                                   tCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGAAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGTTCGTTCATAATCGCTGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGTGGGCt                                                                                                                          >  1:994722/1‑238 (MQ=255)                                                                                                                       gcggcgTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTcgc                                                                                                                      >  2:1255057/1‑238 (MQ=255)                                                                                                                          ggaaagaaaaCCGTCGCCATCGCGGGGGGCGGGGAGAAACG‑GG‑AACCGCGTCCGGGGCCGTTGTTTTTGATGGCCCAGGCGTGAGCGGC‑GGAAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGcc                                                                                                                   <  1:1176277/231‑1 (MQ=255)                                                                                                                             tAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCCTAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAg                                                                                                               >  1:858389/1‑239 (MQ=255)                                                                                                                                     cATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCCGTACCGGGTCGCTGTGTTTCTTGAGGGCTTTGCGCCACAGTTCAATCa                                                                                                       >  1:1447926/1‑239 (MQ=255)                                                                                                                                     cATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACCCCAGTACCGGATCGCTGTGTTg                                                                                                                                    >  2:615938/1‑210 (MQ=255)                                                                                                                                      tccGACATCGCGGGAATGGGCGAGAAACC‑GG‑AAACGGGATCGGGGCCGTGGGTTTTGACGGTCCAGACGTAAGCGGC‑GGCAATCAGTTGGTTTAGTTCCG‑GCCAGTGGGA‑GG‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAg                                                                                                      <  2:68100/237‑1 (MQ=255)                                                                                                                                      gcccgcacccgcgggtgcGGCGAAAACCG‑GG‑CAACGCGTCCGGGGCGGTGGGTTTTG‑TTGCCCAAGCGTTAGGGGCGGGCAATCAGTTGTTTTAGTTCCT‑CCAATTTGAG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAg                                                                                                      <  1:1100658/222‑1 (MQ=255)                                                                                                                                        cGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTc                                                                                                    >  2:1725547/1‑239 (MQ=255)                                                                                                                                              cGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGa                                                                                                                                               >  2:982825/1‑190 (MQ=255)                                                                                                                                                         ggggaaaaCCGGG‑CAACGCGATCCGGGCGGTAGTTTTTGAGGGCCCGGAGTTTGGCGGC‑GGCAATCAGCTGTTTTAGTTCCT‑GCCAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCgg                                                                                    <  2:273035/234‑1 (MQ=255)                                                                                                                                                          ccggaAACC‑GGCCACCGGGACCCGGGCGGTGGGTGTTGAGGGCCCAGCGGTTAGGGGC‑GGCACTGGGCTGGTTTAGTCCCT‑GCCAGTGGGG‑GC‑GAATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCgg                                                                                    <  1:456930/235‑1 (MQ=255)                                                                                                                                                                        cacggtacccGGGCGGTGGGTTTTGAGGGCCCAGGCGTTAGCGGC‑GGCAATCGGCTGTTTTAGTTCCG‑GCCAGTTGGG‑GC‑GGAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGt                                                                      <  1:780492/231‑1 (MQ=255)                                                                                                                                                                           gcgATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCGGCTGGTTTATTTCCT‑GCCAGTTGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAg                                                                   <  2:1103639/239‑1 (MQ=255)                                                                                                                                                                               tCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTCCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAgcgc                                                                <  2:624034/238‑1 (MQ=255)                                                                                                                                                                               gccGGGCTGTAGGTTTGGAGGGGCAAGCCGTGGGCGGC‑GGCACTCAGCGGGTTTATTTCCG‑GCTGGTGGGG‑GC‑GAAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCt                                                               <  1:1111248/238‑1 (MQ=255)                                                                                                                                                                               ccgggtctgTAGTTTTGGATGGCCCAGACGTTAGCGGC‑GCCACTCAGCTTGTTTGGTTCCT‑CCTAGTTGGG‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCt                                                               <  1:459661/231‑1 (MQ=255)                                                                                                                                                                               cccGGGCGTTAGGTTTTGATGGTCCGAACGTTAGCGGC‑GCCAACCGGCT‑GTTTACGTCCTGGCCAGTTGGA‑GC‑GAATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAgcgc                                                                <  2:843621/237‑1 (MQ=255)                                                                                                                                                                                     cggtggttttgGAGGGCCCAGACGTTAGCGGG‑GGCAATCGGTTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTaaa                                                          <  2:27296/231‑1 (MQ=255)                                                                                                                                                                                     ccGTGGGTTTTGATGGTCCAGCCGTTAGCGGC‑GGCAATAAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GG‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTaaa                                                          <  2:940026/238‑1 (MQ=255)                                                                                                                                                                                     ccGTGGGTTTGGATGGTCAAGGCGTTAGCGGC‑GACAATCAGCTGGTTTAGTGCCT‑GCCAGTTGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTaaa                                                          <  2:1544546/238‑1 (MQ=255)                                                                                                                                                                                      ggggggTTTTGGTGTTCAAGCCGTTAGCGGG‑GGCAATGAGGTGGTTGAGTTCCG‑GCAGGTGGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGa                                                        <  2:1517856/235‑1 (MQ=255)                                                                                                                                                                                         aGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATAc                                                     <  2:555000/239‑1 (MQ=255)                                                                                                                                                                                         aGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGAGCTGTAAAGATa                                                      >  2:737393/1‑238 (MQ=255)                                                                                                                                                                                         aGGTTTTGAGGGCCAAGACGTTAGCGGC‑GGCAACCAGCTGGTTTAGTTCCT‑CCCAGTTGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATa                                                      <  2:1143216/238‑1 (MQ=255)                                                                                                                                                                                           gTTTTGATGGTCGACACGTAAGGGGG‑GAGAAGCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAgcg                                                                 <  1:355764/225‑1 (MQ=255)                                                                                                                                                                                            tttgggtGGTCCAGACGTTAGGGGC‑GGCAATCAGTTGGTTTATTCCCT‑GCCAGTTGGA‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCa                                                   <  1:632340/232‑1 (MQ=255)                                                                                                                                                                                                     tCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTATGGCTGTAAAGATACCAGGGGTAACtt                                         >  1:1510957/1‑239 (MQ=255)                                                                                                                                                                                                     tCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGAt                                                       >  2:355764/1‑225 (MQ=255)                                                                                                                                                                                                     tCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACtt                                         >  2:632340/1‑239 (MQ=255)                                                                                                                                                                                                               aGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGCTCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTGTTTCAGGCGGTTAGCGCTGTAACGATACCTGGAGTAACTTGCGcccctgc                               >  1:822271/1‑234 (MQ=255)                                                                                                                                                                                                               aGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGAACTTGTTTGTAGCTCAGGCACTTTTGCGGAACGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGATTCGCGCCACAGTTCACTCAGTCGGTTACGGATGAGCGGGTCTTTCAGGCGGTTAGCGCTGGAAAGAAACCAGGAGTAAATTGCGCGCcgcg                               >  1:1532861/1‑239 (MQ=255)                                                                                                                                                                                                                          aTCAGTTGGTTTGGTTCCG‑GCCAGTTGGG‑GC‑‑GATAAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCCCTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                     <  2:1090988/239‑1 (MQ=255)                                                                                                                                                                                                                          aTCAGCTGGTTTGGTTCCG‑TCCAGTTGGG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTGGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                     <  1:1036023/239‑1 (MQ=255)                                                                                                                                                                                                                            cAGCTGGTTTATTCCCG‑GCGGTTTGGG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCcgcg                     <  1:524778/237‑1 (MQ=255)                                                                                                                                                                                                                                 ctggttgggtcgcggcAGGTGGAGGC‑GGAAAAC‑CCCGCCGCGCCCACGCACTGGTTTGTAGCTCAGGCAGTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGt                 <  1:932853/223‑1 (MQ=255)                                                                                                                                                                                                                                   tttAGTTCCT‑GCCAGGTGGA‑GC‑GAAAAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCAt             <  2:1673381/238‑1 (MQ=255)                                                                                                                                                                                                                                    ttAGTGCCG‑GCCAGTTGGG‑GC‑GGATAAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGa           <  1:640418/239‑1 (MQ=255)                                                                                                                                                                                                                                    tatgTTCCT‑G‑CTGTTGGA‑ACAGGATCAC‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATGCCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGa           <  1:592525/236‑1 (MQ=255)                                                                                                                                                                                                                                              gCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAgg  >  2:592525/1‑239 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         TCCGGTACGTCGGTTTGCTCGCCCCAGGTCATCGGCGAGGCGGGCGGCAGGT‑CGCAATACCAGTCGTAGAAACTTAAACAGG‑TGCCGCCAAGCAGCGAC‑A‑GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACC‑GG‑CAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGC‑GGCAATCAGCTGGTTTAGTTCCT‑GCCAGTTGGA‑GC‑GGATAAA‑CCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTTCAGGCGGTTAGCGCTGTAAAGATACCAGGAGTAACTTGCGCCACGCGGGCAGCCGCGAGGTTCATGATTGGGCAGG  >  NC_000913/1541872‑1542339

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