BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,790,300	A→T	5.9%	intergenic (‑219/+20)	waaH ← / ← tdh	LPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,790,300	0	A	T	5.9%	2537.1 / 98.1	704	intergenic (‑219/+20)	waaH/tdh	LPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand): ref base A (348/314); new base T (21/21); total (369/335)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-02

Only 100 of 712 total aligned reads displayed.

ATATAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCA  >  NC_000913/3790037‑3790564

                                                                     aTCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTAGTGTTAATCCCAGCTCAGAATAACTTTCACGGCCTGGCCCGAACGAATAGCGTcctgtatctt                                                                                                                                                                                           >  1:74092/1‑265 (MQ=255)

                                                                                                        aTTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGAt                                                                                                                                                       <  4:253618/275‑1 (MQ=255)

                                                                                                                         ttattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                                                                                                                                                                                    >  6:794938/1‑149 (MQ=255)

                                                                                                                         ttattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                                                                                                                                                                                    <  8:333748/149‑1 (MQ=255)

                                                                                                                         ttattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                                                                                                                                                                                    <  7:245469/149‑1 (MQ=255)

                                                                                                                           attTTTACAGGTGAGAGAGTCACTTCAGGGTGGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                                                                                                                                                                                  >  6:727174/1‑149 (MQ=255)

                                                                                                                               ttAAAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                                                                                                                                                                                              <  8:1023806/149‑1 (MQ=255)

                                                                                                                                 aCAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                                                                                                                                                                                            <  5:692451/149‑1 (MQ=255)

                                                                                                                                       gagagaGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                                                                                                                                      <  6:522238/149‑1 (MQ=255)

                                                                                                                                          agagTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                                                                                                                                                                                   <  6:1056876/149‑1 (MQ=255)

                                                                                                                                          agagTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                                                                                                                                                                                   <  7:335179/149‑1 (MQ=255)

                                                                                                                                            agTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                                                                                                                                                                                                 >  8:375018/1‑149 (MQ=255)

                                                                                                                                             gTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                <  8:213067/149‑1 (MQ=255)

                                                                                                                                              tCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                                                                                                                                                                                               <  8:1422523/149‑1 (MQ=255)

                                                                                                                                               cACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCTAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCa                                                                                                                 >  2:153170/1‑274 (MQ=255)

                                                                                                                                                aCTTCAGGGTA‑AGAGAATAATGAGAGCAGCTAAATACAAGATCAGCACAATTACGAAACCTTTATTACACGGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATc                                                                                                               <  4:280903/275‑1 (MQ=255)

                                                                                                                                                  ttCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                                                                                                                                                                                           >  5:1129835/1‑149 (MQ=255)

                                                                                                                                                    cAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                                                                                                                                                                                         <  5:761969/149‑1 (MQ=255)

                                                                                                                                                       ggTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                                                                                                                                                                                      <  6:941878/149‑1 (MQ=255)

                                                                                                                                                          agagag‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                   <  7:1085022/149‑1 (MQ=255)

                                                                                                                                                          agagag‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                   >  8:710332/1‑149 (MQ=255)

                                                                                                                                                           gagag‑aTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                                                                                                                                  <  7:1154520/149‑1 (MQ=255)

                                                                                                                                                              ag‑aTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                                                                                                                                                                                               >  8:1284971/1‑149 (MQ=255)

                                                                                                                                                                   aaTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                                                                                                                                  >  8:588479/1‑142 (MQ=255)

                                                                                                                                                                    aTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                                                                                                                                                                                          >  7:645365/1‑149 (MQ=255)

                                                                                                                                                                         agcagcTATTTAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATACTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGGCTgg                                                                                                                                                                                                                     >  6:30745/1‑149 (MQ=255)

                                                                                                                                                                         agcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                                                                                                                                                                                     >  8:1345023/1‑149 (MQ=255)

                                                                                                                                                                          gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGc                                                                                                                                                                                                                    <  6:123536/149‑1 (MQ=255)

                                                                                                                                                                           cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                                                                                                                                                                                   <  7:479194/149‑1 (MQ=255)

                                                                                                                                                                            agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                                                                                                                                                                                                  <  5:570439/149‑1 (MQ=255)

                                                                                                                                                                            agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                                                                                                                                                                                                  >  5:186424/1‑149 (MQ=255)

                                                                                                                                                                            agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                                                                                                                                                                                                  >  7:1247179/1‑149 (MQ=255)

                                                                                                                                                                             gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                                                                                                                                                                                                 <  8:545874/149‑1 (MQ=255)

                                                                                                                                                                               tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGaa                                                                                                                                                                                                               <  8:1221656/149‑1 (MQ=255)

                                                                                                                                                                                 tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                                                                                                                                                                             >  5:822772/1‑149 (MQ=255)

                                                                                                                                                                                  ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGTAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGc                                                                                                                                                                                                            >  7:562179/1‑149 (MQ=255)

                                                                                                                                                                                   taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                                                                                                                                                                                           >  5:377393/1‑149 (MQ=255)

                                                                                                                                                                                   taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                                                                                                                                                                                           >  5:317294/1‑149 (MQ=255)

                                                                                                                                                                                   taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                                                                                                                                                                                           >  5:922963/1‑149 (MQ=255)

                                                                                                                                                                                   taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                                                                                                                                                                                           >  7:1295450/1‑149 (MQ=255)

                                                                                                                                                                                    aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCAt                                                                                                                                                                                                          <  7:1441163/149‑1 (MQ=255)

                                                                                                                                                                                     agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                                                                                                                                                                                         <  6:1401136/149‑1 (MQ=255)

                                                                                                                                                                                     agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                                                                                                                                                                                         >  5:1361470/1‑149 (MQ=255)

                                                                                                                                                                                     agagATCACCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                                                                                                                                                                                         >  7:1345901/1‑149 (MQ=255)

                                                                                                                                                                                       agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAGGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                                                                                                                                                                                       <  6:1029667/149‑1 (MQ=255)

                                                                                                                                                                                         aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgc                                                                                                                        >  1:33440/1‑226 (MQ=255)

                                                                                                                                                                                           cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                                                                                                                                                                       >  3:35619/1‑177 (MQ=255)

                                                                                                                                                                                            aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                                                                                                                                                                                  <  7:1414154/149‑1 (MQ=255)

                                                                                                                                                                                             gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                                                                                                                                                                                              >  5:610888/1‑88 (MQ=255)

                                                                                                                                                                                             gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGACCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                 >  7:964350/1‑149 (MQ=255)

                                                                                                                                                                                             gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCATCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                <  4:115228/102‑1 (MQ=255)

                                                                                                                                                                                               acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGc                                                                                                                                                                                               >  8:691965/1‑149 (MQ=255)

                                                                                                                                                                                               acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                                                                                                                                                                <  4:185694/164‑1 (MQ=255)

                                                                                                                                                                                                       gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                                                                                                                                                                                       >  8:770578/1‑149 (MQ=255)

                                                                                                                                                                                                       gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                                                                                                                                                                                       >  7:379620/1‑149 (MQ=255)

                                                                                                                                                                                                         aCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCTCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                                                                                                                                  >  6:1183578/1‑104 (MQ=255)

                                                                                                                                                                                                          ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGgaa                                                                                                                                                                                    >  5:1039968/1‑149 (MQ=255)

                                                                                                                                                                                                          ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                                                                                                                                                                    >  7:411093/1‑149 (MQ=255)

                                                                                                                                                                                                           tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                                                                                                                                                                   >  6:267304/1‑149 (MQ=255)

                                                                                                                                                                                                             tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca                                                                                                                                                                                 >  5:1318865/1‑149 (MQ=255)

                                                                                                                                                                                                             tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca                                                                                                                                                                                 >  6:151879/1‑149 (MQ=255)

                                                                                                                                                                                                              ttattaCACAGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                                                                                                                                                                <  7:856197/149‑1 (MQ=255)

                                                                                                                                                                                                               tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                                                                                                                                                                               <  8:1401459/149‑1 (MQ=255)

                                                                                                                                                                                                               tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                                                                                                                                                                               <  6:1129835/149‑1 (MQ=255)

                                                                                                                                                                                                                  taCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGAt                                                                                                                                                                            >  7:694/1‑149 (MQ=255)

                                                                                                                                                                                                                    cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCtt                                                                                                                                                                                                                                                              >  6:651835/1‑65 (MQ=255)

                                                                                                                                                                                                                     acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCg                                                                                                                                          >  2:74561/1‑180 (MQ=255)

                                                                                                                                                                                                                       tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCACACATATCACGACCGTAAATACCTTTAATGAACAAGCATTTAAAAATCACTTTgt                                         >  1:11064/1‑274 (MQ=255)

                                                                                                                                                                                                                        gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaa                                                                                                                                                                      >  8:429091/1‑149 (MQ=255)

                                                                                                                                                                                                                         aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa                                                                                                                                                                     >  6:851355/1‑149 (MQ=255)

                                                                                                                                                                                                                           aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACg                                                                                                                                                                   >  5:472732/1‑149 (MQ=255)

                                                                                                                                                                                                                                       atTTAACGGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGAt                                                                                                                                                       <  2:27730/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGAt                                                                                                                                                       >  5:1441257/1‑148 (MQ=255)

                                                                                                                                                                                                                                         ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                                                                                                                                                                                    <  5:291200/102‑1 (MQ=255)

                                                                                                                                                                                                                                          tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                                                                                                                                                                               >  6:922068/1‑122 (MQ=255)

                                                                                                                                                                                                                                           aaCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGc                                                                                                                                                   <  6:1331954/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aaCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGc                                                                                                                                                   <  5:96210/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  tATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGAt                                                                                                                                            <  8:1155643/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATc                                                                                                                                           <  6:822772/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATc                                                                                                                                           <  6:378123/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    tCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCg                                                                                                                                          >  8:283226/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACg                                                                                                                                                                   <  2:374532/123‑1 (MQ=255)

                                                                                                                                                                                                                                                      cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                        >  7:176951/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                        <  8:1112107/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       tCAGATAAAAGGGCTGGAATACTAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTGTCCCGGACTGGTCCGAACGCATAGCGTCAAAGCCCTTCTGTCAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATGGAgg                                                                                                                                       >  8:908953/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                                                                                                                                       >  6:473761/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         aGATAAAAGGGCTGGTATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGGGCCAAAGCCCTTCTGGAACTCATCGAGAGAGACAAGATGGGTAATGATCGGCGAAAGATCGAGGCAAGACTGAAACAGCGCCGCCAACTTGTACCCGGTTTCAAACATCTCACGACCGTAAATACATTTAATGAGCAAGCCTTTAAAGAACACTTTGGTCCAGTCACAAGACTAATCCGACGGCGGAATAc        >  4:64500/1‑274 (MQ=255)

                                                                                                                                                                                                                                                         aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGcc                                                                                                                                     >  8:865413/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                                                                                                                                                                                       <  5:1218227/98‑1 (MQ=255)

                                                                                                                                                                                                                                                          gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGAACGGCGACAGATCGAGGCCAGACTGAATCAgcg                                                                                                                       >  2:264802/1‑162 (MQ=255)

                                                                                                                                                                                                                                                          cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaaa                                                                                                                                               <  7:1126208/133‑1 (MQ=255)

                                                                                                                                                                                                                                                          cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaaa                                                                                                                                               >  8:1126208/6‑138 (MQ=255)

                                                                                                                                                                                                                                                            tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc                                                                                                                                                                    <  5:520176/115‑1 (MQ=255)

                                                                                                                                                                                                                                                            tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGa                                                                                                                                  <  5:851355/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                                                                                                                                                                                  >  6:596677/2‑99 (MQ=255)

                                                                                                                                                                                                                                                                aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGAt                                                                                                                                                                 >  6:972445/1‑114 (MQ=255)

                                                                                                                                                                                                                                                                 gggCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGaa                                                                                                                             >  7:75077/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAACGATCGGCGAAAGATCGAGGCCAGACTGAATCa                                                                                                                          >  8:150320/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCNAAAGATCGAGGCCAGACTGAATCNGCGCNGCCATCtt                                                                                                             >  4:348049/1‑161 (MQ=255)

                                                                                                                                                                                                                                                                       gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                             <  2:191016/95‑1 (MQ=255)

ATATAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAG‑ATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCA  >  NC_000913/3790037‑3790564

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