Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,875,521	0	T→C	28.8%	24.5 ‑11.4 0.0	67	A28A (GCT→GCC)	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 (19/0):  ref base (21/26):  total (40/27)
Fisher's exact test strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACC  >  NC_000913/3875323‑3875711                                                                                                                                                                                                       |                                                                                                                                                                                               gCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAAGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTATTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGctc                                                                                                                                                                                               <  4:523407/200‑1 (MQ=255)           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)                                                                ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGt                                                                                                                                                                                           >  1:550000/1‑141 (MQ=255)                                                                  aaaGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGttt                                                                                                                                                                                         <  1:70022/141‑1 (MQ=255)                                                                     gaagaGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTGCTTGACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTa                                                                                                                           >  4:569941/1‑200 (MQ=255)                                                                        gagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTgg                                                                                                                                                                                            <  2:207413/132‑1 (MQ=255)                                                                        gagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCg                                                                                                                        >  4:26462/1‑200 (MQ=255)                                                                        gagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTATACAAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCg                                                                                                                        >  4:625174/1‑200 (MQ=255)                                                                         aatATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGAGTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGc                                                                                                                       <  4:454827/198‑1 (MQ=255)                                                                          gTATGCTTTCAACGGCTTAGCTACTCGTCTAACGGACCAACTATGAAGTTTAATTTTAAGGGTTTTTTCAAGGCTGCCGGTTTATTCTCACTGGCTCTGATCTTTTCACGCGTTCTCCGGCATGCTCTGGtt                                                                                                                                                                                          >  2:358560/1‑132 (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)                                                                                                                                                                                              ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTa                                                             >  1:283783/1‑141 (MQ=255)                                                                                                                                                                                              ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAcc  >  4:154476/1‑200 (MQ=255)                                                                                                                                                                                              ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAcc  >  4:479618/1‑200 (MQ=255)                                                                                                                                                                                              ctcGTATACTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATAACAGTCAGACAACATCACTGGGCAATCGCAGTCAAAAGATAGGAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAc                                                    <  3:109076/150‑1 (MQ=255)                                                                                                                                                                                                   aTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATg                                               >  3:636222/1‑150 (MQ=255)                                                                                                                                                                                                       |                                                                                                                                                                                               GCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACC  >  NC_000913/3875323‑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.