BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,561,206	A→T	5.9%	intergenic (+308/‑162)	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,206	0	A	T	5.9%	773.6 / 27.5	219	intergenic (+308/‑162)	yffL/yffM	CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (105/101); new base T (8/5); total (113/106)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 222 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)

                                                                                                                  tAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGc                                                                                                                                                                                                                                                                <  8:453322/149‑1 (MQ=255)

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

                                                                                                                    aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                                                                                                                                              >  5:693739/1‑149 (MQ=255)

                                                                                                                        ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaa                                                                                                                                                                                                                                                          <  8:184097/149‑1 (MQ=255)

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

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

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

                                                                                                                            aGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGcc                                                                                                                                                                                                                                                      >  6:1376633/1‑149 (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:564365/1‑149 (MQ=255)

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

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

                                                                                                                              gCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCa                                                                                                                                                                                                                                                    <  6:1465783/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                                                                                                                                            <  5:1289127/149‑1 (MQ=255)

                                                                                                                                         tttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGc                                                                                                                                                                                                                                         <  7:432350/149‑1 (MQ=255)

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

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

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

                                                                                                                                                    ccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAc                                                                                                                                                                                                                              >  7:1069321/1‑149 (MQ=255)

                                                                                                                                                        tcggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCCAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                                                                                                                                          <  8:1120045/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                        gCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                                                                                                                                          >  8:1176386/1‑149 (MQ=255)

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

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

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

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

                                                                                                                                                                                ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcact                                                                                                                                                                                                                                              >  6:136932/1‑105 (MQ=255)

                                                                                                                                                                                  aaaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGt                                                                                                                                                                                                >  8:397987/1‑149 (MQ=255)

                                                                                                                                                                                  aaaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCACTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGt                                                                                                                                                                                                >  6:1182003/1‑149 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                                                                       gaagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGACATTCGTTCATg                                                                                                                                                                                           <  5:1474879/149‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                            attattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAg                                                                                                                                                                                      >  7:1418048/1‑149 (MQ=255)

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

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

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

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

                                                                                                                                                                                                            aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgg                                                                                                                                                                                          <  8:617347/129‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                 tCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtct                                                                                                                                                                          <  6:1415994/140‑1 (MQ=255)

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

                                                                                                                                                                                                                   ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACaaa                                                                                                                                                                                                                <  8:1032997/100‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                                                                             gAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                                                                                                                                     >  8:639066/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)

                                                                                                                                                                                                                                cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                                                                                                                                  >  8:486348/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                        aTGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCgg                                                                                                                                          >  8:333281/1‑149 (MQ=255)

                                                                                                                                                                                                                                              aTTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCAACGAATGTACGTAATCTTGCATTGCGCAGAAACGCTGGGCACTGAGGGCAGTTTCCGAATTAGTTTTGGGGAcac      >  4:310825/1‑273 (MQ=255)

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

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

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

                                                                                                                                                                                                                                                  gcgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa                                                                                                                                <  7:55641/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                      aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCGTTCGGAAGGCTTCGAGCaa                                                                                                                            <  7:855930/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa                                                                                                                            >  5:763413/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                                                                          gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa                                                                                                                                <  5:1439382/140‑1 (MQ=255)

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

                                                                                                                                                                                                                                                          gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGATTCCGTTCTTGGGGGGCCCGGATAATTTCAAAGGGTAATTCTTTATTGTGGGCCCTTCGGAGGGCTCCGACCattccc                                                                                                                        >  6:1129176/2‑144 (MQ=255)

                                                                                                                                                                                                                                                          aTTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGc                                                                                                                        >  7:788553/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                                                                            ttCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa                                                                                                                      <  8:1418048/149‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                                                   aaGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaaa                                                                                                               <  5:506839/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    aGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaaaa                                                                                                              >  6:176561/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` `< 29 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`