BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,561,205	A→T	6.0%	intergenic (+307/‑163)	yffL → / → yffM	CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,561,205	0	A	T	6.0%	766.3 / 27.5	219	intergenic (+307/‑163)	yffL/yffM	CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (105/100); new base T (8/5); total (113/106)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 220 total aligned reads displayed.

TGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTC  >  NC_000913/2560946‑2561462

tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATCGGTCGGCTGGTTGCGCGCCTTACCCATGCTGATTAGCGCAATCATTTCAATTAACCTGGAAACGctcc                                                                                                                                                                                                                                                    >  1:35354/1‑272 (MQ=255)

   atttcggTGGGTTGAAAGTTGGCGGTAGTGGTGGGTTAAAAGTTGGCGTTAAGGTGGGTCAAAAGTTGCCGGTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTTCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTTCAGCAAGCTGGAAACGCCCAtt                                                                                                                                                                                                                                                  <  4:151031/269‑1 (MQ=255)

        ggTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGTTCGGTCATAGAGAGTCGCCTAGAGTTAATCAGATGAATTTTTTCGCCGCCTTGTAGCTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtaactca                                                                                                                                                                                                                                            >  2:19605/1‑275 (MQ=255)

                                                                                                                 ttAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                                                                                                                                 >  6:973409/1‑149 (MQ=255)

                                                                                                                   aaGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                               >  7:80319/1‑149 (MQ=255)

                                                                                                                    aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                                                                                                                                              <  8:80319/149‑1 (MQ=255)

                                                                                                                        ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaa                                                                                                                                                                                                                                                          <  7:976794/149‑1 (MQ=255)

                                                                                                                        ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGaaa                                                                                                                                                                                                                                                          <  5:439717/149‑1 (MQ=255)

                                                                                                                        ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactc                                                                                                                                                                                                                                             >  2:369824/1‑162 (MQ=255)

                                                                                                                           tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCCATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCAAGGCATTCGTACATGGGGAGCGCGGTATATCTCAAAGGGTGCTTTCTTATTGGTGGGCCTTCGGAAGGCTTTGAGCAATTg                                                                                                                         >  4:290975/1‑275 (MQ=255)

                                                                                                                            aGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGc                                                                                                                                                                                                                                         <  4:369823/162‑1 (MQ=255)

                                                                                                                            aGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                                                                                                                                      <  2:338551/181‑1 (MQ=255)

                                                                                                                             ggCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGccc                                                                                                                                                                                                                                                     <  5:1260639/149‑1 (MQ=255)

                                                                                                                              gCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCa                                                                                                                                                                                                                                                    >  5:955407/1‑149 (MQ=255)

                                                                                                                              gCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCa                                                                                                                                                                                                                                                    >  5:564365/1‑149 (MQ=255)

                                                                                                                                gACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCAtt                                                                                                                                                                                                                                                  >  5:774989/1‑149 (MQ=255)

                                                                                                                                   gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                               <  6:896809/149‑1 (MQ=255)

                                                                                                                                   gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                               >  6:1354734/1‑149 (MQ=255)

                                                                                                                                    ggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcact                                                                                                                                                                                                                                              >  8:246585/1‑149 (MQ=255)

                                                                                                                                    ggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcact                                                                                                                                                                                                                                              >  6:359678/1‑149 (MQ=255)

                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                                                                                                                                            <  6:622325/149‑1 (MQ=255)

                                                                                                                                        ttttttAGCCTGCCGGTAGGTCATAGAGAGTCGCCTAGAGTGAAACAGAAGAGTTATTTCGCCGCCTGAGAGATCGGCTGCTTGCGAGCATTAACCATGCTGATTAGCGCATTCATTTCCAGCAAGCTGGAAAAGCCCATTCACTCACg                                                                                                                                                                                                                                          >  8:1418291/1‑149 (MQ=255)

                                                                                                                                          ttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACaaa                                                                                                                                                                                                                >  2:305484/1‑173 (MQ=255)

                                                                                                                                              aGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcat                                                                                                                                                                                                            <  4:305484/173‑1 (MQ=255)

                                                                                                                                                ccTCCCGGGCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCt                                                                                                                                                                                                                                  <  5:1430614/149‑1 (MQ=255)

                                                                                                                                                  tCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatga                                                                                                                                                                                                                                                       >  6:477752/1‑123 (MQ=255)

                                                                                                                                                  tCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatga                                                                                                                                                                                                                                                       <  5:477752/126‑4 (MQ=255)

                                                                                                                                                   cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGaa                                                                                                                                                                                                                               <  6:794191/149‑1 (MQ=255)

                                                                                                                                                   cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGaa                                                                                                                                                                                                                               >  5:690792/1‑149 (MQ=255)

                                                                                                                                                   cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGaa                                                                                                                                                                                                                               <  5:704376/149‑1 (MQ=255)

                                                                                                                                                        tcggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                                                                                                                                          <  5:1422826/149‑1 (MQ=255)

                                                                                                                                                              aTAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATTCTGATTAGCGCAATCATTTACAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaa                                                                                                                                                                                                                    >  7:137486/1‑149 (MQ=255)

                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAc                                                                                                                                                                                                                              >  2:154995/1‑131 (MQ=255)

                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcat                                                                                                                                                                                                            >  8:5662/1‑149 (MQ=255)

                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATgg                                                                                                                                                                                                       >  3:23384/1‑154 (MQ=255)

                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCCGTGATGGTTATCCGCAAca                                                                              >  2:124948/1‑274 (MQ=255)

                                                                                                                                                                        gCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                                                                                                                                          >  6:988397/1‑149 (MQ=255)

                                                                                                                                                                         ccTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                                                                                                                                         >  8:1028181/1‑149 (MQ=255)

                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATTCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                                                                                                                                          <  4:154995/131‑1 (MQ=255)

                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAtt                                                                                                                                                                                                   <  1:23384/154‑1 (MQ=255)

                                                                                                                                                                           tAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATgg                                                                                                                                                                                                       >  7:737848/1‑149 (MQ=255)

                                                                                                                                                                             gagTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                               <  6:206803/91‑1 (MQ=255)

                                                                                                                                                                                   aaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGtt                                                                                                                                                                                               >  5:375024/1‑149 (MQ=255)

                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGTGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                            >  8:1216765/1‑149 (MQ=255)

                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                            >  5:28616/1‑149 (MQ=255)

                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                            >  6:1478576/1‑149 (MQ=255)

                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                            >  6:250602/1‑149 (MQ=255)

                                                                                                                                                                                       gaagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATg                                                                                                                                                                                           <  6:693739/149‑1 (MQ=255)

                                                                                                                                                                                       gaaaaaTGATTTCGCCGCCTGATAGGTCGGCTGGTTGCGGGCATTAACCATGCTGATTCGCGCAAGCATTTCCAGCAAGCTTGAAACGCCCATTCACTCACGCACAACCTCAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAATGGTGATGTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaa                                                                                                                <  4:145913/224‑1 (MQ=255)

                                                                                                                                                                                         agaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATggg                                                                                                                                                                                         >  5:197300/1‑149 (MQ=255)

                                                                                                                                                                                           aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                                                                                                                                       >  6:1035238/1‑149 (MQ=255)

                                                                                                                                                                                           aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                                                                                                                                       >  6:953047/1‑149 (MQ=255)

                                                                                                                                                                                           aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                                                                                                                                       <  6:426362/149‑1 (MQ=255)

                                                                                                                                                                                           aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                                                                                                                                       >  7:65448/1‑149 (MQ=255)

                                                                                                                                                                                             ttattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgc                                                                                                                                                                                     >  7:361091/1‑149 (MQ=255)

                                                                                                                                                                                                tttCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGcgcg                                                                                                                                                                                  <  7:397987/149‑1 (MQ=255)

                                                                                                                                                                                                  tCGCCGCCTGATAGGTCGGGTGGTTGCGAGCATTAACCATCCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCGCTCACCCAAAACCTCAACAGCAGACCAACAAACCACCATGGCATTAGTTCATGGGGAGCGCGGt                                                                                                                                                                                <  6:13746/149‑1 (MQ=255)

                                                                                                                                                                                                     ccgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGtata                                                                                                                                                                             <  7:921964/149‑1 (MQ=255)

                                                                                                                                                                                                      cgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTatat                                                                                                                                                                            <  8:141398/149‑1 (MQ=255)

                                                                                                                                                                                                           gATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaa                                                                                                                                                                       >  7:194458/1‑149 (MQ=255)

                                                                                                                                                                                                            aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgg                                                                                                                                                                                          >  7:617347/1‑129 (MQ=255)

                                                                                                                                                                                                            aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                                                                                                                                      <  8:395130/149‑1 (MQ=255)

                                                                                                                                                                                                             tAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAg                                                                                                                                                                     <  7:603820/149‑1 (MQ=255)

                                                                                                                                                                                                              aGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAgg                                                                                                                                                                    >  8:1076058/1‑149 (MQ=255)

                                                                                                                                                                                                                gTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                                                                                                                                  <  8:375381/149‑1 (MQ=255)

                                                                                                                                                                                                                 tCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtct                                                                                                                                                                          >  5:1415994/1‑140 (MQ=255)

                                                                                                                                                                                                                  cGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGa                                                                                                                                                                <  8:1326163/149‑1 (MQ=255)

                                                                                                                                                                                                                   ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACaaa                                                                                                                                                                                                                >  7:1032997/1‑100 (MQ=255)

                                                                                                                                                                                                                   ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATCTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAt                                                                                                                                                               >  8:480069/1‑149 (MQ=255)

                                                                                                                                                                                                                     cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                            <  6:549520/118‑1 (MQ=255)

                                                                                                                                                                                                                     cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAttt                                                                                                                                                             >  6:1104011/1‑149 (MQ=255)

                                                                                                                                                                                                                       ggTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCt                                                                                                                                                           <  7:538060/149‑1 (MQ=255)

                                                                                                                                                                                                                         ttGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                                                                                                                                         <  6:834914/149‑1 (MQ=255)

                                                                                                                                                                                                                          tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAt                                                                                                                                                        >  5:662681/1‑149 (MQ=255)

                                                                                                                                                                                                                           gCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt                                                                                                                                                       <  5:1182003/149‑1 (MQ=255)

                                                                                                                                                                                                                           gCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt                                                                                                                                                       >  6:1395863/1‑149 (MQ=255)

                                                                                                                                                                                                                             gAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                                                                                                                                     >  8:639066/1‑149 (MQ=255)

                                                                                                                                                                                                                             gAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                                                                                                                                     >  8:1155021/1‑149 (MQ=255)

                                                                                                                                                                                                                               gCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtg                                                                                                                                                   <  6:197300/149‑1 (MQ=255)

                                                                                                                                                                                                                                cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                                                                                                                                       >  6:1236498/1‑112 (MQ=255)

                                                                                                                                                                                                                                cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                                                                                                                                  <  6:712844/149‑1 (MQ=255)

                                                                                                                                                                                                                                     aCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCtt                                                                                                                                             <  5:890650/149‑1 (MQ=255)

                                                                                                                                                                                                                                          gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGaa                                                                                                                                        <  6:174157/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCtt                                                                                                                                   <  5:1376633/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCtt                                                                                                                                   >  7:1323779/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 aGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTc                                                                                                                                  <  5:146655/148‑1 (MQ=255)

                                                                                                                                                                                                                                                  gcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATGGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGACACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTc  <  2:290975/275‑1 (MQ=255)

                                                                                                                                                                                                                                                      gataaTTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCag                                                                                                                            >  8:939397/5‑148 (MQ=255)

                                                                                                                                                                                                                                                       aTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAt                                                                                                                           <  5:1281979/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAt                                                                                                                           <  8:439097/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAtt                                                                                                                          >  7:133746/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        gcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAtt                                                                                                                          >  8:678936/4‑149 (MQ=255)

                                                                                                                                                                                                                                                          gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa                                                                                                                                >  6:1439382/2‑141 (MQ=255)

                                                                                                                                                                                                                                                          gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGc                                                                                                                        >  7:1444816/2‑149 (MQ=255)

                                                                                                                                                                                                                                                           tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                                                                                                                                         <  5:1143904/115‑1 (MQ=255)

                                                                                                                                                                                                                                                           tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                                                                                                                                         >  6:1143904/1‑115 (MQ=255)

                                                                                                                                                                                                                                                           tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGcc                                                                                                                       >  6:510618/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               cAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCt                                                                                                                   >  8:895816/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               cAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCt                                                                                                                   <  7:1216765/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTa                                                                                                                 >  7:1098330/1‑149 (MQ=255)

TGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTC  >  NC_000913/2560946‑2561462

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