BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,418,629	T→C	6.1%	intergenic (+18/‑413)	yhdV → / → yhdW	putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,418,629	0	T	C	6.1%	527.6 / 16.6	165	intergenic (+18/‑413)	yhdV/yhdW	putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand): ref base T (84/70); new base C (3/7); total (88/77)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-02

Only 100 of 168 total aligned reads displayed.

CGTTAGCGTTTATTACTGAGGTAACACCATGAAAAGATTAATTCCTGTGGCATTGCTCACCGCATTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCA  >  NC_000913/3418362‑3418871

ggttggggttTAGTACTGAGGTAACACCATGAAAAGATTAATGCCGTTGGCACGGATAGCCGGATTGTTCGCGCGGTTTGCTCACGATTCGGCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACGCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgcg                                                                                                                                                                                                                                              <  4:309122/267‑1 (MQ=255)

          ttttaCTGAGGTAAAACCCTGAAAAGATTAATGCCTGTGGCATTGCTGCCCGCATTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAAc                                                                                                                                                                                                                                   <  4:402199/273‑1 (MQ=255)

                              gAAAAGATTAATTCCTGTGGCATTGCTCACCGCATTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCagtc                                                                                                                                                                                                               >  2:123196/1‑272 (MQ=255)

                                     ttcttcccTGTGGCATTGCTCCCCGCAGTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATTAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                                                                                                                        <  4:236901/269‑1 (MQ=255)

                                     ttAATTCCTGTGGCATTGCTCACCGCATTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTTCTGTTCCTGGCCTGACGCTGGGGATCATCGCTCTATCGGAATAACTGAACAGGGCCCGGAGTTCCGCGCTCTTTTTAACGATACACTATCCCCCCCtaaaataaaa                                                                                                                                                                                                        >  4:12752/1‑275 (MQ=255)

                                                      gCTCACCGCATTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCCCTATCGAAATAACAGAACAGGGCGCGGAGTTCCGCGCGCTTTTTAACGATACCGTATCTCTCCATAAGAAAAAATATTGCAATATTctgct                                                                                                                                                                                       >  4:198974/1‑275 (MQ=255)

                                                                                      aTTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACAGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACa                                                                                                                                                        <  2:198974/274‑1 (MQ=255)

                                                                                                                        tcgtcTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGtt                                                                                                                                                                                                                                                   >  6:703621/1‑149 (MQ=255)

                                                                                                                           tcTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTc                                                                                                                                                                              >  2:46415/1‑215 (MQ=255)

                                                                                                                            cTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgc                                                                                                                                                                                                                                               <  5:1009364/149‑1 (MQ=255)

                                                                                                                             tGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTc                                                                                                                                                                                                                                                  >  2:340787/1‑145 (MQ=255)

                                                                                                                             tGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGGTGcgcg                                                                                                                                                                                                                                              >  7:817025/1‑149 (MQ=255)

                                                                                                                               gTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctg                                                                                                                                                                          <  4:46415/215‑1 (MQ=255)

                                                                                                                                ttCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctc                                                                                                                                                                                                                                           >  8:237648/1‑149 (MQ=255)

                                                                                                                                ttCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctc                                                                                                                                                                                                                                           >  7:851436/1‑149 (MQ=255)

                                                                                                                                ttCACACCAATACCAGCATGCAAGGGACGACACAGGATAACAGGGAAACGCCTTGCGCTATTGCCGGCGGGGCAGCAACTGTTGCTGGCCTGTCGATGGGGATGATCGCGCTACCGAAATAACTGAACATTGCCCGGGGTTCCGCGctc                                                                                                                                                                                                                                           >  8:191239/1‑149 (MQ=255)

                                                                                                                                 tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgcg                                                                                                                                                                                                                                              <  4:340787/145‑1 (MQ=255)

                                                                                                                                 tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                                                                                                                          <  8:657115/149‑1 (MQ=255)

                                                                                                                                 tCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccct                                                                                                                                                                                                                                          <  8:102843/149‑1 (MQ=255)

                                                                                                                                   acacCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttt                                                                                                                                                                                                                                        >  5:1063739/1‑149 (MQ=255)

                                                                                                                                    cacCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCtttt                                                                                                                                                                                                                                       <  5:653632/149‑1 (MQ=255)

                                                                                                                                     acCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGTGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                                                                                                                      >  5:760989/1‑149 (MQ=255)

                                                                                                                                     acCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctc                                                                                                                                                                                                                      >  1:190257/1‑165 (MQ=255)

                                                                                                                                         aTACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCa                                                                                                                                                                                                                  <  3:190257/165‑1 (MQ=255)

                                                                                                                                         aTACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATAATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACg                                                                                                                                                                                                                                  >  5:191131/1‑149 (MQ=255)

                                                                                                                                          tACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGa                                                                                                                                                                                                                                 <  5:556563/149‑1 (MQ=255)

                                                                                                                                          tACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACtt                                                                                                    >  2:393937/1‑274 (MQ=255)

                                                                                                                                           aCCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGAt                                                                                                                                                                                                                                <  6:1087417/149‑1 (MQ=255)

                                                                                                                                            ccTGTATGAAAGGCACGACACAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGGACAGGGCGCGGAGTTCCGCGCTCTTTTTAGCGATa                                                                                                                                                                                                                               <  6:961962/149‑1 (MQ=255)

                                                                                                                                               gTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAg                                                                                                                                                                                                                            >  8:445437/1‑149 (MQ=255)

                                                                                                                                                tATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTc                                                                                                                                                                                                                                                  >  4:691077/1‑126 (MQ=255)

                                                                                                                                                   gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATa                                                                                                                                                                                                                               <  4:737517/142‑1 (MQ=255)

                                                                                                                                                   gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                                                                                                                        >  7:918807/1‑149 (MQ=255)

                                                                                                                                                   gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                                                                                                                        >  6:809723/1‑149 (MQ=255)

                                                                                                                                                   gAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtc                                                                                                                                                                                                                        >  7:764017/1‑149 (MQ=255)

                                                                                                                                                      aGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctct                                                                                                                                                                                                                     >  7:1032830/1‑149 (MQ=255)

                                                                                                                                                      aGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctct                                                                                                                                                                                                                     <  7:239229/149‑1 (MQ=255)

                                                                                                                                                          aCGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaa                                                                                                                                                                                                          >  4:738498/1‑156 (MQ=255)

                                                                                                                                                              cTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGa                                                                                                                                                                                                             <  6:997470/149‑1 (MQ=255)

                                                                                                                                                              cTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAata                                                                                                                                                                                                      <  2:738498/156‑1 (MQ=255)

                                                                                                                                                               tCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGAt                                                                                                                                                                                                            >  6:545760/1‑149 (MQ=255)

                                                                                                                                                                     tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAata                                                                                                                                                                                                      >  8:706373/1‑149 (MQ=255)

                                                                                                                                                                     tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                                                                                                                          >  5:678940/1‑113 (MQ=255)

                                                                                                                                                                     tAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                                                                                                                                                                          <  6:678940/113‑1 (MQ=255)

                                                                                                                                                                      aaCTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAatat                                                                                                                                                                                                     >  7:648247/1‑149 (MQ=255)

                                                                                                                                                                        cTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTg                                                                                                                                                                                                   <  8:52729/149‑1 (MQ=255)

                                                                                                                                                                           ggAGACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATTGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCatta                                                                  <  1:75182/275‑1 (MQ=255)

                                                                                                                                                                           ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                                                                                                                        <  8:778407/141‑1 (MQ=255)

                                                                                                                                                                           ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                                                                                                                                >  7:432149/1‑149 (MQ=255)

                                                                                                                                                                           ggAAACGGCTGAGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                                                                                                                                >  7:614190/1‑149 (MQ=255)

                                                                                                                                                                            gAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                                                                                                                        >  8:302580/1‑140 (MQ=255)

                                                                                                                                                                             aaaCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTc                                                                                                                                                                                                                                                  >  4:678642/1‑97 (MQ=255)

                                                                                                                                                                               aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATAtt                                                                                                                                                                                            >  8:992185/1‑149 (MQ=255)

                                                                                                                                                                                 ggCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgcg                                                                                                                                                                                                                                              <  2:678642/97‑1 (MQ=255)

                                                                                                                                                                                   cTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgc                                                                                                                                                                                        <  7:566767/149‑1 (MQ=255)

                                                                                                                                                                                    tGGGGCTATTGCCGGCGGGGCAGCTACTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGCCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgct                                                                                                                                                                                       <  7:191239/149‑1 (MQ=255)

                                                                                                                                                                                      gggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTgg                                                                                                                                                                                     <  6:1035412/149‑1 (MQ=255)

                                                                                                                                                                                                cggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCct                                                                                                                                                                           <  5:1080959/149‑1 (MQ=255)

                                                                                                                                                                                                  gcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgc                                                                                                                                                                         <  5:746672/149‑1 (MQ=255)

                                                                                                                                                                                                    ggggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAGCTGAACAGGGCTCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                                                                                                                                                                       <  8:1034117/149‑1 (MQ=255)

                                                                                                                                                                                                       gCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAtt                                                                                                                                                                    >  7:1030135/1‑149 (MQ=255)

                                                                                                                                                                                                                gTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCtttt                                                                                                                                                                                                                                       >  4:468217/1‑73 (MQ=255)

                                                                                                                                                                                                                    cTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAAc                                                                                                                                                                                                                                   <  2:468217/73‑1 (MQ=255)

                                                                                                                                                                                                                     tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCAt                                                                                                                                                              >  5:196200/1‑141 (MQ=255)

                                                                                                                                                                                                                     tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCAt                                                                                                                                                              <  6:196200/141‑1 (MQ=255)

                                                                                                                                                                                                                     tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATc                                                                                                                                                      >  5:294249/1‑149 (MQ=255)

                                                                                                                                                                                                                         cTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt                                                                                                                                                  >  7:995224/1‑149 (MQ=255)

                                                                                                                                                                                                                           gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCAt                                                                                                                                                                                                                 >  4:721672/1‑84 (MQ=255)

                                                                                                                                                                                                                           gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTa                                                                                                                                                >  5:947793/1‑149 (MQ=255)

                                                                                                                                                                                                                            aCGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTAt                                                                                                                                               >  5:883163/1‑149 (MQ=255)

                                                                                                                                                                                                                               aTGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGa                                                                                                                                                                                                             <  2:721672/84‑1 (MQ=255)

                                                                                                                                                                                                                                  gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                                                                                                                         <  8:57911/149‑1 (MQ=255)

                                                                                                                                                                                                                                     atcatcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgct                                                                                                                                      >  7:678683/1‑149 (MQ=255)

                                                                                                                                                                                                                                      tcatcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCtt                                                                                                                                     >  6:48168/1‑149 (MQ=255)

                                                                                                                                                                                                                                        atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                                                                                                                                                                               <  6:315088/73‑1 (MQ=255)

                                                                                                                                                                                                                                        atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTcc                                                                                                                                                                             >  6:690160/1‑107 (MQ=255)

                                                                                                                                                                                                                                        atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTcc                                                                                                                                                                             <  5:690160/107‑1 (MQ=255)

                                                                                                                                                                                                                                        atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTc                                                                                                                                   >  7:12444/1‑149 (MQ=255)

                                                                                                                                                                                                                                           gCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCa  <  2:12752/275‑1 (MQ=255)

                                                                                                                                                                                                                                            ctctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATGGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGaaa                                                                                                                               <  6:760989/149‑1 (MQ=255)

                                                                                                                                                                                                                                             tctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAAt                                                                                                                              <  7:202276/149‑1 (MQ=255)

                                                                                                                                                                                                                                              ctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTaa                                                                           >  2:330706/1‑199 (MQ=255)

                                                                                                                                                                                                                                               tATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                                                                                                                            >  5:353246/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                                                                                                                            >  7:458926/1‑149 (MQ=255)

                                                                                                                                                                                                                                                aTCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCtata                                                                                                                                                                                              <  3:40983/82‑1 (MQ=255)

                                                                                                                                                                                                                                                aTCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAAc                                                                                                                           >  7:790730/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                                                                                                                          >  8:362121/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                                                                                                                          >  7:726837/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCAt                                                                                                                       >  8:494095/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      aTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCa                                                                                                                     >  5:228918/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatt                                                                                                                  <  8:842111/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              aaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaatt                                                                                                             >  5:518772/6‑149 (MQ=255)

                                                                                                                                                                                                                                                                agagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTg                                                                                                                                                                                                   <  8:808397/58‑1 (MQ=38)

                                                                                                                                                                                                                                                                aGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAAtttt                                                                                                           >  8:814033/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  ggCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAata                                                                                       >  1:118481/1‑167 (MQ=255)

                                                                                                                                                                                                                                                                  ggCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGACTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTca                                                                                                         <  7:623542/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    cgcgGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAAt                                                                                                                                                                                >  2:725101/1‑76 (MQ=255)

                                                                                                                                                                                                                                                                       ggAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACtt                                                                                                    >  8:556867/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           ttCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt                                                                                                                                                  >  6:971023/1‑99 (MQ=255)

                                                                                                                                                                                                                                                                           ttCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt                                                                                                                                                  <  5:971023/99‑1 (MQ=255)

CGTTAGCGTTTATTACTGAGGTAACACCATGAAAAGATTAATTCCTGTGGCATTGCTCACCGCATTGCTGGCGGGCTGTGCTCACGATTCGCCCTGCGTTCCGGTTTATGACGATCAGGGTCGTCTGGTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCA  >  NC_000913/3418362‑3418871

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