Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	2,823,830	0	A→G	24.7%	50.6 ‑0.2 ‑0.1	85	intergenic (‑61/+19)	recA/pncC	DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor/nicotinamide‑nucleotide amidohydrolase; NMN amidohydrolase
Reads supporting (aligned to +/- strand):  new base (12/9):  ref base (33/31):  total (45/40)
Fisher's exact test strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.84e-01

AACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCG  >  NC_000913/2823637‑2824024                                                                                                                                                                                                  |                                                                                                                                                                                                   aaCCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTAt                                                                                                                                                                                              <  2:160321/200‑1 (MQ=255)     ggTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATAc                                                                                                                                                                                                   <  2:35889/191‑1 (MQ=255)         gagaTGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATAc                                                                                                                                                                                      <  2:553262/200‑1 (MQ=255)         gagaTGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                                                      <  2:65799/200‑1 (MQ=255)                 ttCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCaa                                                                                                                                                                              >  2:249217/1‑200 (MQ=255)                                   gTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                                                            >  2:349255/1‑200 (MQ=255)                                   gTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTg                                                                                                                                                                                                 >  2:77491/1‑163 (MQ=255)                                       tCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGcc                                                                                                                                                        >  2:684706/1‑200 (MQ=255)                                                ggcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTGTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAAGACCGGTTAGTCAATATGTTCTGTTGAAGCAATTATACTg                                                                                                                                                                                                 <  3:21506/149‑1 (MQ=255)                                                cgcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTg                                                                                                                                                                                                 <  3:130417/150‑1 (MQ=255)                                                cgcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTg                                                                                                                                                                                                 <  3:420607/150‑1 (MQ=255)                                                cgcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTg                                                                                                                                                                                                 <  3:77491/150‑1 (MQ=255)                                                  cATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTa                                                                                                                                                                                               <  3:630210/150‑1 (MQ=255)                                                         gAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTg                                                                                                                                                                                                 <  4:404832/141‑1 (MQ=255)                                                           gCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATAc                                                                                                                                                                                      <  3:173620/150‑1 (MQ=255)                                                           gCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATCTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTa                                                                                                                                                                                               >  4:61925/1‑141 (MQ=255)                                                           gCCTGTACCAAATTGTTACTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCAGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTa                                                                                                                                                                                               <  4:21753/141‑1 (MQ=255)                                                                tACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATAc                                                                                                                                                                                                   <  1:96433/132‑1 (MQ=255)                                                                tACCAAATTGTTTCTCAATCGGGCCCAGTGCTGCCGCCAACGCTGTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTGATGCCGGGTAATACCGGATAGTCAATATGTACTGTTGAAGCAATTATAc                                                                                                                                                                                                   <  1:69050/132‑1 (MQ=255)                                                                 aCCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTa                                                                                                                              >  2:354397/1‑200 (MQ=255)                                                                     aaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgtg                                                                                                                                                                           <  2:138128/151‑1 (MQ=255)                                                                     aaTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgt                                                                                                                                                                            >  3:138128/1‑150 (MQ=255)                                                                      aTTGTTGCTCTATCTGGCCCAGTGCTGCCGCCAACTCTTTCTGTTTGTTTTCGTCGATATCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgtn                                                                                                                                                                           <  3:649521/150‑2 (MQ=255)                                                                               cAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAg                                                                                                                                                                                    <  1:297158/132‑1 (MQ=255)                                                                               cAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAgtg                                                                                                                                                                           >  4:66384/1‑141 (MQ=255)                                                                                 aTCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTa                                                                                                                                                                                  >  1:232921/1‑132 (MQ=255)                                                                                 aTCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGaaa                                                                                                                                                                >  3:284480/1‑150 (MQ=255)                                                                                     ggCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATc                                                                                                          >  2:237077/1‑200 (MQ=255)                                                                                      gaccAGTGCTGCCGCCAACGCTTTCTGTTTGTTTACGTCGATAGCCATTTTTACTCATGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTGGAAGCAATTATGCTGTATGCTCATACAGTATCAAg                                                                                                                                                                             <  1:337867/130‑1 (MQ=255)                                                                                        ccAGTGCTGCCGCCAACGGTTTCTGTTTGTTTTCGTCGATAGCCATTTTTGCTCCTGTCATGCCGGGGAATACCGGATAGTCAATATGTTCTGTTGGGGCAATTATAGTGTATGCTGATACAGAATCAAgtg                                                                                                                                                                           <  1:187012/132‑1 (MQ=255)                                                                                          aGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGt                                                                                                     <  2:646830/200‑1 (MQ=255)                                                                                          aGTGCTGCCGCCAACGCTTTCTGTTTGTTTACGTCGATAGCCATTTTTACTCCTGTCATGCCGCGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGt                                                                                                     <  2:253342/200‑1 (MQ=255)                                                                                           gTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCcac                                                                                                                                                      >  3:221657/1‑150 (MQ=255)                                                                                           gTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTc                                                                                                    >  2:527951/1‑200 (MQ=255)                                                                                                      aaCGCGTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCaa                                                                                                                                           <  3:487453/150‑1 (MQ=255)                                                                                                        cGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAgg                                                                                                                                                  <  4:118407/141‑1 (MQ=255)                                                                                                         gttttCTGGTTGTTTTCGTCGATATCCATTTTTACTCCTGTCAAGCCGGGTAATACCGGATAGGCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTGTTGTAGAAATTGTTGCCACAAGGTCTGCAATGc                                                                                                                                        <  3:587878/148‑1 (MQ=255)                                                                                                            ttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAg                                                                                                                                                                                    >  2:752518/1‑103 (MQ=255)                                                                                                            ttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAg                                                                                                                                                                                    <  3:752471/103‑1 (MQ=255)                                                                                                                 tttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAg                                                                                                                                <  3:704427/150‑1 (MQ=255)                                                                                                                      tttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGcc                                                                                                                           <  3:632920/150‑1 (MQ=255)                                                                                                                      tttCGTCGATAGCCATTTTTACGCCTGTCCTGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGcc                                                                                                                           <  3:758294/150‑1 (MQ=255)                                                                                                                        tcgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGc                                                                                                                                  <  4:386107/141‑1 (MQ=255)                                                                                                                         cgtcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTga                                                                                                                        <  3:349234/150‑1 (MQ=255)                                                                                                                           tcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGc                                                                                                                                        <  1:244859/132‑1 (MQ=255)                                                                                                                                   cATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAAGTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAg                                                                                                                                <  1:416166/132‑1 (MQ=255)                                                                                                                                     tttttACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgacga                                                                                                                     <  4:107059/141‑1 (MQ=255)                                                                                                                                     tttttAATCCTGTCATGCCGGGTAATACCGGAACGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGCACTTGTTGCCACAAGGTCTGCAATGCATACGCAGTATCCTGACGACGTACCgac                                                                                                            >  3:627831/1‑148 (MQ=255)                                                                                                                                      ttttACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTgacgac                                                                                                                    >  4:375145/1‑141 (MQ=255)                                                                                                                                      ttttACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCa                                                         >  2:75829/1‑200 (MQ=255)                                                                                                                                           cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGc                                                                                                                                        <  4:6144/116‑1 (MQ=255)                                                                                                                                           cTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGc                                                                                                                                        >  1:6144/1‑116 (MQ=255)                                                                                                                                              cTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCaaa                                                 >  2:79217/1‑200 (MQ=255)                                                                                                                                                gTCATCCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTcgcc                                                                                                 >  3:268041/1‑150 (MQ=255)                                                                                                                                                      ccGGGTAATACCGGATAGTCAATATGTTCGGTTGAAGCAATTATACGGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGc                                                                                                             <  1:351512/132‑1 (MQ=255)                                                                                                                                                        gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCCCCGCATCACGGTCGCCGCTGAAGc                                                                                         >  3:24913/1‑150 (MQ=255)                                                                                                                                                        gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGc                                                                                         >  3:750634/1‑150 (MQ=255)                                                                                                                                                        gggTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAg                                       >  2:170296/1‑200 (MQ=255)                                                                                                                                                             aTACCGGATAGTCAATATGTTCTGTTGAAGCACTTATACTGTATGCTCATACAGTATCAAGTGTTCTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCTCgg                                                                                                      >  1:227193/1‑132 (MQ=255)                                                                                                                                                                ccGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGcc                               >  2:702194/1‑200 (MQ=255)                                                                                                                                                                       gtccctcttttCTTTTTAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTccc                                                                                   >  4:71688/10‑141 (MQ=255)                                                                                                                                                                       gTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACTCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTaa                                                                          >  3:760539/1‑150 (MQ=255)                                                                                                                                                                       gTCAATATGTTCTGTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTccc                                                                                   >  4:343566/1‑141 (MQ=255)                                                                                                                                                                             atGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCaa                                                                                                                                           <  3:662520/79‑1 (MQ=255)                                                                                                                                                                             atGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCaa                                                                                                                                           >  2:662558/1‑79 (MQ=255)                                                                                                                                                                             atGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCAtt                                                                                      >  1:34358/1‑132 (MQ=255)                                                                                                                                                                               gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTct                                                                                    >  1:58696/1‑131 (MQ=255)                                                                                                                                                                               gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTa                                                                           >  4:327321/1‑141 (MQ=255)                                                                                                                                                                               gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCa                                                                  >  3:339106/1‑150 (MQ=255)                                                                                                                                                                               gTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTc                <  2:544667/200‑1 (MQ=255)                                                                                                                                                                                ttCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGTGTAATGCCTTCAc                                                                 <  3:6812/150‑1 (MQ=255)                                                                                                                                                                                  cTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATg                                                                        >  4:9095/1‑141 (MQ=255)                                                                                                                                                                                  cTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACt             <  2:235379/200‑1 (MQ=255)                                                                                                                                                                                   tGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACcgc                                                              >  3:573390/1‑150 (MQ=255)                                                                                                                                                                                    gTTGAAGCAATTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGcc                                                                      >  4:56589/1‑141 (MQ=255)                                                                                                                                                                                           cAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCTACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGCAATg                                                                        >  1:43156/1‑132 (MQ=255)                                                                                                                                                                                           cAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTg                                                      >  3:96512/1‑150 (MQ=255)                                                                                                                                                                                           cAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATc    >  2:682300/1‑200 (MQ=255)                                                                                                                                                                                            aaTTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGc                                                                       <  1:332649/132‑1 (MQ=255)                                                                                                                                                                                             aTTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGcc                                                                      >  1:48059/1‑132 (MQ=255)                                                                                                                                                                                             aTTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGcc                                                                      >  1:96875/1‑132 (MQ=255)                                                                                                                                                                                             aTTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCCGTGGc                                                    >  3:527837/1‑150 (MQ=255)                                                                                                                                                                                             aTTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGc                                                    >  3:203701/1‑150 (MQ=255)                                                                                                                                                                                             aTTATGCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATccg  >  2:352560/1‑200 (MQ=255)                                                                                                                                                                                             aTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACcgcg                                                             >  4:74113/1‑141 (MQ=255)                                                                                                                                                                                                ataCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAgg                                                                                                                                                  >  1:95475/1‑53 (MQ=255)                                                                                                                                                                                                ataCTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAgg                                                                                                                                                  <  4:95477/53‑1 (MQ=255)                                                                                                                                                                                                  |                                                                                                                                                                                                   AACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCG  >  NC_000913/2823637‑2824024

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