Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,974	C→T	6.7%	intergenic (+448/+113)	narI → / ← tpr	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/protamine‑like protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,286,974	0	C→T	6.7%	2.6 ‑0.4 ‑0.0	77	intergenic (+448/+113)	narI/tpr	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/protamine‑like protein
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (26/42):  total (30/47)
Fisher's exact test strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.01e-01

GGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286677‑1287271                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                          ggAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                                                                                                                                                                                         <  1:604994/300‑1 (MQ=255)    gTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGCGCGGTTCGCCCTTAGGGTCGTTGCCTGCGGC‑AACGCTCTCGCGGTGACGCTCGAGCGGaag                                                                                                                                                                                                                                                                                                      >  1:440142/1‑299 (MQ=255)            ttttGATTATTTAAGGAAGGGATGGTGGTGGGGGAAGGATTTCTCAGCGCTGGGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCg                                                                                                                                                                                                                                                                                              <  1:589471/300‑1 (MQ=255)                gATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGACTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGAAAGCGCGGCGCGGTTCGCGGTATGGGTGGAGGCCGGCGGC‑AAGGCTCTCTCGCTAACGTCCGAGCGGATCTTTAGtagtat                                                                                                                                                                                                                                                                                           >  1:85169/1‑294 (MQ=14)                    attTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGGCGGACCTTAGTCGAAGCTtct                                                                                                                                                                                                                                                                                      >  1:552948/1‑300 (MQ=255)                       tAATGAAGGGATGGTAGTGGGGGAAGGATTACTCATCGCTGCGCGCT‑TCGCCCTTCGGGTCATTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATc                                                                                                                                                                                                                                                                                  <  1:440945/301‑1 (MQ=255)                        aaGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATc                                                                                                                                                                                                                                                                                  <  1:8285/300‑1 (MQ=255)                           gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATCTGTGACTTGTATGCGAATTCCGGGGCTTATGAACGTTACGGTGTTGGGCTAAGGGTTAGTCGGCGCTGTGAGCCTCACCCTTCGGTTCGTTGCCTCGGGC‑CACGCTATCTCGCTCGAGCTCGAGCCGAACCTTGGTAGAAGCTTCACATACtt                                                                                                                                                                                                                                                                               >  2:643020/1‑300 (MQ=14)                            aaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                                                                                                                                                                                                              <  1:968272/300‑1 (MQ=255)                              gCGATGGTGGTGGGGGCAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                                                                                                                                                                                            <  1:450386/300‑1 (MQ=255)                                       gtggGGGAAGGATTACTCAGCGCTGCGCTCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                                                                                                                                                                                   <  1:837806/300‑1 (MQ=255)                                       gtggGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                                                                                                                                                                                   <  1:419173/300‑1 (MQ=255)                                            ggAAGGATTACTCGGCGCTGCG‑GCTGGCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCGCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGTTTCTCATCGGTCCCCGCGTGGGCAGAATATTTGATTGGGGATGCGTATGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                                                                                                                                                                              <  2:783527/300‑1 (MQ=14)                                            ggAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                                                                                                                                                                              <  1:891188/300‑1 (MQ=255)                                                ggATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCCCTCTCGCTGACGCCCCAAGCGAACCTTAGTCGAATCTTCTCATCCTTCCCCGCTTGTGCAGATTATTTg                                                                                                                                                                                                                                                         >  1:897305/1‑301 (MQ=255)                                                 gATTACTCAGCGCTGAGCGCT‑TCGGCGTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTAAAGTCGAACCGTGGTCGAAGCTTGTCATACTTCCCCGAGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGGGATGTTGGTGGGGGAAGGATTATTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGC‑AACTCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                                                                                                                                                                                         <  1:426015/300‑1 (MQ=14)                                                  aTTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                                                                                                                                                                                       <  1:643044/301‑1 (MQ=255)                                                  aTTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                                                                                                                                                                                       <  1:745164/301‑1 (MQ=255)                                                    tACTCAGCGCTGCGCGTG‑GCGCCCTTCGGGTCGTGGCCTGCGGCACAGCTCTCTCGCTGGCGCGCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGTTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACGCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                                                                                                                                                                                                     <  1:802786/301‑1 (MQ=255)                                                     aCTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACCCTCGAGCCGAACCTTAGTCGAAGCCTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTgg                                                                                                                                                                                                                                                    >  1:885379/1‑300 (MQ=255)                                                          gcgcTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGTTGCGCGTTTCGCCCTTCGGGTCGGCGGCGGAGGCAAACGC‑CACTGGTTGGCGCTCGAGCCGAACCTGTGGCGAAGCTACTAATTCTTCCCCCGAAGGGCAGGATGTTGGAttcggggag                                                                                                                                                                                                                                               >  2:539707/1‑294 (MQ=14)                                                                          cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCACGGGCAGAATTTTTAATTTGTGATTCGCTTGAGTGTTcgg                                                                                                                                                                                                                                >  1:930840/1‑299 (MQ=255)                                                                             ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                                                                                                                                                             <  1:189159/301‑1 (MQ=255)                                                                             ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                                                                                                                                                             <  1:797546/301‑1 (MQ=255)                                                                              tCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCACCCTTCCCCGCTTGGGCAGAATATTAGATTGCGGATTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                                                                                                                                                             <  1:168474/300‑1 (MQ=255)                                                                                ggggtggtGCCTGAGGAAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTGGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACt                                                                                                                                                                                                                           <  1:449378/293‑1 (MQ=255)                                                                                          cTGCGGCAACGCTCTCTCGCCGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATCTTTGCTTGCGGAATCGCTTGCGAGTTCAGGGCCTTTTGAAAGTGa                                                                                                                                                                                                                >  1:130194/1‑301 (MQ=255)                                                                                                            ccTGGCGCTCAGGTCCAACCTTGGTCGAAGCTTCTCATCCTGTCCCGCTTGGGCAGAATATTTGTGTGCGGATACGTTTGAGAATTTCGGGGCTTTAGAAAGTGATGGTGGTGGGAGAAGGATTACTCAGCGCTGCGCGCTTCGCCATTCGAGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                                                                                                                              <  1:876856/300‑1 (MQ=14)                                                                                                                    tCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc                                                                                                                                                                                      <  1:539716/301‑1 (MQ=255)                                                                                                                      aaGTCGAACCTTGGTCGAATCGTATCAACCGTCACCGCTGGGGCAGAATATTTGATTGCGGAGTCGTTTGAGAATTTCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    <  1:866050/301‑1 (MQ=255)                                                                                                                      aaGCCGAACCTTGGTCGAAGCGTTGCATCCGTCCACGTGTGGGCAGAATATGTGATTGCGGAGTAGTATGAGGTTTCCGGGGCTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    <  1:930301/301‑1 (MQ=14)                                                                                                                      aaGCCGAACCTTGACCGAAGGTTCGGATCGGTGCACGCTTGGGCAGAGGATTTGGTGTCGGTGTCGTTTGAGAATTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    <  1:353122/301‑1 (MQ=14)                                                                                                                       aGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGGAGAATATTTTATTGCGGATTCGCTTGACAGTTCAGGGACTTTTGAAAGTGGTGGTGGTGGGGGAAGGATTACTCAgag                                                                                                                                                                                     >  1:528300/1‑297 (MQ=255)                                                                                                                            aaCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGTGCAGAATATTTGATTGCGGATTCGTGTGAGAATTCCGGGTGTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                                                                                                                                                                               <  1:591310/300‑1 (MQ=255)                                                                                                                             aCCTTGTTCGAATCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATGCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGTTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                                                                                                                                                             <  2:899888/253‑1 (MQ=255)                                                                                                                             aCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                                                                                                                                                                             >  1:899929/1‑253 (MQ=255)                                                                                                                                          cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCTGACTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGGCCTGCGCGCTTCGCCCTTCggg                                                                                                                                                                >  1:286/1‑301 (MQ=255)                                                                                                                                                 tacTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTTGCCCTTCGGGTAGtgg                                                                                                                                                          >  1:246548/3‑298 (MQ=255)                                                                                                                                                    ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATCTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGTAGGATTACTCCGCACTGCGCTCTTCGCCGTTCTGCTCCATGCCt                                                                                                                                                       >  1:805256/1‑300 (MQ=255)                                                                                                                                                    ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGAGTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGTGCTGCGCGCTTCGCCCTTCGGGTTGTTGCCtt                                                                                                                                                      >  1:823303/1‑300 (MQ=255)                                                                                                                                                        ccGCTTGGGTAGAATATTTGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGc                                                                                                                                                  <  1:884334/301‑1 (MQ=255)                                                                                                                                                             tGGGCAGAATATTTGCTTGCGGATTGGGTTGAGAATTCCGGGGCTTTGGAAAGTGATGGTGGTGGGGGAAGGATTACGCCGCGCTGGGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                                                                                             <  1:577889/301‑1 (MQ=255)                                                                                                                                                              gggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGAATCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGGTTCGCCGTTCAGTTCGTTGCCTCCGGCAACGc                                                                                                                                             >  1:643851/1‑300 (MQ=255)                                                                                                                                                                   gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATGTCGGATTCGCTTGTGGGTTCATGCACGTTATCAACTGTTGGTGGTTGTGGACGGTTTACTCATCGCTGCACAATCCGCCGTACGCGTGGTCTCCTGCTGCATCActcttt                                                                                                                                        >  2:896571/1‑298 (MQ=14)                                                                                                                                                                   gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATACGCATCAGACTTTGTCCCGGATTCGCTTGGTGTTTCCGTGGCTGTTCAAAGTGTTGGTGGTGGCAGAAGGATTTCTTCACGCTTCGCGCTTCTCCTTGTGGGTCGttcacttcgctgacgcctct                                                                                                                                         >  2:677997/1‑281 (MQ=14)                                                                                                                                                                     atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                                                                                                                                                                                                    <  1:348716/142‑1 (MQ=16)                                                                                                                                                                     atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                                                                                                                                                                                                    >  2:348712/1‑142 (MQ=16)                                                                                                                                                                       atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    <  2:93347/252‑1 (MQ=255)                                                                                                                                                                       atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    >  1:93350/1‑252 (MQ=255)                                                                                                                                                                                                 tttgtcgtctgtggaaGGGATGGTGGTGTGGTAAGGAGTACTCCTCGGTGAGCGCTTCTCCCGTCGGGTCGGTGCGTGCTGA‑AACACTCTCTGGCGGCCGCTCGAGACGAACCTTAGTCGATGCTTATCATCATATGTCGAATAGGCAGAGTATTAGATTGCGGAGTAGCTTGAGGGTTCGGTGATTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                                                                                                         <  2:612235/287‑1 (MQ=9)                                                                                                                                                                                                          tttGAAAGGGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTCCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  1:678021/301‑1 (MQ=255)                                                                                                                                                                                                                 gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                                                                                                                                                                                                     >  1:598344/1‑145 (MQ=38)                                                                                                                                                                                                                 gTGATGGGGGTGGGTGAAGGATTCCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                                                                                                                                                                                                     <  2:598329/145‑1 (MQ=255)                                                                                                                                                                                                                       gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCCCCCTTCGGTTCGTTGCCTGCCGCCACGCTCTCTCGCTGGCGCTCGCATCGAACCTTGGTTTAACCTTCTCATCCTTCCCCGCATgg                                                                                   >  1:613612/1‑301 (MQ=255)                                                                                                                                                                                                                          gtggGGGAGAGATTACTCGGCGCTGTGCCCTTCGCCCTTCGGGTCGTGGCCTGCGGC‑AACGCTCCCGCGCTGACGCTCCAGGCGAACCTAGGTGGAAGCTTCTCACCCGTACCCGTATTGGCAGAAGATTTGATTGCGGCTTCGCTTGAGAGTTCAAGTATTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                 <  1:152174/300‑1 (MQ=14)                                                                                                                                                                                                                          gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGACCTTCGGGTCGTTGCCTGCGGC‑ATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGc                                                                                 <  1:579204/300‑1 (MQ=255)                                                                                                                                                                                                                               ggAAGGATGACGCAGCGCTGCGCGCTTGGCCCTTCGGGTCGTGGCATGCGGC‑AACGCTCGCGCGCTGACGCTCGAGGCGAACCTGGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCACAAGATTTGATTGCGGATTCGGTTGAGAGTTGAGAGACTATTGAAAGTGATGGTAGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                            <  1:176978/300‑1 (MQ=14)                                                                                                                                                                                                                               acctgggttACTCTGCGCTGCTCGCCTCGCCCTTCGGGTCTTTGCTCGTGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           <  2:774426/294‑1 (MQ=255)                                                                                                                                                                                                                                   ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGCata                                                                           >  1:774460/1‑297 (MQ=255)                                                                                                                                                                                                                                   ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCCCATGGGCAGAATAtttt                                                                       >  1:946136/1‑300 (MQ=255)                                                                                                                                                                                                                                    gATTACTCGGGGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTa                                                                       <  1:872713/300‑1 (MQ=255)                                                                                                                                                                                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTa                                                                       <  1:896612/300‑1 (MQ=255)                                                                                                                                                                                                                                     atatgtatgcccttggggcttgctcttgtgggTCGGTGACTGCGGC‑GATGCTCTCTCGCTGGCGCTCCGGCCGAGTCTTAGTCGAATCTTCTCATTCGTCCCCGCAGTGGCAGAAGACTTGATTGGCTCTTCGTTTGAGAGTTTAGGTACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAt                                                                     <  2:554300/272‑1 (MQ=14)                                                                                                                                                                                                                                     aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAt                                                                     <  1:525030/301‑1 (MQ=255)                                                                                                                                                                                                                                      tgcacaggggtttcgcggtcagCCCTCTGGGCCGTTGCCTGCGGC‑AACGCTCTGTCGCTGGCGCTCAAGCCGAACCTTGGTCGAAGTTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTGCCGGGACTTTTGAAAGTGATGGTGGTGGGCGAAGGAGTACTCAGCGCAGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAt                                                                     <  2:930792/279‑1 (MQ=11)                                                                                                                                                                                                                                       tACTCAGCGCAGCGCTGTGCGCCCGTCGGGTCGTTTCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGATGCTTCTCATCCTTCCCCTCATGGGCAGATTATTTGATTGCGGATTCAATTGATAGTTCAGGGTCTTTTGACAGTGATGGTGTGGGTGGAAGGATTACTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAt                                                                     <  1:630188/299‑1 (MQ=14)                                                                                                                                                                                                                                        cagatgtggatcggagacggcctgtCGGGGCGTGGCCAGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGATTCGCGTGAGAGTGCAGGGACTGTTGAAAGTGGTGGTGGTCGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                   <  2:491108/276‑1 (MQ=14)                                                                                                                                                                                                                                        aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    <  2:522893/187‑1 (MQ=255)                                                                                                                                                                                                                                        aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                                    >  1:522902/1‑187 (MQ=255)                                                                                                                                                                                                                                                            cccTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGCGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTACGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCCGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGCATATTTAATTGctgtcccttataccactctc                                               >  1:491116/1‑281 (MQ=14)                                                                                                                                                                                                                                                                 cGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAGGTTCAggggat                                         >  1:623603/1‑298 (MQ=255)                                                                                                                                                                                                                                                                         ttccTGAGGC‑ATCGCTCTCTCGCTGGCGCTCGAGCCGGACCTTGGTCTAAGCTTCTCATCCTTCTCCGCATGTGCAGAATATTTTATTGCCGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGTGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 <  1:707190/299‑1 (MQ=14)                                                                                                                                                                                                                                                                                  c‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTTAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                        <  1:957248/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTGGAAAGTGATGGTGGTggg                    >  1:23975/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGCAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           >  1:551149/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                      cTCGAGTCGAACCTTAGTCGAAGCGTCTCGTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATACGCTTGAGAGTTCATGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAcc     <  1:400964/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                          aGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGGTTCGCGGTTCGGGTCGTGGCCTGCGGCCACGCTCTCTCGCTGGCGCCAGAGTCGGAACTTGGTCGACGCTTCTCATCCTTCCCCGCCTGGGCAGAATAGTTAAGTGAGGGTTAGTTGGGAAGTTCAGGGCCTTGGTATAGGTGGGGTGGTGGGGAGAGGATTCGACCTTTc  >  2:158582/1‑300 (MQ=17)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                          GGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286677‑1287271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.