Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,965	A→G	9.3%	intergenic (+439/+122)	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,965	0	A→G	9.3%	3.5 ‑0.3 ‑0.0	76	intergenic (+439/+122)	narI/tpr	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/protamine‑like protein
Reads supporting (aligned to +/- strand):  new base (4/5):  ref base (25/40):  total (29/47)
Fisher's exact test strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

AAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT  >  NC_000913/1286666‑1287262                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                          aaaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGAc                                                                                                                                                                                                                                                                                                          <  1:266328/301‑1 (MQ=255)            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:745164/301‑1 (MQ=255)                                                             aTTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                                                                                                                                                                              <  1:643044/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                                                                                                                                                                                                                                                                                           >  2:348712/1‑142 (MQ=16)                                                                                                                                                                                atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                                                                                                                                                                                           <  1:348716/142‑1 (MQ=16)                                                                                                                                                                                  atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                           >  1:93350/1‑252 (MQ=255)                                                                                                                                                                                  atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                                                                           <  2:93347/252‑1 (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)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                          AAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCT‑TCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGC‑AACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT  >  NC_000913/1286666‑1287262

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.