BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,825,265	0	C→T	10.0%	14.1 ‑44.6 0.0	560	E121K (GAA→AAA)	spy	periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand): new base (56/0): ref base (212/292): total (268/292)
Fisher's exact test strand distribution p-value = 4.21e-20
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 564 total aligned reads displayed.

TAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGACCGAACTTGCCTTTGTGGTG  >  NC_000913/1825017‑1825511

  gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGGATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTTCTTTTTCTGTTTCGGCGTCAGGAATTTGTAAATTTTGTTTTGGGTTTTCAtg                                                                                                                                                                                                                                                    >  1:1100950/1‑251 (MQ=255)

  gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGAATATTCAGCAGGTGCAGGCATTTTAACTTTTGCCGCTGGACGTTCTGGCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTTTGGGGTTCCAtg                                                                                                                                                                                                                                                    >  1:260383/1‑251 (MQ=255)

                     tacattttttgaagaagaaaaAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCtt                                                                                                                                                                                                                                 <  1:3255912/239‑1 (MQ=255)

                          ttttgagaaaagaaaaAAGAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgcgc                                                                                                                                                                                                                            <  1:5029502/242‑1 (MQ=255)

                                    aaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCa                                                                                                                                                                                                                  <  2:3926037/247‑1 (MQ=255)

                                     gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTTTTcc                                                                                                                                                                                                                   >  2:5392372/1‑249 (MQ=255)

                                       aaaaaaaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgcgc                                                                                                                                                                                                                            <  1:5585496/237‑1 (MQ=255)

                                                 ccgccgATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGGTTATTCCGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTTAGACGCTTCTCAAAATTAACATTTAAATTCCTTTTCTGTTCCCGCGGCCGGATGGTGTTAATTTTTTTCTTGGTTTTCATGTGCGCCCGCATGTTAGCTTTGCGCTGTTCTTTCCTTTTTTCgtt                                                                                                                                                                                                        >  1:2256869/1‑246 (MQ=255)

                                                       aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCCGCAGTTGCAGGCCTTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAATAGCATTTAATTTCTTTTTCTTTTCCGGCGGCCGGATGTTGTAAATTTTGTTCTGGGTTTTCCTTTGCTCCAGCCTGTTAGCTTTTCGCTGTTCTTCCCTTTTTGCGATCTGCGCtt                                                                                                                                                                                                >  1:132926/1‑250 (MQ=255)

                                                       aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCCGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTTAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTTCGATCTGCGCtt                                                                                                                                                                                                >  1:3102176/1‑250 (MQ=255)

                                                       aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCCGGCAATTTACCTTTTGCCGCTGGACGTTCTTTCCGACGCTTCTCAAAATTAGCATTTAATTGGTTTTTCTTTTCCGGCGTCCGGGTGGTTTAAAATTTTTTCTTGGGTTTCAtg                                                                                                                                                                                                                                                    >  1:5712758/1‑198 (MQ=255)

                                                       aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGcc                                                                                                                                                                                                                                              <  2:2963467/204‑1 (MQ=255)

                                                       aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAAATTCTTTTTTTGTTCCGGCGTCCGGGTGTTGTAAATTTTGTTTTGGGGTTCCATGTGCGCCAGCATGTTATGTTTGCGCGTTTCTTTCATTTTTGCGATCTGCGCtt                                                                                                                                                                                                >  1:5663042/1‑250 (MQ=255)

                                                            gTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTCCTGCGTCAGGATGTTGTAAATTTTGTTTTGGGGTTTCCTGTGCGCCCGCCTGTTAGTCTTTCGCTGTTCTTCCa                                                                                                                                                                                                                  >  2:5976818/1‑227 (MQ=255)

                                                            gTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGGTGTTGTAAATTTTGTTCTGGGTTTTCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTTCGATCTgggc                                                                                                                                                                                                  >  2:4341814/1‑243 (MQ=255)

                                                            gTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTTAATTGCTTTTTTTTTTCCGGCGTCAGGATGTTGTAAATTTTTTTTTTGGTTTTCATTTGCGCCAGCATGTTAGCTTTTCGCTGTTCTTCCATTTTTTTGATTTTCGCTTCAGCtt                                                                                                                                                                                          >  2:3651314/1‑251 (MQ=255)

                                                            gTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTgcgc                                                                                                                                                                                                  <  1:4342576/243‑1 (MQ=255)

                                                              ccACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTTTTCCGGGGTCAGGGTGTTGTAAATTTTTTTTTTGGTTTTCATTTTGGCCAGCATGTTTGc                                                                                                                                                                                                                                   >  2:3026287/1‑208 (MQ=255)

                                                              ccACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGc                                                                                                                                                                                                                                   <  1:3026781/208‑1 (MQ=255)

                                                                 ccGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCCGTTGCAGGCCTTTTACCTTTTTCCGCTGGACGTTTTGTCAGACGCTTCTCCAAATTAGCATTTAATTGCTTTTTTTGTTTCCGCGTCCGGGTGTTGTAAATTTTGTTTTTGGTTTTCATTTGCGCCAGCATGTTTGCTTTGCGCTGTTTTTTCATTTTTGCGGTCTGTGCTTTAGCTTTTatt                                                                                                                                                                                     >  2:2600459/1‑249 (MQ=255)

                                                                 ccGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACt                                                                                                                                                                                     <  2:5819494/251‑1 (MQ=255)

                                                                        aCTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGGATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTCCGGCGTCAGGATGTTGTAAATTTTTTTCTTGGGTTTCATGTTCGCCCGCATGTTAGCTTTTCGCTGTTTTTTCATTTTTTCGATCTGCGCTTCCGCTTTTTCTTTATCGa                                                                                                                                                                              >  2:2251801/1‑251 (MQ=255)

                                                                         cTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTGCCGCTGGACGTTTTGTCAGACGCCTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTTCGGCGGCAGGGTGTTGGAAATTTTGTTCTGGGGTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCCTTTTTTCGATCTGCGCTTCAGCTTTTTCTTTTTCGaa                                                                                                                                                                             >  1:2408498/1‑251 (MQ=255)

                                                                           gCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCAttt                                                                                                                                                                                                               <  1:2154353/215‑1 (MQ=255)

                                                                            cAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGt                                                                                                                                                                          <  1:3728081/251‑1 (MQ=255)

                                                                                gTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTTAAAAATTATTTAGCAGTTGCAGGCATTTTAACTTTTTCCGCTGGACGGTCTGTCAGAAGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTTTTTCGGCGTCAGGATGTTGTAAATTTTGTTTTTGGTTTTCATTTGCGCCAGCATGTTAGc                                                                                                                                                                                                                                   >  2:3864393/1‑190 (MQ=255)

                                                                                gTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTTCTTTTTTTGTTTCCGCGTCCGGGTGTTGTAAATTTTTTTTTGGGTTTCCATTTGCGCCCGCATGTTAGCTTTGCGCTTTTCTTTCATTTTTTCGGTCTTCGCTTCAGCTTTTTCTTTTTCGAAAGTATCg                                                                                                                                                                      >  2:3289675/1‑251 (MQ=255)

                                                                                gTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTTTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTCCCGCGTCAGGGTGTTGTAAATTTTTTTCTTGGTTTCCATGTTCGCCAACATGTTTGCTTTGCGCTGTTTTTCCATTTTTTCGATCTGCGCTTCAGCTTTTTCTTTATTGAAGGTATCg                                                                                                                                                                      >  1:4436243/1‑251 (MQ=255)

                                                                                gTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTTTCAGACGCTTCTCAAAATTAGCATTTAATTGCCTTTTTTTTTCCCGCGTCAGGGTGGTTTAAATTTTTTTCTGGGGTTCCATGTGCCCCAGGCTGTTTGCTTTGCGCTGTTCTTCCCTTTTTGCGGTCTGCGCCTCTGCTTTTTCTTTTTCGAAGGttttg                                                                                                                                                                      >  2:5459085/1‑247 (MQ=255)

                                                                                    tGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCt                                                                                                                                                                    <  1:3868625/249‑1 (MQ=255)

                                                                                     ggACAAGACCGGCGGTCTTTAGTTTTTTGGCTGGAAGATTTTTCAGCAGTTGCAGGCCTTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTTAATTGCTTTTTCTGTTTCCGCGTCAGGATGTTGTAAATTTTTTTTTTGGGTTTCATTTGCGCCAGCATGTTTGCTTTGCGCTGTTCTTCCCTTTTTGCGATCTGCGCTTCAGCTTTTTCTTTTTCGAAGGTATCg                                                                                                                                                                      >  2:2432345/1‑246 (MQ=255)

                                                                                     ggACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGGTGTTGTAAATTTTGTTCTGGGTTTCCCTGTGCGCCAGCATGTTAGCTTTGGGCTGGTTTTCCATTTTTTGGATCTTCGCTTCCGCTTTTACttttttgtaaggatcg                                                                                                                                                                      >  1:4461251/1‑233 (MQ=255)

                                                                                     ggACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCg                                                                                                                                                                               <  1:5901611/237‑1 (MQ=255)

                                                                                     ggACAAGACCGGCGGCCTTAAGCTTTTTGGCTGAAAGATTATTCATCCGTTGCAGGCCTTTTTCCTTTTGCCGCTGGACGCTTTGTCAGACGCCTCTCAAAATTAGCATTTAATTTCTTTTTCTGTTCCGGCGGCAGGGTGTTGTAAATTTTGTTTTGGGGTTTCATGTGCGCCAGCATGTTatg                                                                                                                                                                                                                                   >  2:1925308/1‑183 (MQ=255)

                                                                                        cAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCAtg                                                                                                                                                                                                                                                    <  2:1212527/165‑1 (MQ=255)

                                                                                           gACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTTTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATTTGCGCCCGCATGTTAGCTTTGCGCTGTTTTTCCATTTTTTCGATCTGCGCTTCAGCTTTTTCTTTATCGAAGGTATCGCTTGCAatgat                                                                                                                                                           >  1:3964349/1‑251 (MQ=255)

                                                                                           gACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGc                                                                                                                                                                                                                                   <  2:4473749/179‑1 (MQ=255)

                                                                                             ccGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGt                                                                                                                                                         <  2:1302264/251‑1 (MQ=255)

                                                                                             ccGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGt                                                                                                                                                         <  2:2839781/251‑1 (MQ=255)

                                                                                             ccGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGt                                                                                                                                                         <  2:836362/251‑1 (MQ=255)

                                                                                                   gTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCa                                                                                                                                                                                                                                          <  2:4519011/164‑1 (MQ=255)

                                                                                                   gTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCa                                                                                                                                                                                                                                          >  1:4519780/1‑164 (MQ=255)

                                                                                                    tCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAGTTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCAt                                                                                                                                                  <  2:1447098/251‑1 (MQ=255)

                                                                                                    tCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCAt                                                                                                                                                  <  1:4137994/251‑1 (MQ=255)

                                                                                                    tCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCAt                                                                                                                                                  <  2:2020196/251‑1 (MQ=255)

                                                                                                    tCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCAt                                                                                                                                                  <  2:1536516/251‑1 (MQ=255)

                                                                                                           ttttttGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCgcgg                                                                                                                                           >  1:4379975/1‑251 (MQ=255)

                                                                                                           ttttttGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCACGCATTgcgc                                                                                                                                             <  1:1993275/249‑1 (MQ=255)

                                                                                                                   cTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCa                                                                                                                                                                                                                  >  2:2401342/1‑172 (MQ=255)

                                                                                                                   cTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGttcttc                                                                                                                                   <  1:5454579/251‑1 (MQ=255)

                                                                                                                   cTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGttcttc                                                                                                                                   <  2:4435475/251‑1 (MQ=255)

                                                                                                                       aaGTTTATTCAGCAGTTGCAGGCATTTTACTTTTTGCCGCTGGTCGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATNNNNNNNNNNNNNNNCNATGATGTCATGCATTGCGCGGCGTTCTTCCAgc                                                                                                                               <  1:5273780/251‑1 (MQ=255)

                                                                                                                        aGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttcttc                                                                                                                                                                                                                    <  1:1552965/165‑1 (MQ=255)

                                                                                                                        aGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATGAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGc                                                                                                                                                                                                                                   <  1:2946736/150‑1 (MQ=255)

                                                                                                                             attCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACtt                                                                                                                                                                                    <  1:4695821/192‑1 (MQ=255)

                                                                                                                             attCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGACGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGcggcgg                                                                                                                          >  2:2636911/1‑250 (MQ=255)

                                                                                                                                     gTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAg                                                                                                                                                                                                                                    >  2:5363718/1‑136 (MQ=255)

                                                                                                                                     gTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGc                                                                                                                                                                                                                                   <  2:5615119/137‑1 (MQ=255)

                                                                                                                                     gTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGc                                                                                                                                                    >  2:4767624/1‑216 (MQ=255)

                                                                                                                                     gTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGc                                                                                                                                                    <  1:4768448/216‑1 (MQ=255)

                                                                                                                                      ttGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATc                                                                                                                <  1:3479030/251‑1 (MQ=255)

                                                                                                                                        gCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTctt                                                                                                                 >  2:1778989/1‑246 (MQ=255)

                                                                                                                                            gCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCtt                                                                                                                                                                                                                                 >  2:5883882/1‑132 (MQ=255)

                                                                                                                                             cATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCa                                                                                                                                                                                                                                                      >  1:763524/1‑110 (MQ=255)

                                                                                                                                             cATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttcttc                                                                                                                                                                                                                    >  1:1008762/1‑144 (MQ=255)

                                                                                                                                             cATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttcttc                                                                                                                                                                                                                    <  2:1008554/144‑1 (MQ=255)

                                                                                                                                                   aCCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCt                                                                                                               <  1:4240761/239‑1 (MQ=255)

                                                                                                                                                    ccTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAatgatg                                                                                                                                                          >  2:4660657/1‑195 (MQ=255)

                                                                                                                                                         tGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCa                                                                                             <  2:2185922/251‑1 (MQ=255)

                                                                                                                                                           ccGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCAt                                                                                                                                                  >  2:1618134/1‑196 (MQ=255)

                                                                                                                                                              cTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCa                                                                                                                                                                                                                  <  2:4110890/129‑1 (MQ=255)

                                                                                                                                                                ggACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgcgc                                                                                                                                                                                                                            >  1:1847022/1‑117 (MQ=255)

                                                                                                                                                                       cTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATc                                                                                                                <  1:3184801/218‑1 (MQ=255)

                                                                                                                                                                       cTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTTGCGGAtt                                                                               >  2:4013786/1‑250 (MQ=255)

                                                                                                                                                                         gTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGt                                                                                                                                                         <  1:5834257/175‑1 (MQ=255)

                                                                                                                                                                                              gCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGc                                                                                                       <  2:3207916/204‑1 (MQ=255)

                                                                                                                                                                                                 ttAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGtcaggtt                                                     <  1:288019/251‑1 (MQ=255)

                                                                                                                                                                                                   aaaTTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCa                                                                                                                                                                                                                                                      <  1:2567509/56‑1 (MQ=255)

                                                                                                                                                                                                       tGCTTTTTCTTTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTc                                               <  1:5789034/251‑1 (MQ=255)

                                                                                                                                                                                                            tttCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGa                                          <  2:4379211/251‑1 (MQ=255)

                                                                                                                                                                                                                 gTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgcgc                                                                                                                                                                                                                            >  2:5746904/1‑68 (MQ=255)

                                                                                                                                                                                                                     cGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGt                                 <  1:891601/251‑1 (MQ=255)

                                                                                                                                                                                                                        cGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGt                                                                                                                                                         >  1:4831890/1‑128 (MQ=255)

                                                                                                                                                                                                                        cGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGt                                                                                                                                                         <  2:4831055/128‑1 (MQ=255)

                                                                                                                                                                                                                         gTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCa                                                                                                                                                                                                                  <  1:5161782/70‑1 (MQ=255)

                                                                                                                                                                                                                         gTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCa                                                                                                                                                                                                                  >  2:5160909/1‑70 (MQ=255)

                                                                                                                                                                                                                          tCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCaa                                                                                                                                                               <  2:3386230/120‑1 (MQ=255)

                                                                                                                                                                                                                          tCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCaa                                                                                                                                                               >  1:3386749/1‑120 (MQ=255)

                                                                                                                                                                                                                              gATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTcc                               <  1:1452835/244‑1 (MQ=255)

                                                                                                                                                                                                                              gATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTcc                               <  2:4927111/244‑1 (MQ=255)

                                                                                                                                                                                                                              gATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATg                          <  2:3088578/249‑1 (MQ=255)

                                                                                                                                                                                                                              gATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATg                          <  1:11969/249‑1 (MQ=255)

                                                                                                                                                                                                                                 gttgtAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTgcg                                                                  <  1:5393007/206‑1 (MQ=255)

                                                                                                                                                                                                                                        aTTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCAt                                                                                                                                                  >  1:407201/1‑119 (MQ=255)

                                                                                                                                                                                                                                        aTTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGACCGAACtt              >  1:1736758/1‑251 (MQ=255)

                                                                                                                                                                                                                                           ttGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGACCGAACTTGcc           <  2:5675566/251‑1 (MQ=255)

                                                                                                                                                                                                                                           ttGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGACCGAACTTGcc           <  1:2760798/251‑1 (MQ=255)

                                                                                                                                                                                                                                           ttGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGACCGAACTTGcc           <  1:1244013/251‑1 (MQ=255)

                                                                                                                                                                                                                                                  gggTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGGCCGAACTTGc            >  2:4231445/1‑243 (MQ=255)

                                                                                                                                                                                                                                                    gTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGc                                                                                                                                                    <  2:1116287/105‑1 (MQ=255)

TAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCCATTTTTGCGATCTGCGCTTCAGCTTTTACTTTATCGAAGGTATCGCTGGCAATGATGTCATGCATTGCGCGGCGTTCTTCCAGCGGCGGACGTTTCATCTGGTCACGCTGGCCTTTCATGATTTCGCGGATCTGCTGTTTCTGCGCGTCGGTCAGGTTCAGGTCTTTGAACATCATGTCCTGATGCGGACCGAACTTGCCTTTGTGGTG  >  NC_000913/1825017‑1825511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 8 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 39 ≤` `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.