BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,344	0	C	T	38.4%	616.8 / 241.1	400	S62F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand): ref base C (15/230); new base T (153/0); total (170/230)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 410 total aligned reads displayed.

GATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAG  >  NC_000913/4604097‑4604591

                 gTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTACTTTTTTTTTTTTTTGCCATgtgatatg                                                                                                                                                                                                                                      >  1:5016850/1‑247 (MQ=255)

                    cAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTTCCTGAGAagagagg                                                                                                                                                                                                                                   >  1:804430/1‑250 (MQ=255)

                                         ttGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGcc                                                                                                                                                                                                              <  2:671124/250‑1 (MQ=255)

                                           gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGGGAAGTGAGAGCCGAGAgggggtatt                                                                                                                                                                                                                   >  2:4616900/1‑238 (MQ=255)

                                             ttATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                                                                                                                          <  2:1250672/251‑1 (MQ=255)

                                                        ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTTCCCTGTGAAGTGAGGGCCGAGAgggggtatttcttttttttctaataat                                                                                                                                                                                                 >  1:5759483/1‑225 (MQ=255)

                                                                              aTTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                                                                         <  1:1164051/251‑1 (MQ=255)

                                                                                                                  gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTTCCTTACTGGACCGGCg                                                                                                                                                                                                >  2:4530337/1‑192 (MQ=255)

                                                                                                                    cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                                                                                                                                 <  2:2557474/237‑1 (MQ=255)

                                                                                                                       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATTAGAAGTGAGGGCag                                                                                                                                                                                                                               >  1:1903716/1‑156 (MQ=255)

                                                                                                                       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGTGAAGTGAGCGCGGGGAGGGATTATTTTTGCTGACTGGACGGGCGATTAAGTTTCCGCGTACCCGGCGGTTGTTTGTTTTGGGTGATGTTATTGagaggcggccggttg                                                                                                                                >  2:5094913/1‑237 (MQ=255)

                                                                                                                       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGaaagg                                                                                                                                                                                                                                      >  1:2939476/1‑146 (MQ=255)

                                                                                                                       gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGGCGAGAggggt                                                                                                                                                                                                                       >  1:585098/1‑162 (MQ=255)

                                                                                                                        agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt                                                                                                                                                                                                                                                        >  2:1684767/1‑130 (MQ=255)

                                                                                                                        agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                                                                                                                                        <  1:1685072/130‑1 (MQ=255)

                                                                                                                         gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                                                                                                            <  1:4801841/237‑1 (MQ=255)

                                                                                                                            gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTATTTTGCCTGATTGAACCGGCGATTTAGGTTTCGCGTACCCGGCGTTTTGTTATTGTGGATGTTGATATTTTAGGTCGGCT‑ggttggtttt                                                                                                                          >  2:5264089/1‑249 (MQ=255)

                                                                                                                            gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGAGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGGACCCGGCGGTTTGTTATTTCGGGTGGTGATATTTAAGGACGGCT‑ggttggtt                                                                                                                            >  1:353143/1‑249 (MQ=255)

                                                                                                                                 ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                                                                                                            <  2:3322642/229‑1 (MQ=255)

                                                                                                                                       tagagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATGCCAGTCAGGAGTCATTTCAGGATGCCATGGCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATc                                                                                                               <  1:5308922/250‑1 (MQ=255)

                                                                                                                                            tgtATTAAGAAAGTTTGGGGGTAACCCTTATAAGTGCCATACATATGCCGGTGAGAGGTCTTTGCAGGAGGCAAGGTGGGGGATGTGGTTTGGGGGGGTAATTTTTTCTTTTTGTGCCAGGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCGGGCGGAACTGGGGATAAAGTTCCCGCGTCCCCGGCGTTTAGTTATTGCGGAGGATGATATTGAAGCGCGGCT‑GATTGGTTCATTGTCGCCATCACCGc                                                                                                          <  1:2581404/249‑1 (MQ=255)

                                                                                                                                                             cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGATTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  1:3722367/251‑1 (MQ=255)

                                                                                                                                                             cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  2:4330445/251‑1 (MQ=255)

                                                                                                                                                             cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  2:2476918/251‑1 (MQ=255)

                                                                                                                                                             cGGGTAACCCTTATAAGCTACATACCTATACCAGTCAGGAGTAATTTCAGGATGCCATGTCGCTGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAg                                                                                          <  1:4749442/250‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATCCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGACGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  2:3980798/244‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGGTATCTTTCCTGACTGAAATgggg                                                                                                                                                                                                >  2:5446241/1‑140 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGGGGGGGGTTATCTTGCCTGGCTGAAATGGGGGTTAAGTTTTCGCGGAACCGGCGGTTTGTTATTg                                                                                                                                                              >  1:618766/1‑176 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGGCTGTACTGGCGATTAAGTTTCCGCGGACCCGGCGGTTTGTTATTGGGGATGATGATATTTAAGCTCGGCT‑GGTTGGTTTATTGTTTCCATCACCGCTGGGAGGTGGATTATGGTAAGCGt                                                                                  >  2:3520392/1‑251 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCGATTAAATTTTCGCGTACCCCGCGTTTAGTTATTTCGGGTGATGATATTTTAGTTCCGGT‑GGTTTGGTTATTTTTGCCATCCCCCCTGGGCGGGGGAtttaatt                                                                                        >  2:3660290/1‑240 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGGGAGGGGTTATCTTGCCTGTCTTAACTGGCGATTAAGATTTCGCGTAACCGGCGGTTAGGTATTTTGGGTGGTGATATTTAAGCTCGGCT‑GGTTGGGTTATTTTTTCCATCACCGCTGGACGGGGGATTTAGTAAAGAGt                                                                                  >  1:3161584/1‑251 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGGTTTCTTTCCTGGCTGAACTgggg                                                                                                                                                                                                >  1:1588439/1‑140 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                                                                                                                                                                                                  <  1:4670044/124‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCg                                                                                                                                                                                                <  2:5638054/142‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                                                                         <  1:178115/165‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga                                                                                                                                                    <  1:3296137/186‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGAc                                                                                                     <  2:5822099/232‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                       <  1:4973431/246‑1 (MQ=255)

                                                                                                                                                                    cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGACGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         >  1:3981457/1‑244 (MQ=255)

                                                                                                                                                                     ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTTAGCGCCGGGAGGGGTTTTTTTTCCTGACTGAAATGGCGGTTAAGTTTTCGCGTACCCGGCGGTTTGTTTTTGCGGGTGGTGGTATTGAAGATCGGCC‑TGTTTGTTCCTTGTTGCCCTTCCCGCTggggcgggtatttagg                                                                                         >  2:2858675/1‑229 (MQ=255)

                                                                                                                                                                     ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaggag                                                                                                                                                                                                                                        >  2:4882306/1‑101 (MQ=255)

                                                                                                                                                                     ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                                                                                                                        <  1:4883141/101‑1 (MQ=255)

                                                                                                                                                                     ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  1:2859131/243‑1 (MQ=255)

                                                                                                                                                                     ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCCTTTTTTCTTTTTTTGCCCTGACGAGGGAGGGCCGAGGGGGGTTATCTTGCCTGACTGTACTTGCGGATAAGTTTCCCCGGCCCCGGCGTTTTGTTTTTTGGGATGATGATATTTAAGGTCGGGTGGTTTGGTGCATT‑TCGTCATCGCCTCTCGCCggcgggttttggttgtgttt                                                                                 >  2:3118930/1‑233 (MQ=255)

                                                                                                                                                                      cTTATAAGCTCCATACCTATGCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                                                                                                                      <  1:1385645/198‑1 (MQ=255)

                                                                                                                                                                      cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCTTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTgg                                                                                                                              <  2:1739335/205‑1 (MQ=255)

                                                                                                                                                                      cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGGGGGCGCCGGGGGGGGGTATCTTGCCTGGCTGAACTGGGGGTTTAGTTTTCCCGTATCCGGCGGTTTGTTATTTTGGGTGGTGATATTGAAGCTCGGGT‑GGTTGGGTTTTTTTTGCCATCACCGCTGGGCGGGGGTGTTAGTTAAGCGTc                                                                                 >  2:5571788/1‑250 (MQ=255)

                                                                                                                                                                      cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGGGCGCAGGGAGGGGGTATTTTGCCTGGATGGAATGGGGATTAAGGTTTCGCGGACCCGGCGGTTTGGTATTGCGGGGGATGGTATTTAAAAGCGGCT‑GGTTgggtt                                                                                                                           >  1:1304377/1‑205 (MQ=255)

                                                                                                                                                                      cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                                                                                                                        <  2:5939411/100‑1 (MQ=255)

                                                                                                                                                                      cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGcc                                                                                                                                                                                                              <  2:5740406/126‑1 (MQ=255)

                                                                                                                                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGTAGGGCAGAGAGGGGTTTTTTTTCCTGTCTGAACTGGCGGTTAAGTTTTCGCGGACCCCGCGTTTTTTTTTTTGGGGTGTTGATATTTAAGCTCGGGT‑GGTTGGGTTGTTTTCGCCATCACCCCttggccgggtgtttatg                                                                                         >  2:4368100/1‑219 (MQ=255)

                                                                                                                                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAg                                                                                                                                                                                          <  1:1830183/139‑1 (MQ=255)

                                                                                                                                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCATGAGGAGTGAGCGCGGCGGGGGGTTATTTTTCCTGTCTGAAATGGCGGTTAAGTTTTCGCGGAACCGGCGGTTTGTTTTTGTGGGTGGTGATATTGAAAGACGGCT‑GATTTGTTCATTGTTGCGCTCAACGCCGGACGGGGTGTTAAGTGAAGCGTCAACg                                                                             >  2:4736039/1‑247 (MQ=255)

                                                                                                                                                                             gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTTCCATGGGGAGGGAGCGGCGAGAGGGGGTTTCTTTCCTGGCTTGACTGGCGGTTAAATTTTCCCGGAACCCGCGGTTTTTTTTTTTGGGTGATGATAATTAAGATCGGGTGGTTTGGTTCTTTTTGGCCacaaccccggga                                                                                                      >  1:3462835/1‑212 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCTGGATTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                       <  1:2228119/236‑1 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGTTCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAGGGGGTTTCTTTCCTGTCTGAACTGGCGGTTAAGTTTCCGCGGACCCGGCGGTTTGTTATTTGGGGTGGTGATATTGTAGCTCGGGT‑GTTTGGGTTTTTTTTGCCATCACCCCTGGATGGGGTGTTTAGTAAAACGTTAACgtgggg                                                                        >  2:3406595/1‑246 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGGGAAAGGAGCGCAGGGGGGGGTTATTTTTCCTTGCTTGACTGGGGGTTTAATTTCCGCGGAACCCGCGGTTTGGTTTTTCGGAGGGTGATATTGAAGcgtggcgggtttg                                                                                                                               >  2:4111598/1‑185 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCCGAGAGGGGTTTTCTTTCCTGTCTGAACTGGCGGTTAAGTTTTCCCGGACCCCGCGGTTTGTTTTTTCGGGTGATGATATTTAAACTCGGCC‑GATTTGTTTAATGTTTCCATTACCGCTGGACGGGGGCTTTAGTAAAACGTGCAAGCTggg                                                                        >  1:5714037/1‑250 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                                                                                                                                                                                                  <  1:2504267/114‑1 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTc                                                                                                                           <  1:59289/200‑1 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGAc                                                                                                     <  1:3905793/222‑1 (MQ=255)

                                                                                                                                                                              cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCNCGNNTGCGGCGNNNNTTNTTTCTTTTNCNGCCATGAGAAGTGAGCGCAGAGAGGGATTAtn                                                                                                                                                                                                                   <  2:7817/113‑2 (MQ=255)

                                                                                                                                                                               tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                                                                                                                            <  1:767826/183‑1 (MQ=255)

                                                                                                                                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCATGAGGAGGGAGCGGGGGGGGGGGGTTTTTTTTCTGGCCGAACTGGGGGTTAAGTTTTcgcggaa                                                                                                                                                                               >  2:1802918/1‑141 (MQ=255)

                                                                                                                                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                                                                                                                          <  2:1599802/85‑1 (MQ=255)

                                                                                                                                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                                                                                                                                           <  2:2772459/100‑1 (MQ=255)

                                                                                                                                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                                                                                                                                                                                            <  2:2169062/99‑1 (MQ=255)

                                                                                                                                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAg                                                                                          <  2:4824340/228‑1 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGGTATTTTGCCTGACTGAAATGGGGGGTAAAGTTTCGCGGAACCCGCGGTTTGTTATTTCGGGTGATGATATTTGAGCTCGGGC‑GGTTGGTTTATTTTTGCCCTCACCGCTTGGCGGGGTGTTTAGTAAAACGGCAAAGCTGGGGAtttt                                                                  >  1:524467/1‑251 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGGGAGGGGTTTTCTTTCCTGACTGAACTGGCGGTTAAGGTTTCGCGTTCCCCGCGGTTTGTTTTTTCGGGTGTTTGGATTTAAGCTCGGCC‑GATTTGGTTTTTGTCTCCCTTACCCCTGGGCGGGGTAt                                                                                              >  2:4974359/1‑223 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGGCATGAGAAAGGAGGGCcgaga                                                                                                                                                                                                                            >  1:2775932/1‑98 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                                                                                                                        <  1:5514266/86‑1 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                                                                                                                                                                                                  <  1:421644/108‑1 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGt                                                                                                                      <  1:4361380/199‑1 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATTGAAGCTCGGCT‑GATTGGTTCATTTTCGTCATCACCGCTGGGCGGTGTATTTAGTaa                                                                                       >  2:1281255/1‑230 (MQ=255)

                                                                                                                                                                                    aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGAAAGGAGCGCAGAGAGGGGGTATCt                                                                                                                                                                                                                  >  2:421569/1‑108 (MQ=255)

                                                                                                                                                                                     ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAgcgc                                                                                                                                                                                                                                 >  1:1266909/1‑92 (MQ=255)

                                                                                                                                                                                     ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCCGAGAGGGATTAt                                                                                                                                                                                                                    >  2:4176744/1‑105 (MQ=255)

                                                                                                                                                                                     ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                                                                                                                       <  2:4433880/102‑1 (MQ=255)

                                                                                                                                                                                     ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                                                                                                                    <  1:4177462/105‑1 (MQ=255)

                                                                                                                                                                                     ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCGGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGATCTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCGCGGCT‑GAtt                                                                                                                                <  1:4315529/188‑1 (MQ=255)

                                                                                                                                                                                      cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatgat                                                                                                                                                   <  1:4392895/169‑1 (MQ=255)

                                                                                                                                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAggggt                                                                                                                                                                                                                       >  2:5889712/1‑97 (MQ=255)

                                                                                                                                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                                                                                                                   <  2:1850805/151‑1 (MQ=255)

                                                                                                                                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                       <  1:992274/226‑1 (MQ=255)

                                                                                                                                                                                           ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAg                                                                                                                                                                                          <  1:3380834/125‑1 (MQ=255)

                                                                                                                                                                                             aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCCGGGAGGGTTTATTTTGTCTGACTGAACTGGGGGTTTAAATTTCGCGGAACCGGCGGTTTGGTTTTGCGGGTGATGATATTGAAGCCCCGGT‑GATTGGTTTATTTTCGCCATTACCGCTGGACGGGGTATTTAgg                                                                                         >  1:2375803/1‑218 (MQ=255)

                                                                                                                                                                                             aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTc                                                                                                                           <  1:5629244/185‑1 (MQ=255)

                                                                                                                                                                                             aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  2:2375396/219‑1 (MQ=255)

                                                                                                                                                                                                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                                                                                                                       <  1:1975110/88‑1 (MQ=255)

                                                                                                                                                                                                       cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                                                                                                                        <  2:1979743/67‑1 (MQ=255)

                                                                                                                                                                                                       cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                                                                                                                       <  2:3066429/84‑1 (MQ=255)

                                                                                                                                                                                                         tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGAGGGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTTCGCGGACCCCGCGGTTTGTTTTTTCGGGTGGTGATATTTTAGCTCGGCT‑GGTTGGGTTCTTGTCGCCCTCTCCCCGGGGCGGGGGATTAAGTAAAACGTCAAAGCTGGAGATTTTTTATCAAGAACTCtttttttt                                             >  1:4518781/1‑248 (MQ=255)

                                                                                                                                                                                                                aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGGACCCGGCGGTTTGTTATTGCGGGTGATGATATTTAAGCTCGGCT‑GAtttggtt                                                                                                                           >  1:4150032/1‑161 (MQ=255)

                                                                                                                                                                                                                         cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt                                                                                         <  2:1512346/191‑1 (MQ=255)

                                                                                                                                                                                                                             gATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGAc                                                                                                     <  2:375444/175‑1 (MQ=255)

                                                                                                                                                                                                                               tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcg                       <  1:4570915/251‑1 (MQ=255)

                                                                                                                                                                                                                                  cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTc                                                                                                                           <  2:2910054/148‑1 (MQ=255)

                                                                                                                                                                                                                                     ttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgac                 <  2:3877151/251‑1 (MQ=255)

                                                                                                                                                                                                                                                   ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGGCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCa   <  1:3075734/251‑1 (MQ=255)

GATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCT‑GATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAG  >  NC_000913/4604097‑4604591

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 6 ≤` `ATCG/ATCG` `< 9 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`