Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,875,529	0	C→G	23.5%	22.1 ‑9.6 0.0	68	S31C (TCC→TGC)	yidX	putative lipoprotein
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 (16/0):  ref base (24/28):  total (40/28)
Fisher's exact test strand distribution p-value = 6.51e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACC  >  NC_000913/3875333‑3875711                                                                                                                                                                                                     |                                                                                                                                                                                       ataaCATACGTAGGGAAGAAATCAGAAAAGGGATCGGAGGCGCTACACAGAAGGAAAAGCAGTGTATGCTTTCGACGGCATAGCTACTCGTTTAAAGACGTCATCGTGAAGTTGTAATTGAAGGGATTTTTAAAGGCTGCTGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCa                                                                                                                                                                                     <  4:116509/197‑1 (MQ=255)         agacGGGCAGAAATCTGAAAAGGGATATCAGACGCTGAACAGGAGGGAAAAAAGAGTATGCATTCAACAGCTTAGCTACGCGTTTAAAAGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCaa                                                                                                                                                                             <  4:649495/197‑1 (MQ=255)                                        gcgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAGCAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTc                                                                                                                                              >  4:715229/1‑200 (MQ=255)                                        gcgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTATACAGCAAAGGATAGTTCAGGAAAGTATCAATCGCAGTc                                                                                                                                              >  4:817105/1‑200 (MQ=255)                                                cAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCTATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCa                                                                                                                                      >  4:33505/1‑200 (MQ=255)                                                    aggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCTTTTAAAGGATTAATCATGAATTTGAATTTTAAGGGATTTTTTAAGTCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAt                                                                                                                                                                                    <  3:700113/150‑1 (MQ=255)                                                     gggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCTTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACt                                                                                                                                 >  4:234413/1‑200 (MQ=255)                                                     gggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCTTTTGCTTGAAAAAAAAAGATAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCaac                                                                                                                                 >  4:402084/1‑198 (MQ=255)                                                      ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCTTCTGCTTGAAAAAAAAAGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACCATCACTg                                                                                                                                >  4:643949/1‑200 (MQ=255)                                                           gaagaGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTGCTTGACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTa                                                                                                                           >  4:569941/1‑200 (MQ=255)                                                              gagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCg                                                                                                                        >  4:26462/1‑200 (MQ=255)                                                              gagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTATACAAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCg                                                                                                                        >  4:625174/1‑200 (MQ=255)                                                               aatATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGAGTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGc                                                                                                                       <  4:454827/198‑1 (MQ=255)                                                                    gCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTccc                                                                                                                                                                                      >  2:415133/1‑132 (MQ=255)                                                                      tttCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTGTAAGGCTGCCGGTTTAGTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAt                                                                                                                                                                                    <  2:594205/132‑1 (MQ=255)                                                                       ttCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAg                                                                                                                                                                 >  3:412611/1‑150 (MQ=255)                                                                              gCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTGCCATAACACAAAGAGTAGTTc                                                                                                                                                                   >  1:578115/1‑141 (MQ=255)                                                                               gttAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCGATCGCAGGCAAAAGATAGTaa                                                                                                       <  4:436173/199‑1 (MQ=255)                                                                                  aGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATAaa                                                                                                                                                      >  3:501414/1‑150 (MQ=255)                                                                                  aGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCTTTAAAGAAAAGAGTAGTTCAGGa                                                                                                                                                               >  1:580797/1‑141 (MQ=255)                                                                                   gCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                              <  1:43428/141‑1 (MQ=255)                                                                                   acTACTCGTTTGAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTAAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCGGGCTGTATCTCGTATGCGCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                              <  1:372286/140‑1 (MQ=255)                                                                                       cTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTc                                                                                                                                                                   <  2:38124/132‑1 (MQ=255)                                                                                          gTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACAAAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAg                                                                                            >  4:216154/1‑200 (MQ=255)                                                                                          gTTTAAAGGATTAATCATGAAGTTAAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCTAGACAGCAAAAGTTAGTTCAGGAAAGTATCAATTGCAGTTAGACACCATCACTGGTCTATCGCAGGAAAAAGATAGTAATGGAACAATAg                                                                                            >  4:789587/1‑200 (MQ=255)                                                                                            ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGaa                                                                                                                                                              <  2:53926/132‑1 (MQ=255)                                                                                            ttAAAAGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                                                          <  4:762841/200‑1 (MQ=255)                                                                                             tAAAGGATAAACAATGGAGTTAGAATTTAAGGAATTTATTAAGGCTGCTAGTCGATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAGCAATc                                                                                                                                                    <  1:246768/141‑1 (MQ=255)                                                                                                 ggATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACAAAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                     >  4:358672/1‑200 (MQ=255)                                                                                                  gATAAATCATGAAGTTGAATTTAAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATATCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                                                                    <  4:188907/200‑1 (MQ=255)                                                                                                       aTCATGAAGATGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGCAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACt                                                                                                                                 >  3:464518/1‑150 (MQ=255)                                                                                                            gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTAGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTc                                                                                                                                              <  2:327841/132‑1 (MQ=255)                                                                                                               gTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAt                                                                                                                                                         >  2:74193/1‑118 (MQ=255)                                                                                                               gTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCAc                                                                                                                                  >  1:350046/1‑141 (MQ=255)                                                                                                                     ttttAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAACTACATCAGGAAAGTATCCACCGCAGTCAGCCACCATTACAGGACTATCGAAGGCAAAAGATAGAAATGGAACAAAAGGATATGTTTATGTAGGGGAATCGCTgg                                                                 >  4:844562/1‑200 (MQ=255)                                                                                                                        tAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCAc                                                                                                                                  <  2:184329/132‑1 (MQ=255)                                                                                                                            ggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                                                          >  4:51818/1‑200 (MQ=255)                                                                                                                            ggATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTACAACAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                                                          >  4:841445/1‑200 (MQ=255)                                                                                                                             gATTTTTGAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGc                                                                                                                             <  2:114094/132‑1 (MQ=255)                                                                                                                              aTTTTATAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCATATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGc                                                                                                                   <  1:537203/141‑1 (MQ=255)                                                                                                                                 ttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATc                                                                                                                                    >  4:118703/1‑121 (MQ=255)                                                                                                                                 ttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATc                                                                                                                                    <  3:118701/121‑1 (MQ=255)                                                                                                                                              ggTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGACAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                                                          >  3:158123/1‑150 (MQ=255)                                                                                                                                              ggTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTGAAAACAAAGAGTAGGTCAGGAAAGTATCAAACGCAGTCAGACACCATCAATGGGCTATCGCAGGCAAAACATAGCAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCGGGATTACCTTATCACTGATGGTGCCGa                                        >  4:426147/1‑200 (MQ=255)                                                                                                                                              ggTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTCTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGa                                        >  4:192792/1‑200 (MQ=255)                                                                                                                                                   cttCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGt                                                                                                                                                                        <  3:680388/66‑1 (MQ=255)                                                                                                                                                   aTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCTTTAAAAAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCATGCCAAAACTAGTAATGGAACaat                                                                                              >  1:272309/1‑140 (MQ=255)                                                                                                                                                   aTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCCATACCGCAAAGGGt                                                                                                                                                                        >  4:680432/1‑67 (MQ=255)                                                                                                                                                        cACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGaa                                                                                                  >  2:260630/1‑132 (MQ=255)                                                                                                                                                        cACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCCATAAAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCCAAAGATAGTAATGGAACAAAAGGCt                                                                                         >  1:195938/1‑141 (MQ=255)                                                                                                                                                        cACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCTTGAAAAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGtt                              >  4:832886/1‑200 (MQ=255)                                                                                                                                                         aCTGTCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAAc                                                                                                 <  2:172903/132‑1 (MQ=255)                                                                                                                                                            ggCGCTGATGCTTTCATGCTGTATCTCGTATGCTCTGGTTTCACATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGa                          <  4:808447/200‑1 (MQ=255)                                                                                                                                                             gcgcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGTTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAggg                                                                           <  3:417134/150‑1 (MQ=255)                                                                                                                                                                      cTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGTAAAGTATCAATCGCAGTCAGACACCTTCCCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGATTCGCTGGATTACCTTATCACTGATGGTGCTTATGACATCTTTAAGATGCTCAATGa                >  4:615306/1‑200 (MQ=255)                                                                                                                                                                          cAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTg                                                                                <  2:411950/132‑1 (MQ=255)                                                                                                                                                                           aGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAg           <  4:232006/200‑1 (MQ=255)                                                                                                                                                                                gTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAgggg                                                                          >  2:383652/1‑132 (MQ=255)                                                                                                                                                                                    ttcttattgtcggggttCCCATACCACGAAGGTTAGGTCTGGAAAGTATTAATCGCAGTCAAAAACCAGCACTGGGATATCACAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                    <  3:294238/135‑1 (MQ=255)                                                                                                                                                                                    ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTa                                                             >  1:283783/1‑141 (MQ=255)                                                                                                                                                                                    ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAcc  >  4:479618/1‑200 (MQ=255)                                                                                                                                                                                    ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAcc  >  4:154476/1‑200 (MQ=255)                                                                                                                                                                                    ctcGTATACTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATAACAGTCAGACAACATCACTGGGCAATCGCAGTCAAAAGATAGGAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                    <  3:109076/150‑1 (MQ=255)                                                                                                                                                                                         aTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATg                                               >  3:636222/1‑150 (MQ=255)                                                                                                                                                                                              ctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAATAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGcc                                          >  3:510430/1‑150 (MQ=255)                                                                                                                                                                                              ctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGcc                                          >  3:456316/1‑150 (MQ=255)                                                                                                                                                                                                ggTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                                                          <  2:282773/132‑1 (MQ=255)                                                                                                                                                                                                  gttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAGCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTGTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATg                                      <  3:168478/149‑1 (MQ=255)                                                                                                                                                                                                    tCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATc                                                      >  2:553979/1‑132 (MQ=255)                                                                                                                                                                                                     |                                                                                                                                                                                       ATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACC  >  NC_000913/3875333‑3875711

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