Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,520,500	0	A→.	10.4%	11.5 ‑1.1 ‑0.1	48	intergenic (‑176/+171)	yjhU/yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (22/21):  total (23/25)
Fisher's exact test strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.89e-01

GAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGC  >  NC_000913/4520292‑4520707                                                                                                                                                                                                                 |                                                                                                                                                                                                                 gAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAc                                                                                                                                                                                                                 <  2:201595/210‑1 (MQ=255)     aGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAACCCGAAGACATG‑c                                                                                                                                                                                                            >  2:77345/1‑210 (MQ=255)        ctgctgTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAGT‑CTcc                                                                                                                                                                                                         >  2:238109/1‑210 (MQ=255)            tgtttgttTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCt                                                                                                                                                                                                     <  2:517007/210‑1 (MQ=255)             gtttgttgGATCTCTAGCCACTTTTGTATACCTGTATTAAACCGGCGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCGCATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTCGAGCGATTGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                                                                                                                                    <  2:247528/210‑1 (MQ=255)              tttgtttgATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                                                                                                                                   >  1:368993/1‑210 (MQ=255)                 gtttgATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                                                                                                                                <  1:408743/210‑1 (MQ=255)                  tttgATCTCTATCCACTTTGGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                                                                                                                               <  1:87312/210‑1 (MQ=255)                   ttgATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGt                                                                                                                                                                                              <  2:514145/210‑1 (MQ=255)                    tgATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                                                                                                                             <  2:582280/210‑1 (MQ=255)                     gATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTAt                                                                                                                                                                                            <  2:301743/210‑1 (MQ=255)                       tctcTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGTCTTTACAGTAttt                                                                                                                                                                                          >  1:488443/1‑209 (MQ=255)                        ctctATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                                                                                                                                         <  1:414452/210‑1 (MQ=255)                               cACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAAt                                                                                                                                                                                  >  2:100097/1‑210 (MQ=255)                                aCTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATa                                                                                                                                                                                 <  1:183447/210‑1 (MQ=255)                                  ttttGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                                                                                                                               <  2:215905/210‑1 (MQ=255)                                     tGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                                                                            >  1:277654/1‑210 (MQ=255)                                                   aaaCCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CCCCGGCTTTACAGTATTACTGAGAATACACCCATGCAGCATAATCg                                                                                                                                                              >  1:468410/1‑210 (MQ=255)                                                    aaCCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCTCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATcta                                                                                                                                                             >  2:514585/1‑208 (MQ=255)                                                      ccGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                                                                           >  1:555571/1‑210 (MQ=255)                                                                     tCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcccc                                                                                                                                            >  2:395781/1‑210 (MQ=255)                                                                               gTTTACACGAAGAGTTGCTGGAGATTGTCCGTGGCATGATTTAATACAGATATTTTATTCTTATACATGAAATAACCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCAACATACTGTAAAGCCGGAGACATG‑CGCCGGATTTACAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                                                                                                                                  <  2:123491/210‑1 (MQ=255)                                                                                  tACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTACGGCTTTACAGTATTACTGAGAATACAACCATCAAGCATAATCAAAGGTTCAATATCGCCCCTTTTTGATCGCCa                                                                                                                               >  1:496959/1‑210 (MQ=255)                                                                                          gagTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg                                                                                                                       <  2:316149/210‑1 (MQ=255)                                                                                           agTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt                                                                                                                      <  2:284619/210‑1 (MQ=255)                                                                                            gTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc                                                                                                                     <  2:61633/210‑1 (MQ=255)                                                                                                         gCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACgg                                                                                                        >  2:419881/1‑210 (MQ=255)                                                                                                         gCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATAGCCCGCTGTTCCCCCATAATGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATAAACCCATCAAGCATAATCGAAGGTTCAATTTCGCCCATTTTTGATCGCCAGCGATTAGTCACCGTGCAGACgg                                                                                                        >  2:531727/1‑210 (MQ=255)                                                                                                                gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGt                                                                                                 >  1:162334/1‑210 (MQ=255)                                                                                                                 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtc                                                                                                >  1:362875/1‑210 (MQ=255)                                                                                                                            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAAt                                                                                     >  2:217524/1‑210 (MQ=255)                                                                                                                                   ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCTTCACCGTGAAGCCGGTTGCCGTACCCCCCTAAATT‑Tggggg                                                                             >  2:392021/1‑210 (MQ=255)                                                                                                                                         tCTTTTACATCAAATAAACCAGCAAACCATTGGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACTCCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGt                                                                        <  1:298051/210‑1 (MQ=255)                                                                                                                                              tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaaca                                                                   <  1:338036/210‑1 (MQ=255)                                                                                                                                              tACATGAAATAAACCAGCATACCATTTGAGCGATCGCACGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaaca                                                                   <  2:182623/210‑1 (MQ=255)                                                                                                                                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaacaac                                                                 >  2:343324/1‑210 (MQ=255)                                                                                                                                                            cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCAc                                                     >  1:376533/1‑210 (MQ=255)                                                                                                                                                                aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGc                                                 >  1:417408/1‑210 (MQ=255)                                                                                                                                                                 tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCa                                                <  1:501612/210‑1 (MQ=255)                                                                                                                                                                     aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGaa                                            >  2:140239/1‑210 (MQ=255)                                                                                                                                                                                 gcccgatctcaccccTTTCTGTAGAGCCGGAGACACG‑CTCCGCGTATACAGTATTACTCATACTACACCCATCAAGCATTATCAAACGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATaaa                                <  1:59120/199‑1 (MQ=255)                                                                                                                                                                                                 tACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGCGCa                >  2:460193/1‑210 (MQ=255)                                                                                                                                                                                                    tGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGaa             >  1:466540/1‑210 (MQ=255)                                                                                                                                                                                                      tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATa           <  2:137666/210‑1 (MQ=255)                                                                                                                                                                                                        aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATaaa         <  1:395745/210‑1 (MQ=255)                                                                                                                                                                                                          gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATaaaaa       <  1:460152/210‑1 (MQ=255)                                                                                                                                                                                                          gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATaaaaa       <  2:466583/210‑1 (MQ=255)                                                                                                                                                                                                             ggAGACATG‑CTCCGGCTTTACAGTATTACTGAGCATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCa    <  2:432334/210‑1 (MQ=255)                                                                                                                                                                                                               agaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGc  <  1:391985/210‑1 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                                 GAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGC  >  NC_000913/4520292‑4520707

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