BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,520,504:1	+T	8.1%	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	8.1%	728.6 / 65.3	185	intergenic (‑180/+167)	yjhU/yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand): ref base . (83/87); new base T (7/8); total (90/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

Only 100 of 187 total aligned reads displayed.

GTCTGGCTTTTGCCTTCCTGATAGTACAGCCATGTTGCATGAAGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTGGTGGTCAGGCCAAAA  >  NC_000913/4520233‑4520776

gTCTGGCTTTTGCCTTCCTGATAGTACAGCCATGTTGCATGAAGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                               <  4:310805/275‑1 (MQ=255)

                                                                                                                              aTCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTCTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                               <  6:888128/149‑1 (MQ=255)

                                                                                                                                   aaCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCt                                                                                                                                                                                                                                                                          <  5:190710/149‑1 (MQ=255)

                                                                                                                                   aaCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                                                                                                                                                           <  7:145613/148‑1 (MQ=255)

                                                                                                                                    aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                                                                                                                                                                                                         <  5:699417/149‑1 (MQ=255)

                                                                                                                                     cccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACGGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                                                                                                                                                                                                        >  6:296898/1‑149 (MQ=255)

                                                                                                                                     cccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                                                                                                                                                                                                        >  5:467199/1‑149 (MQ=255)

                                                                                                                                     cccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                                                                                                                                                                                                        >  5:846459/1‑149 (MQ=255)

                                                                                                                                     cccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                                                                                                                                                                                                        <  6:846459/149‑1 (MQ=255)

                                                                                                                                      ccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                                                                                                                                                                                                        <  8:760186/148‑1 (MQ=255)

                                                                                                                                      ccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                                                                                                                                                          <  5:618033/146‑1 (MQ=255)

                                                                                                                                      ccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                                                                                                                                                          >  6:618033/1‑146 (MQ=255)

                                                                                                                                          tttACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                                                                                                                                                                                                      >  2:742090/1‑146 (MQ=255)

                                                                                                                                          tttACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAg                                                                                                                                                                                                                                                                    <  8:743059/148‑1 (MQ=255)

                                                                                                                                          tttACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGt                                                                                                                                                                                                                                                                   >  5:734349/1‑149 (MQ=255)

                                                                                                                                             aCCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                                                                                                                                                                                                     >  5:93127/1‑144 (MQ=255)

                                                                                                                                             aCCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                                                                                                                                                                                                     <  6:93127/144‑1 (MQ=255)

                                                                                                                                              cccGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                                                                                                                                                                                                  <  4:742090/146‑1 (MQ=255)

                                                                                                                                                cGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                                                                                                                                                                                                              >  8:904115/1‑148 (MQ=255)

                                                                                                                                                       ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGGAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAAt                                                                                                                                                                                                                                                       >  7:716404/1‑148 (MQ=255)

                                                                                                                                                          aTGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATacac                                                                                                                                                                                                                                                   >  6:518106/1‑149 (MQ=255)

                                                                                                                                                            gctgATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                                                                                                                                                 >  6:706924/4‑149 (MQ=255)

                                                                                                                                                            gCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                                                                                                                                                                                                  >  8:20495/1‑148 (MQ=255)

                                                                                                                                                                   gCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                                                                                                                                                                                                          <  6:467199/149‑1 (MQ=255)

                                                                                                                                                                     cAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                                                                                                                                                                                                        <  6:758395/149‑1 (MQ=255)

                                                                                                                                                                            aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                                                                                                                                                                                                 >  6:381164/1‑149 (MQ=255)

                                                                                                                                                                             tGATTTTATACAGATATTTTATTCTTTTACATGGAACAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                                                                                                                                                                                                 <  8:716404/148‑1 (MQ=255)

                                                                                                                                                                                ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGtt                                                                                                                                                                                                                             >  6:855334/1‑149 (MQ=255)

                                                                                                                                                                                ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGtt                                                                                                                                                                                                                             <  5:809196/149‑1 (MQ=255)

                                                                                                                                                                                 tttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGtt                                                                                                                                                                                                                             <  8:161237/148‑1 (MQ=255)

                                                                                                                                                                                       cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                                                                                                                                                                                                       <  8:127259/148‑1 (MQ=255)

                                                                                                                                                                                         gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATa                                                                                                                                                                                                                                       <  7:241648/130‑1 (MQ=255)

                                                                                                                                                                                         gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATa                                                                                                                                                                                                                                       >  8:241648/1‑130 (MQ=255)

                                                                                                                                                                                          atatTTTATTCTTTTACATGAAATAAACCAGCATACCGTTTGAGCGATCGCCAGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGc                                                                                      <  2:156950/274‑1 (MQ=255)

                                                                                                                                                                                          atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactc                                                                                                                                                                 >  2:783424/1‑199 (MQ=255)

                                                                                                                                                                                           tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                                                                                                                                                                                                   <  8:978153/148‑1 (MQ=255)

                                                                                                                                                                                            atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                                                                                                                                                                                                  >  7:545580/1‑148 (MQ=255)

                                                                                                                                                                                              tttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                                                                                                                                                           <  8:478499/89‑1 (MQ=38)

                                                                                                                                                                                              tttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                                                                                                                                                                                                           >  7:478499/1‑89 (MQ=38)

                                                                                                                                                                                              tttATTCTTTTACATGAAATAAACCAGCAGACCATTGGNGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCTTCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTa                                                                                                                                                             <  4:783424/199‑1 (MQ=255)

                                                                                                                                                                                               ttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                                                                                                                                                                                              <  5:381164/149‑1 (MQ=255)

                                                                                                                                                                                               ttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtt                                                                                                                                                                                                               <  7:590825/148‑1 (MQ=255)

                                                                                                                                                                                                tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                                                                                                                                                                                              <  8:624916/148‑1 (MQ=255)

                                                                                                                                                                                                tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                                                                                                                                                                                              <  6:809820/148‑1 (MQ=255)

                                                                                                                                                                                                tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTACCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACGGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                                                                                                                                                                                                              >  5:809820/1‑148 (MQ=255)

                                                                                                                                                                                                     ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                                                                                                                                                                                                        >  6:440605/1‑149 (MQ=255)

                                                                                                                                                                                                           aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                                                                                                                                                                                                  >  6:1039556/1‑149 (MQ=255)

                                                                                                                                                                                                            tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCTCCTTTTTGATCGCCAGc                                                                                                                                                                                                  >  7:874462/1‑148 (MQ=255)

                                                                                                                                                                                                               aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGGTCAATATCGCCCCTTTTTGCTCGCCAGCGAt                                                                                                                                                                                               >  7:142113/1‑148 (MQ=255)

                                                                                                                                                                                                                aTAAACCAGCATACCATTTGAGTGATTGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                                                                                                                                                                                             >  6:839490/1‑149 (MQ=255)

                                                                                                                                                                                                                aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                              <  7:859369/148‑1 (MQ=255)

                                                                                                                                                                                                                aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                              <  7:508333/148‑1 (MQ=255)

                                                                                                                                                                                                                aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                                                                                                                                                                                             <  6:758219/149‑1 (MQ=255)

                                                                                                                                                                                                                 tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                                                                                                                                                                                             >  7:358292/1‑148 (MQ=255)

                                                                                                                                                                                                                    aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCCAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGAt                                                           >  4:538990/1‑275 (MQ=255)

                                                                                                                                                                                                                        gCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                                                                                                                                                                                                           >  8:384823/1‑127 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAGAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAaga                                                       >  2:636544/1‑274 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                              <  3:108180/139‑1 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                              >  1:108180/1‑139 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactcac                                                                                                                                                               >  2:464750/1‑170 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCAATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCTTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGtgtg                                                  >  1:15564/1‑279 (MQ=255)

                                                                                                                                                                                                                          aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAgt                                                     <  2:538990/275‑1 (MQ=255)

                                                                                                                                                                                                                           tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa                                                                                                                                                                                  >  5:1067557/1‑149 (MQ=255)

                                                                                                                                                                                                                             ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTa                                                                                                                                                                                                                                                               <  7:388984/70‑1 (MQ=255)

                                                                                                                                                                                                                             ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAg                                                                                                                                                                                 >  8:875060/1‑148 (MQ=255)

                                                                                                                                                                                                                             ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTaaa                                                                                                                                                           <  4:464750/170‑1 (MQ=255)

                                                                                                                                                                                                                              cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc                                                                                                                                                                               >  6:867440/1‑149 (MQ=255)

                                                                                                                                                                                                                                 ttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACgg                                                                                                                                                                             >  7:630768/1‑148 (MQ=255)

                                                                                                                                                                                                                                  tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGt                                                                                                                                                                      >  4:784908/1‑154 (MQ=255)

                                                                                                                                                                                                                                     gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGc                                                                                                                                                                         <  7:20495/148‑1 (MQ=255)

                                                                                                                                                                                                                                     gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGc                                                                                                                                                                         <  8:549966/148‑1 (MQ=255)

                                                                                                                                                                                                                                         tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtca                                                                                                                                                                    <  5:150237/149‑1 (MQ=255)

                                                                                                                                                                                                                                          cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACACATGGATAAGAGTGTGGACATCAGCGTcc                                     >  4:332817/1‑275 (MQ=255)

                                                                                                                                                                                                                                             ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactca                                                                                                                                                                <  6:184435/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTTACCGTGCAGACGGTTGCCGTCACTCACTAAAt                                                                                                                                                          >  8:633151/1‑148 (MQ=255)

                                                                                                                                                                                                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTAt                                                                                                                                                       >  6:59006/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAACATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATg                                                                                                                                                      <  7:736757/148‑1 (MQ=255)

                                                                                                                                                                                                                                                         cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                                                                                                                                                     >  8:516171/1‑148 (MQ=255)

                                                                                                                                                                                                                                                             cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGc                                                                       >  4:442756/1‑222 (MQ=255)

                                                                                                                                                                                                                                                              tGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATg                                                                                                                                                      >  5:639827/1‑142 (MQ=255)

                                                                                                                                                                                                                                                              tGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATg                                                                                                                                                      <  6:639827/142‑1 (MQ=255)

                                                                                                                                                                                                                                                                tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGt                                                                                                                                             <  5:296898/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCCCTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTtggt               >  2:698121/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTg                                                                                                                                            <  5:739142/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCAc                                                                   <  2:442756/222‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc                                                                                                                                                                                          <  7:919028/103‑1 (MQ=38)

                                                                                                                                                                                                                                                                 aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGt                                                                                                                                             >  7:956162/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                   aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTggtggt            <  4:636544/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                    gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa                                                                                                                                          >  8:393879/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                     ccGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaac                                                                                                                                         >  8:470844/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                      cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaaca                                                                                                                                        <  7:666378/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                        gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactcac                                                                                                                                                               >  4:594048/1‑123 (MQ=255)

                                                                                                                                                                                                                                                                        gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACNGTGAACAACAGGCTGGCTCACTGGCAAGAACCCGATGATAAAATTAAGCTGGCAATCAGAATAAAAATCAGCCCNCTGATGGAAAAGAGTGTGGACATAAGCGTCNNAGATAAACGAGCTTCTTTGGTGGTTAGGCCana   >  3:91921/1‑279 (MQ=255)

                                                                                                                                                                                                                                                                        gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAg                                                                                                                                              <  8:24034/140‑1 (MQ=255)

                                                                                                                                                                                                                                                                        gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGccc                                                                     >  2:195399/1‑213 (MQ=255)

                                                                                                                                                                                                                                                                            cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTaaa                                                                                                                                                           <  2:594048/123‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGtt                                                                                                                                 <  5:518106/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTg                                                                 <  4:195399/213‑1 (MQ=255)

                                                                                                                                                                                                                                                                              tG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTgg                                                                                                                               <  5:440605/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                              tG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTg                                                                                                                                >  7:583052/1‑148 (MQ=255)

GTCTGGCTTTTGCCTTCCTGATAGTACAGCCATGTTGCATGAAGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTGGTGGTCAGGCCAAAA  >  NC_000913/4520233‑4520776

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