BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,520,500	Δ1 bp	8.6%	intergenic (‑176/+171)	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,500	0	A	.	8.6%	603.4 / 45.8	174	intergenic (‑176/+171)	yjhU/yjhF	putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand): ref base A (81/78); new base . (7/8); total (88/86)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 181 total aligned reads displayed.

GTCTGGCTTTTGCCTTCCTGATAGTACAGCCATGTTGCATGAAGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTGGTGGTCAGGCCAAA  >  NC_000913/4520233‑4520775

                                      aTGAAGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                                                                                                                                                                                                   >  2:372552/1‑233 (MQ=255)

                                          aGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                              <  4:372552/233‑1 (MQ=255)

                                                                                                                             cATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑c                                                                                                                                                                                                                                                                                >  7:921901/1‑148 (MQ=255)

                                                                                                                             cATCTTAACCCTGTTTACCCGGAGAGCTGATGAAGATTGCCATTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCTGCATACCATATCAGCGATAGCCGGCTGTTCTCCCATACTGTAAAGACGGAGACATG‑c                                                                                                                                                                                                                                                                                >  8:128149/1‑148 (MQ=255)

                                                                                                                              aTCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTc                                                                                                                                                                                                                                                                              >  6:266525/1‑149 (MQ=255)

                                                                                                                                  tAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCgg                                                                                                                                                                                                                                                                           >  7:671885/1‑148 (MQ=255)

                                                                                                                                  tAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCgg                                                                                                                                                                                                                                                                           <  8:408145/148‑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)

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

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

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

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

                                                                                                                                      ccTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                                                                                                                                                                                                         <  5:618033/146‑1 (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)

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

                                                                                                                                               ccGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                                                                                                                                                                                                             >  5:758219/1‑149 (MQ=255)

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

                                                                                                                                                 ggAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTa                                                                                                                                                                                                                                                                      >  8:634441/1‑138 (MQ=255)

                                                                                                                                                         gATGCAGATTGCCAGTCGCATGATTTTATACAGATATCTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCCATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                                                                                                                                                                                                   <  6:862471/149‑1 (MQ=255)

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

                                                                                                                                                                 ttGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                                                                                                                                                                                            >  7:483943/1‑148 (MQ=255)

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

                                                                                                                                                                    ccAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                                                                                                                                                                                                        <  6:402411/149‑1 (MQ=255)

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

                                                                                                                                                                            aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                                                                                                                                                                                                 <  8:18638/148‑1 (MQ=255)

                                                                                                                                                                              gATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                                                                                                                                                                                                              >  5:567774/1‑149 (MQ=255)

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

                                                                                                                                                                                     taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                                                                                                                                                                                                       <  6:87583/149‑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)

                                                                                                                                                                                              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)

                                                                                                                                                                                                  ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                                                                                                                                                                                                 >  8:274611/1‑126 (MQ=255)

                                                                                                                                                                                                  ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                                                                                                                                                                                                 <  7:274611/126‑1 (MQ=255)

                                                                                                                                                                                                   tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                                                                                                                                                                                                         <  5:182048/149‑1 (MQ=255)

                                                                                                                                                                                                     ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                                                                                                                                                                                                       <  6:59353/149‑1 (MQ=255)

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

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

                                                                                                                                                                                                           aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                                                                                                                                                                                                 <  6:519289/149‑1 (MQ=255)

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

                                                                                                                                                                                                             gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                                                                                                                                                                                                        >  7:995483/1‑140 (MQ=255)

                                                                                                                                                                                                               aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTAt                                                                                                                                                                                                                                                                >  4:252699/1‑82 (MQ=255)

                                                                                                                                                                                                               aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCTCCCCTTTTTGATCGCCAGCGAtt                                                                                                                                                                                             <  6:686682/149‑1 (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:636447/149‑1 (MQ=255)

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

                                                                                                                                                                                                                aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCAATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTAtg                                                                                                                                                                                                                                                               >  8:477776/1‑81 (MQ=21)

                                                                                                                                                                                                                aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCAATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTAtg                                                                                                                                                                                                                                                               <  7:477776/82‑2 (MQ=21)

                                                                                                                                                                                                                   aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc                                                                                                                                                                                         <  6:909045/149‑1 (MQ=255)

                                                                                                                                                                                                                    aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                                                                                                                                                                                                       <  5:136978/134‑1 (MQ=255)

                                                                                                                                                                                                                        gCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                                                                                                                                                                                                          <  7:384823/127‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATaa                                                                              >  4:519166/1‑250 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATaaaga                                                     >  4:448845/1‑272 (MQ=255)

                                                                                                                                                                                                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAagag                                                     >  4:620740/1‑275 (MQ=255)

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

                                                                                                                                                                                                                             ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTa                                                                                                                                                                                                                                                              >  8:388984/1‑70 (MQ=255)

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

                                                                                                                                                                                                                             ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAAt                                                                          <  2:519166/250‑1 (MQ=255)

                                                                                                                                                                                                                              cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGa                                                                                                                                                                               >  8:736757/1‑148 (MQ=255)

                                                                                                                                                                                                                                tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtcactc                                                                                                                                                                >  2:540780/1‑161 (MQ=255)

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

                                                                                                                                                                                                                                 ttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATaaa                                                                                                    >  2:660843/1‑220 (MQ=255)

                                                                                                                                                                                                                                    aGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTa                                                                                                                                                            <  4:540780/161‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                     gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTa                                                                                                <  4:660843/220‑1 (MQ=255)

                                                                                                                                                                                                                                        aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATaa                                                                                                                                                                                                                                     >  2:772657/1‑84 (MQ=255)

                                                                                                                                                                                                                                          cGCCCGCTGTTCCCCCATACTGTAAATCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAACGAAGGTTCAATATCGCCCCTTGTGGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAg                                                                                      <  2:711515/225‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                                     cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAAtt                                                                                                                                                        <  8:233120/148‑1 (MQ=255)

                                                                                                                                                                                                                                                        ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                                                                                                                                                    >  5:1051576/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                          aTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATggg                                                                                                                                                   <  8:483943/148‑1 (MQ=255)

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

                                                                                                                                                                                                                                                               gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACa                                                                                                                                              >  7:732663/1‑148 (MQ=255)

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

                                                                                                                                                                                                                                                                 aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGt                                                                                                                                            >  7:581232/1‑148 (MQ=255)

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

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

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

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

                                                                                                                                                                                                                                                                        gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAg                                                                                                                                             >  7:24034/1‑140 (MQ=255)

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

GTCTGGCTTTTGCCTTCCTGATAGTACAGCCATGTTGCATGAAGTACGGGATCATCATTGAACAGGCTGCTGTTTGTTTGATCTCTATCCACTTTTGTATACCTGTATTAAACCGGGGAAATGCTCATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGGTTGGCTCACTGGCAAGAACCCGATTATAAAACTAAGCTGGCGAGCAGAATAAAAATCAGCCCACTGATGGATAAGAGTGTGGACATCAGCGTCCAGGATAACAGCGTTTCTTTGGTGGTCAGGCCAAA  >  NC_000913/4520233‑4520775

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