Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,781,152	0	A→G	14.9%	10.4 ‑28.0 ‑0.0	97	intergenic (+135/‑63)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/14):  ref base (75/5):  total (77/20)
Fisher's exact test strand distribution p-value = 6.67e-13
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.85e-01

GGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCG  >  NC_000913/3780963‑3781342                                                                                                                                                                                                   |                                                                                                                                                                                                gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAgg                                                                                                                                                                                      >  4:321843/1‑200 (MQ=255) gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAgg                                                                                                                                                                                      >  4:12759/1‑200 (MQ=255)  ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                                                                                                                     <  4:208229/200‑1 (MQ=255)  ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                                                                                                                     >  4:781458/1‑200 (MQ=255)       agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑ACAAGGCCGCAc                                                                                                                                                                                >  4:263547/1‑200 (MQ=255)        gagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑ACAAGGCCGCACt                                                                                                                                                                               >  4:501944/1‑200 (MQ=255)                                           gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAgg                                                                                                                                                                                               >  2:664913/1‑150 (MQ=255)                                             aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAgggg                                                                                                                                                                                             >  2:864852/1‑150 (MQ=255)                                              ggcgccatgtcctttgccccccccCCACCCTACCCCTC‑CCCCCCGGGGGGGGGGAGCGGTCGGGGCTTTTTTTTACCCCCCCCTCGGAAACGACCACCCTCTTTTTTTTTCCTCCCCCCCCCCGC‑G‑GGGGGGG‑GGGCGGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                                                                                                                         <  4:95292/181‑1 (MQ=255)                                                  gcgtAGTCGTTTGCGCCGCCTCACCCTACCCCTCTCCCTCAGGAAGAGGGGACCCTTCGGCGCTGTATGTACTACCTCACTCTGAACCGAGACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGCGA‑GGGACGGGGT‑GAgcgacaaag                                                                                                                                                                                        >  2:559819/4‑142 (MQ=255)                                                  gCATAGTAGATTGAGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCAGCATGTACACCCTCACTCTGAAACGACAAAGCACTCTTTTTGT‑CTCCCTCGCCCTTC‑C‑GGGGAGA‑GGTCCGGGGT‑GAGCGG‑AAAAg                                                                                                                                                                                       >  2:168663/1‑150 (MQ=255)                                                      aGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAgggg                                                                                                                                                                                             <  3:231747/141‑1 (MQ=255)                                                             tGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                                                                                                                     <  3:227522/141‑1 (MQ=255)                                                                   cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTCTGGACTCCCCCACCCTGACACGACACCGCCCTCTTTTCTT‑CTCCCTCGCCCCTC‑C‑GGCGACA‑GGTCCGGGGT‑CCGGGG‑aaca                                                                                                                                                                                        >  1:440883/1‑130 (MQ=25)                                                                   cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCTCTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                                                                                                                    >  4:663570/1‑200 (MQ=255)                                                                   cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑aaac                                                                                                                                                                                        >  1:799515/1‑131 (MQ=255)                                                                   cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACt                                                                                                                                                                               >  3:790413/1‑141 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑TAGGGG‑AAAAGGCCGCACTGCTccc                                                                                                                                                                         >  3:626240/1‑139 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGc                                                                                                                                                                                  >  1:456010/1‑132 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCAc                                                                                                                                                                         >  3:772434/1‑141 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCAc                                                                                                                                                                         >  3:641605/1‑141 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                                                                                >  2:738886/1‑150 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCTGCACTAAGGGGGCAGTCTGCTAAACATCGTACTTTa                                                                                                              >  4:21835/1‑200 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                                                                                                              >  4:531344/1‑200 (MQ=255)                                                                         cccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGG‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                                                                                                              >  4:257037/1‑200 (MQ=255)                                                                          ccTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTCCTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GTGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                                                                                                             <  4:257377/200‑1 (MQ=255)                                                                               cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GATGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                                                                                                                          >  2:1007775/1‑150 (MQ=255)                                                                                 ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                                                                                                                        >  2:235517/1‑150 (MQ=255)                                                                                     ccctgggggAGGGGGGCCCTTCGGGCTCTGTTTGTCCCCCCCCCCTCGAAAACGCCCCCACCCT‑TTTTTTT‑CCCCCCCCCCCCCC‑C‑GGGGGGG‑GGGGCGGGGG‑GGGGGGAAAAAGGCCGCACTGCTCACCaa                                                                                                                                                                      <  1:481203/126‑1 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                                                                                             >  3:442574/1‑141 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                                                                                             >  3:258161/1‑141 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                                                                                                             >  3:747149/1‑141 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                                                                    >  2:282038/1‑150 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                                                                    >  2:360650/1‑150 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                                                                                                                    >  2:513780/1‑150 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt                                                                                                  >  4:253707/1‑200 (MQ=255)                                                                                     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt                                                                                                  >  4:560511/1‑200 (MQ=255)                                                                                      ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAtt                                                                                                 >  4:727847/1‑200 (MQ=255)                                                                                       cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCaa                                                                                                                                                                      >  3:549444/1‑130 (MQ=255)                                                                                         cAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                                                                                                                         >  3:791369/1‑141 (MQ=255)                                                                                           gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                                                                                                              >  2:596419/1‑150 (MQ=255)                                                                                           gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                                                                                                              >  2:1048742/1‑150 (MQ=255)                                                                                              agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaaa                                                                                         >  4:735161/1‑200 (MQ=255)                                                                                                 ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGCCCCCACTCCCCACCATTTTTCTCCTTTc                                                                                                                                                          >  1:519903/1‑132 (MQ=38)                                                                                                        gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAc                                                                               >  4:901094/1‑200 (MQ=255)                                                                                                             gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                                                                                                                     >  3:728744/1‑141 (MQ=255)                                                                                                             gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                                                                                                            >  2:1076007/1‑150 (MQ=255)                                                                                                                cTGTTTGTACCCCCTCCCTCGAAACCGACCCCGCCCTTTTTTTTC‑CCCCCCCCCCCCTC‑G‑GGGGGGG‑GGGGGGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                                                                                                                         <  2:387934/150‑1 (MQ=255)                                                                                                                cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAAt                                                                                        >  4:135941/1‑183 (MQ=255)                                                                                                                  gtattaccccccTAACTCAAAACCGACACCGC‑CTCTTTTTTT‑C‑CCCCCGCGCCCC‑G‑GGGGGGAGGGGGCGGGGG‑GAGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     <  4:368571/191‑1 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCTCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                                                                                                                       >  2:663764/1‑150 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCCCTACgg                                                                                                                                         >  1:186006/1‑132 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                                                                                                                         >  1:617380/1‑132 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                                                                                                >  3:298761/1‑141 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                                                                                                >  3:400786/1‑141 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAACCAt                                                                                                                       >  2:157748/1‑150 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                                                                                                                       >  2:787845/1‑150 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGACATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     >  4:187611/1‑200 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     >  4:734396/1‑200 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     >  4:625470/1‑200 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGG‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCATAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     >  4:419654/1‑200 (MQ=255)                                                                                                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CACCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     >  4:549652/1‑200 (MQ=255)                                                                                                                    atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                                                                                                                     >  3:656030/1‑134 (MQ=255)                                                                                                                    atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                                                                                                                     >  2:1029334/1‑150 (MQ=255)                                                                                                                    atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                                                                                                                     >  2:766232/1‑150 (MQ=255)                                                                                                                    atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCtt                                                                   >  4:819100/1‑200 (MQ=255)                                                                                                                    atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGC‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                                                                                                                     >  2:258334/1‑150 (MQ=255)                                                                                                                      gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTg                                                                 >  4:152749/1‑200 (MQ=255)                                                                                                                         cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCACTTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                                                                                                                  >  1:618839/1‑132 (MQ=255)                                                                                                                                tcttaaaacccacacCACATTTTTTTTCC‑CCCCCCCCCCCCTC‑G‑GGGGAGC‑GGGCGGGGTG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                                                                                                                  <  3:66003/130‑1 (MQ=255)                                                                                                                                 ctctGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc                                                                                                        >  2:657771/1‑150 (MQ=255)                                                                                                                                 ctctGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCg                                                                       >  4:934476/1‑183 (MQ=255)                                                                                                                                 ctctGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCg                                                      >  4:1042046/1‑200 (MQ=255)                                                                                                                                   ctGAAACGCCACCTCACTCTTTTTTT‑CTCCCTCGCCCTTC‑C‑GGGGACA‑CGGCCGGGGT‑GTGGGG‑AAAACGCCTCACTGCTCACCAATTCTCCGCTCTCCTCCCCCCGCACTTACGGGGCACTATGCCACACa                                                                                                                        >  1:49723/1‑132 (MQ=255)                                                                                                                                       aaCGCCCCCGCCCTTTTTTTTC‑CCCCCCCCCCCCTC‑G‑GGGGGGG‑GGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt                                                                                                  <  2:185812/150‑1 (MQ=255)                                                                                                                                         gcacacAGCATTTTTTTTCC‑CCCCCCCCCCCCCC‑G‑GGGGAGG‑GGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                                                                                                                  <  1:442555/130‑1 (MQ=255)                                                                                                                                                ccctctttttttt‑CCCCCCCCCCCCCCGC‑GGGGAGG‑GGCCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCat                                        <  4:910161/194‑1 (MQ=255)                                                                                                                                                  ccttttttttt‑CCCCCCCCCCCCCC‑CGGGGGGGG‑GGGCGGGGGT‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAAt                                                                                        <  2:135938/146‑1 (MQ=255)                                                                                                                                                     ctTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcatc                                  >  4:925303/1‑200 (MQ=255)                                                                                                                                                     ctTTTTTT‑CCCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGTCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACACCGTACTTTACTAACCAg                                                                                                      >  1:688331/1‑132 (MQ=255)                                                                                                                                                       tttttt‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaa                                                                                                    >  1:566574/1‑132 (MQ=255)                                                                                                                                                        ttttt‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaa                                                                                          >  3:155751/1‑141 (MQ=255)                                                                                                                                                               tcCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAACTACTg                                                                                    >  3:178909/1‑141 (MQ=255)                                                                                                                                                                   cccccccCC‑G‑GGGGGGG‑GGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCg                                                                       <  2:934404/147‑1 (MQ=255)                                                                                                                                                                      ccccGC‑G‑GGGGGGG‑GGGCGGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATAc                                                                                      <  1:355233/132‑1 (MQ=255)                                                                                                                                                                        ccTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCttt                                                                  >  2:732861/1‑150 (MQ=255)                                                                                                                                                                          tC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCACATACTGGc                                                                                  >  1:587080/1‑132 (MQ=255)                                                                                                                                                                           c‑c‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTgg             >  4:38605/1‑200 (MQ=255)                                                                                                                                                                           c‑c‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTgg             >  4:927990/1‑200 (MQ=255)                                                                                                                                                                                ggggggggggggggggG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCt                                                             <  2:970880/138‑1 (MQ=255)                                                                                                                                                                                 gggggggggCGGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAca                                                                              <  1:747113/126‑1 (MQ=255)                                                                                                                                                                                 ggggggggcggggggG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                     <  3:352152/130‑1 (MQ=255)                                                                                                                                                                                  gaga‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCtt                                                                   >  3:29409/1‑141 (MQ=255)                                                                                                                                                                                   gggggggGCGGGGTGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCt                                                                    <  3:48083/137‑1 (MQ=255)                                                                                                                                                                                    ga‑gGGCGGGGGT‑GAGGAG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATAc                                                        <  2:253698/150‑1 (MQ=255)                                                                                                                                                                                       gggCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGTCGATAAg    >  4:426711/1‑200 (MQ=255)                                                                                                                                                                                         gCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAgcg  >  4:263154/1‑200 (MQ=255)                                                                                                                                                                                                   aGGGG‑AAAAGGCCGCACTGCTCACAAAGTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGc                                                    <  3:587101/141‑1 (MQ=255)                                                                                                                                                                                                   |                                                                                                                                                                                                GGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTT‑CTCCCTCGCCCCTC‑C‑GGGGAGA‑GGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCG  >  NC_000913/3780963‑3781342

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ 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.