BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,520,504:1	+T	10.6%	intergenic (‑180/+167)	yjhU ← / ← yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,520,504	1	.	T	10.6%	437.8 / 56.8	123	intergenic (‑180/+167)	yjhU/yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand): ref base . (53/57); new base T (7/6); total (60/63)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-01

Only 100 of 128 total aligned reads displayed.

GCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTGGTGGTCAGGCCAAAAAA  >  NC_000913/4520298‑4520778

gctgctGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                                                                                                                                                                                                >  4:45357/1‑275 (MQ=255)

                      tCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCGGTTTACCCGGAGAGTTTATGCAGATTGCCAGTGGCATGATTTTATACGGATATTTTATTCTTTTACATGAAATAAGCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAGGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc                                                                                                                                                                                          <  2:407968/275‑1 (MQ=255)

                                tataCCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGAAATG‑CTCCGGCTCTACAGTATTACTGAGAATACACCCATAAAGCATAATCCACGCATCAATATCGCCCCTTTTTTATCGCCAGCGATTCGTCACCGTGCCGAc                                                                                                                                                                                 >  1:144030/1‑274 (MQ=255)

                                                  gggAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAGCGCCAGCGCTCCGGCACCGTGCAGACGGTTGCCGTCCCTCACTa                                                                                                                                                               >  1:150342/1‑274 (MQ=255)

                                                        tGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATg                                                                                                                                                        <  2:45358/275‑1 (MQ=255)

                                                           tCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATg                                                                                                                                                                                                                                                                                     >  6:112949/1‑148 (MQ=255)

                                                            cATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑c                                                                                                                                                                                                                                                                                   >  5:607208/1‑148 (MQ=255)

                                                             cccTTAACCTTGTTTACCCTGATAGTTGCNGCAGATGACAAGTGGCATGATTTTCNACCGATCTTTTATTCGTTTACATGCAATAAACCAGGTTACCATTTGTGCGATCGCCCGCTGTTCCCCCATAATGTAAATCNGGAGACATG‑CGCCGGCGTTACAGTATTACTGGGAATACACCCATCAAGCATAATCGAATGTTCAATATGGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGa                                                                                                                                                   <  3:150341/273‑1 (MQ=255)

                                                              tCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                                 <  5:749202/148‑1 (MQ=255)

                                                              tCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                                 <  2:776185/148‑1 (MQ=255)

                                                               cTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTcc                                                                                                                                                                                                                                                                                >  6:265362/1‑148 (MQ=255)

                                                                 tAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                                                                                                                                                             <  7:701857/149‑1 (MQ=255)

                                                                   aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCt                                                                                                                                                                                                                                                                            >  5:277512/1‑148 (MQ=255)

                                                                    cccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                                                                                                                                                            <  5:803427/147‑1 (MQ=255)

                                                                     ccTGTTTACCCGGGGCGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAACGCCGGAGACATG‑CTCCGGCTTTACAGTAATACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaac                                                                                                                                           <  2:170168/275‑1 (MQ=255)

                                                                      cTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                                                                                                                                                                                                        >  8:219866/1‑149 (MQ=255)

                                                                       tGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                                 <  3:197969/139‑1 (MQ=255)

                                                                         tttACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                                                                                                                                                                                                             >  2:18690/1‑141 (MQ=255)

                                                                             cccGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                                                                                                                                                                                                         <  4:18690/141‑1 (MQ=255)

                                                                                ggAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                                                                                                                                                                                                              >  7:252930/1‑149 (MQ=255)

                                                                                 gagaGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTga                                                                                                                                                                                                                                                             <  8:50942/149‑1 (MQ=255)

                                                                                      ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAAAATa                                                                                                                                                                                                                                                        >  8:658319/1‑149 (MQ=255)

                                                                                                    cAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                                                                                                                                                                                                           >  5:508897/1‑148 (MQ=255)

                                                                                                     aGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGACCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                                                                                                                                                                                                          >  6:443110/1‑148 (MQ=255)

                                                                                                      gTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACAGAGAATACACCCATCAAGCATaa                                                                                                                                                                                                                                        >  7:507023/1‑149 (MQ=255)

                                                                                                          cATGATTTTATACAGATCTTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                                                                                                                                                                                                     >  5:344889/1‑148 (MQ=255)

                                                                                                           aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                                                                                                                                                                                                   <  8:43021/149‑1 (MQ=255)

                                                                                                           aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                                                                                                                                                                                                    >  6:754068/1‑148 (MQ=255)

                                                                                                            tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                                                                                                                                                  <  7:219866/149‑1 (MQ=255)

                                                                                                            tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                                                                                                                                                             <  5:78591/106‑1 (MQ=39)

                                                                                                            tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                                                                                                                                                             >  6:78591/1‑106 (MQ=39)

                                                                                                             gATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                                                                                                                                                                                                  >  6:306829/1‑148 (MQ=255)

                                                                                                                    taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                                                                                                                                                                                                          >  8:190395/1‑149 (MQ=255)

                                                                                                                        gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCac                                                                                                                                                                                                                                              >  7:577456/1‑124 (MQ=255)

                                                                                                                        gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                                                                                                                                                                                                              <  8:577456/125‑1 (MQ=255)

                                                                                                                        gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAt                                                                                                                                                                                                                                       <  7:91458/132‑1 (MQ=255)

                                                                                                                        gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACa                                                                                                                                                 <  4:173478/218‑1 (MQ=255)

                                                                                                                            ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                                                                                                                                                                                                                  >  7:6778/1‑149 (MQ=255)

                                                                                                                                aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                                                                                                                                                                                                               <  5:265362/148‑1 (MQ=255)

                                                                                                                                  tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                                                                                                                                                                                                            <  7:458681/149‑1 (MQ=255)

                                                                                                                                  tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                                                                                                                                                                                                            <  8:770577/149‑1 (MQ=255)

                                                                                                                                    ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                                                                                                                                                                                                           <  5:306829/148‑1 (MQ=255)

                                                                                                                                    ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                                                                                                                                                                                                          >  7:119055/1‑149 (MQ=255)

                                                                                                                                     tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                                                                                                                                                                                                          >  5:720674/1‑148 (MQ=255)

                                                                                                                                      ttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTAt                                                                                                                                                                                                                                                                   >  4:460505/1‑90 (MQ=255)

                                                                                                                                       tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATa                                                                                                                                                                                                                                                        <  1:179518/100‑1 (MQ=255)

                                                                                                                                           tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                                                                                                                                                                                                   >  8:905358/1‑149 (MQ=255)

                                                                                                                                            gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                                                                                                                                                                                                  >  7:376609/1‑149 (MQ=255)

                                                                                                                                               aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                                >  5:450007/1‑148 (MQ=255)

                                                                                                                                               aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                                <  6:630172/148‑1 (MQ=255)

                                                                                                                                               aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                                                                                                                                                                                               <  7:474335/149‑1 (MQ=255)

                                                                                                                                                  aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                                                                                                                                                                                                           >  5:810596/1‑134 (MQ=255)

                                                                                                                                                  aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                                                                                                                                                                                                           <  6:810596/134‑1 (MQ=255)

                                                                                                                                                   aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                                                                                                                                                                                                                 >  4:183490/1‑95 (MQ=255)

                                                                                                                                                    ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTTAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc                                                                                                                                                                                          <  8:376609/149‑1 (MQ=255)

                                                                                                                                                      aGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCg                                                                                                                                                                                        <  7:658319/149‑1 (MQ=255)

                                                                                                                                                       gCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                                                                                                                                                                                                             <  2:183491/95‑1 (MQ=255)

                                                                                                                                                        cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAGAGCCGGAGACATG‑CTTCGGCTTTACAGTATTACTGAGAATATACCCATCCAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACCGTTGCCGTCACTCACTAAATTATTGGGGACAGTGAACAACAGGTTGGTCCACTGGCAAGAACCCGATTCTAAAACTAAGCTGGCGAGCAGAATAAAAACAAGCCCACTGATGGaata                                                           >  4:681139/1‑269 (MQ=255)

                                                                                                                                                        cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACCCCCATCACGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGCCGGTTGCCGTCACTGACCAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTACGCTGGCCAGCCGAATATAAATCAGCCCACTGATGGATAagc                                                         >  2:288352/1‑273 (MQ=255)

                                                                                                                                                             cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATaa                                                                                                                                                                                                                                        >  2:1513/1‑94 (MQ=255)

                                                                                                                                                               tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACgg                                                                                                                                                                               >  8:198930/1‑149 (MQ=255)

                                                                                                                                                                    gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGcc                                                                                                                                                                          <  7:905358/149‑1 (MQ=255)

                                                                                                                                                                      gATCGCCCGCTGTGCCCCCATACTGGAAAGCCGGAGACATG‑CTCAGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATTATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCAGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTGACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCg                                          <  2:325194/275‑1 (MQ=255)

                                                                                                                                                                       aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACCCCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATaa                                                           >  2:136124/1‑257 (MQ=255)

                                                                                                                                                                       aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGt                                                                                                                                                                        >  5:647438/1‑148 (MQ=255)

                                                                                                                                                                           cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAgt                                                       <  4:136123/257‑1 (MQ=255)

                                                                                                                                                                             cGCTGTTCCCCCATACCGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                                                                                                                                                  >  5:399754/1‑148 (MQ=255)

                                                                                                                                                                              gCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                                                                                                                                                                                                >  7:713497/1‑101 (MQ=255)

                                                                                                                                                                                 gTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTaa                                                                                                                                                              <  6:450007/148‑1 (MQ=255)

                                                                                                                                                                                      cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTAt                                                                                                                                                         <  6:607208/148‑1 (MQ=255)

                                                                                                                                                                                      cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTAt                                                                                                                                                         <  6:277512/148‑1 (MQ=255)

                                                                                                                                                                                      cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATg                                                                                                                                                        <  8:6778/149‑1 (MQ=255)

                                                                                                                                                                                        cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATggg                                                                                                                                                      <  8:119055/149‑1 (MQ=255)

                                                                                                                                                                                         aTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATggg                                                                                                                                                      <  6:720674/148‑1 (MQ=255)

                                                                                                                                                                                               tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCCCTTAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCg                                                                                                               >  2:777787/1‑181 (MQ=255)

                                                                                                                                                                                               tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                                                                                                                                                       >  4:780260/1‑141 (MQ=255)

                                                                                                                                                                                               tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAgg                                    >  2:196416/1‑256 (MQ=255)

                                                                                                                                                                                               tAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAg                                                                                                                                                >  5:396764/1‑148 (MQ=255)

                                                                                                                                                                                                 aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                                                                                                                                                                                                           >  4:88693/1‑71 (MQ=255)

                                                                                                                                                                                                 aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTtggtgg               >  2:622902/1‑275 (MQ=255)

                                                                                                                                                                                                 aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTtggtgg               >  2:538276/1‑275 (MQ=255)

                                                                                                                                                                                                  aGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                                                                                                                                                       >  8:227765/1‑138 (MQ=255)

                                                                                                                                                                                                  aGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGa                                                                                                                                             >  6:709236/1‑148 (MQ=255)

                                                                                                                                                                                                   gCCGNAGACAGG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAGCAAGCATAATCGAAGGTTCAATATCGCCCCTTTATGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTGGGCTCACTGGCAAGAACCCGATta                                                                                                           <  4:777785/181‑1 (MQ=255)

                                                                                                                                                                                                   gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGa                                                                                                                                                   <  2:780262/141‑1 (MQ=255)

                                                                                                                                                                                                   gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAAt                                                                             <  1:102581/211‑1 (MQ=255)

                                                                                                                                                                                                   gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATaa                                <  4:196415/256‑1 (MQ=255)

                                                                                                                                                                                                   gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa                                                                                                                                            >  6:875222/1‑148 (MQ=255)

                                                                                                                                                                                                   gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa                                                                                                                                            >  5:876140/1‑148 (MQ=255)

                                                                                                                                                                                                    ccGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTgg                                                                                             >  2:744629/1‑194 (MQ=255)

                                                                                                                                                                                                     cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa                                                                                                                                            >  7:317836/1‑146 (MQ=255)

                                                                                                                                                                                                     cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa                                                                                                                                            <  8:317836/146‑1 (MQ=255)

                                                                                                                                                                                                       gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcac                                                                                                                                                                     <  5:452902/119‑1 (MQ=255)

                                                                                                                                                                                                        agaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAacaaca                                                                                                                                       >  6:458082/1‑148 (MQ=255)

                                                                                                                                                                                                           cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGt                                                                                                                                    <  6:344889/148‑1 (MQ=255)

                                                                                                                                                                                                           cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGt                                                                                                                                    >  5:30769/1‑148 (MQ=255)

                                                                                                                                                                                                           cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTc                                        <  2:465456/240‑1 (MQ=255)

                                                                                                                                                                                                            aTG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTg                                                                                                                                  >  8:90440/1‑149 (MQ=255)

                                                                                                                                                                                                              g‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                                                                                                                                                       >  1:119478/1‑126 (MQ=255)

                                                                                                                                                                                                              g‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGTACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTACAAAACTAAGCTGGCGAGCAGAATAAAACACAGCCCACTGATGGATAAGAGTGTGGAAATAAGCGTCAAGGATAACAGCGTTTATTTGGTGGTCAGGCCaaaaaa  >  4:660062/1‑275 (MQ=255)

GCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTGGTGGTCAGGCCAAAAAA  >  NC_000913/4520298‑4520778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`