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%	680.8 / 22.5	195	intergenic (+307/‑163)	yffL/yffM	CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (98/85); new base T (4/8); total (102/93)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

Only 100 of 197 total aligned reads displayed.

TGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTCCCTGCTTCTGGAGCAA  >  NC_000913/2560946‑2561478

tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGCAGAATTATTTCGCCGCCTGATAGGTCGTCTGGTTGCGAGCATTAACCATTCTGATTAGCGCAATAATTTCCACCAAGCTGGAAACgacaa                                                                                                                                                                                                                                                                    >  4:499898/1‑271 (MQ=255)

           aggTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGGTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGACTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGc                                                                                                                                                                                                                                                         <  2:424168/274‑1 (MQ=255)

                                                 cGTTAAGGTGTGTTAAAAGTTGCCGTTAAGGTGGGTTAAAAGTTGTCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAAGTATTTCGCCGCCTGATAGGTCAGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACACACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAtt                                                                                                                                                                                                                   <  4:602860/275‑1 (MQ=255)

                                                                                  ggTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                                                                                                                                                  >  2:409762/1‑275 (MQ=255)

                                                                                    tCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATTACTCAAAGGGGGa                                                                                                                                                                                >  4:411737/1‑275 (MQ=255)

                                                                                                               taTTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCa                                                                                                                                                                                                                                                                                   >  6:477513/1‑149 (MQ=255)

                                                                                                               taTTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCa                                                                                                                                                                                                                                                                                   >  5:528964/1‑149 (MQ=255)

                                                                                                                 ttAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACTATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACAAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGCGATTTACTATTGGTGCGCCTTCGGAAGGCtt                                                                                                                                                   >  1:155357/1‑275 (MQ=255)

                                                                                                                 ttAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                                                                                                                                                 <  7:379523/149‑1 (MQ=255)

                                                                                                                   aaGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               >  8:120525/1‑149 (MQ=255)

                                                                                                                       aTTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGaa                                                                                                                                                                                                                                                                           <  6:159821/149‑1 (MQ=255)

                                                                                                                        ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTAGaaa                                                                                                                                                                                                                                                                          <  8:526656/149‑1 (MQ=255)

                                                                                                                         tCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaat                                                                                                                                                                                                                                                                         <  5:483184/149‑2 (MQ=255)

                                                                                                                          cTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgggg                                                                                                                                                                                                        >  4:482197/1‑213 (MQ=255)

                                                                                                                               cGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatgattgc                                                                                                                                                                                                                                                                   <  6:864642/149‑8 (MQ=255)

                                                                                                                               cGACGGGAGTATTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                                                                                                                                                          <  2:542007/174‑1 (MQ=255)

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

                                                                                                                                gACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCAAGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAAGGGGAGCGCGGTATATCCCAAAGGGTGATTTCTTAGTGGTGGGCCTTGGGAAGGCTTCGAGCAATTGCCAGc                                                                                                                                    >  4:209799/1‑275 (MQ=255)

                                                                                                                                   gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                                               <  7:932792/149‑1 (MQ=255)

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

                                                                                                                                     gAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtccctc                                                                                                                                                                                                                                                             <  8:575972/149‑1 (MQ=255)

                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                                                                                                                                                            <  7:826176/149‑1 (MQ=255)

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

                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                                                                                                                                                    <  2:704169/157‑1 (MQ=255)

                                                                                                                                               gCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAcc                                                                                                                                                                                                                                                   <  7:393065/149‑1 (MQ=255)

                                                                                                                                                ccTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               <  8:314341/120‑1 (MQ=255)

                                                                                                                                                ccTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               >  7:314341/1‑120 (MQ=255)

                                                                                                                                                  tCCCGGTCGGTCATAGAGAGTCGCCTAAAGTTAAACAGCAGAATTATTTCGCCGCCTGATATGTCGGCTGGTTGAGAGCATTAACCATGCTGATTAGCGCCATCCTTTCCAGCAAGCTGTACAGACCCATTCACTCACGCCCAACCTGAACAGCAGACAAAAAAACCATAATGGAATACGTTCATGGGGACCGCGGTATATGTCAAAGGGTGATTTGTTATTGTTGGGCCTTCGGATTGCTTGNNNGAATGTCCAGATGAAAAGACGGGcatg                                                                                                                    >  3:190908/1‑270 (MQ=255)

                                                                                                                                                   cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               >  5:64158/1‑117 (MQ=255)

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

                                                                                                                                                         cggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagcag                                                                                                                                                                                                                                         <  7:449149/149‑1 (MQ=255)

                                                                                                                                                            tcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                                                                                                                                                      >  5:17892/1‑149 (MQ=255)

                                                                                                                                                            tcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                                                                                                                                                      <  6:434287/149‑1 (MQ=255)

                                                                                                                                                            tcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                                                                                                                                                      >  5:1092274/1‑149 (MQ=255)

                                                                                                                                                                agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACCCACAACCTg                                                                                                                                                                                                                                                 <  7:850235/134‑1 (MQ=255)

                                                                                                                                                                    agTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAAcc                                                                                                                                                                                                                              <  7:100129/149‑1 (MQ=255)

                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGTAAACGCCCATTCACTCACGCACAACCTGAACCGCAGACCAACAAACCCTCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTTTTCGTGGCCCTTCGGAAGGCTTCGAGCAATTGCCAAGTTAAAGGACGGGCTTACCAGTGATTGTTATCCGacacc                                                                                              >  4:52610/1‑274 (MQ=255)

                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCCTTCGTTCCTGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGa                                                                                                                            >  2:87861/1‑245 (MQ=255)

                                                                                                                                                                        gCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                                                                                                                                                          >  5:108339/1‑149 (MQ=255)

                                                                                                                                                                         ccTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                                                                                                                                                         >  6:955998/1‑149 (MQ=255)

                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                                                                                                                                                    <  4:301882/121‑1 (MQ=255)

                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                                                                                                                                                        >  5:884230/1‑149 (MQ=255)

                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGtat                                                                                                                                                                                              <  4:706731/175‑1 (MQ=255)

                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTCCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTAGCTTGAAACGCCCATACACTCACGCACAACCAGAACAGCAAACCAACAAACCATCATGGCATTCGTTCATGGGGGGCGCGGTATATCTCAAAGGGTGATGTCTTATTTGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGTGTATCCGCCACGaagg                                                                                          <  2:43922/275‑3 (MQ=255)

                                                                                                                                                                               gTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAtt                                                                                                                                                                                                                   >  5:538652/1‑149 (MQ=255)

                                                                                                                                                                                    aCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTc                                                                                                                                                                                                              >  8:90870/1‑149 (MQ=255)

                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGTGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            <  7:511873/149‑1 (MQ=255)

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

                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            >  8:1019209/1‑149 (MQ=255)

                                                                                                                                                                                      agaagaAGTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAGCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            <  8:652380/149‑1 (MQ=255)

                                                                                                                                                                                       gaagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATg                                                                                                                                                                                                           <  8:773187/149‑1 (MQ=255)

                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgggg                                                                                                                                                                                                        >  2:712258/1‑151 (MQ=255)

                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGcgcg                                                                                                                                                                                                  >  4:127785/1‑157 (MQ=255)

                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCa                                                                                                                                                                                        >  4:745060/1‑167 (MQ=255)

                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATattgca                                                                            <  2:499898/275‑1 (MQ=255)

                                                                                                                                                                                         agaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATggg                                                                                                                                                                                                         <  6:606160/149‑1 (MQ=255)

                                                                                                                                                                                          gaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                                               <  5:83968/94‑1 (MQ=255)

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

                                                                                                                                                                                              tattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcg                                                                                                                                                                                                    >  5:1032028/1‑149 (MQ=255)

                                                                                                                                                                                                   cgccgcCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGta                                                                                                                                                                                               <  5:244165/149‑1 (MQ=255)

                                                                                                                                                                                                   cgccgcCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCAAGGAATTCGTTCATGGGGAGCGCGGta                                                                                                                                                                                               <  8:70477/149‑1 (MQ=255)

                                                                                                                                                                                                       gccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtc                                                                                                                                                                                           <  6:1021021/149‑1 (MQ=255)

                                                                                                                                                                                                       gccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtc                                                                                                                                                                                           <  6:836164/149‑1 (MQ=255)

                                                                                                                                                                                                        ccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCt                                                                                                                                                                                                                                                  >  4:432462/1‑93 (MQ=255)

                                                                                                                                                                                                                  cGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCCTTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGa                                                                                                                                                                                >  5:327987/1‑149 (MQ=255)

                                                                                                                                                                                                                   ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAt                                                                                                                                                                               <  8:170993/149‑1 (MQ=255)

                                                                                                                                                                                                                   ggCTGGTTGCGAGCATTAACAATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGa                                                                                                                                                                                                                                                >  4:701687/1‑84 (MQ=255)

                                                                                                                                                                                                                     cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAttt                                                                                                                                                                             <  7:420031/149‑1 (MQ=255)

                                                                                                                                                                                                                      tGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTc                                                                                                                                                                            <  7:174685/149‑1 (MQ=255)

                                                                                                                                                                                                                          tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactc                                                                                                                                                                                                                                                             >  1:53298/1‑64 (MQ=255)

                                                                                                                                                                                                                          tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAt                                                                                                                                                                        <  5:960103/149‑1 (MQ=255)

                                                                                                                                                                                                                           gCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt                                                                                                                                                                       <  6:257319/149‑1 (MQ=255)

                                                                                                                                                                                                                            cGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                                                                                                                                                      <  8:946027/149‑1 (MQ=255)

                                                                                                                                                                                                                                cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                                                                                                                                                  <  8:21794/149‑1 (MQ=255)

                                                                                                                                                                                                                                 aTTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTggg                                                                                                                                                                 >  5:420697/1‑149 (MQ=255)

                                                                                                                                                                                                                                  ttAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTGCATGGGGAGCGCCGTATATCTCAAAGGGTGATTTCTTATTGGTGGGc                                                                                                                                                                <  6:327987/149‑1 (MQ=255)

                                                                                                                                                                                                                                      ccATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                                                                                                                                                            >  6:121622/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                          gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTgggg                          >  2:101460/1‑275 (MQ=255)

                                                                                                                                                                                                                                           cTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGa                         <  2:246216/275‑1 (MQ=255)

                                                                                                                                                                                                                                            tGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                                                                                                                                                      >  5:15162/1‑149 (MQ=255)

                                                                                                                                                                                                                                            tGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                                                                                                                                                      >  5:358990/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   cgcAATCATTTCCAGCAAGCTGGAAACGCTCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATCTTGCATTCCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTTGTTTT‑GGGACCCTTTccc                >  2:631570/1‑275 (MQ=255)

                                                                                                                                                                                                                                                    gcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                                                                                                                                                         >  8:384998/1‑74 (MQ=255)

                                                                                                                                                                                                                                                    gcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGc                                                                                                                                              >  5:1053221/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      gggcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa                                                                                                                                            >  8:496611/6‑149 (MQ=255)

                                                                                                                                                                                                                                                      aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGa                                                                                                                                                                                                                                        <  6:283438/57‑1 (MQ=255)

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

                                                                                                                                                                                                                                                       aTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTggg                           <  2:163916/261‑1 (MQ=255)

                                                                                                                                                                                                                                                        tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCCTTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa                                                                                                                                      >  1:144774/1‑153 (MQ=255)

                                                                                                                                                                                                                                                         cATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTg                                                                                                                                         >  6:531215/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                           tttCGAGCTTGCTGGAAACGCACATTCACCCACGCACAACCTTATCCGCAGACCAACACAACATCATGGCACTCGTCCATGGGGAGCGCGGTATAACTCAAAGGGTGATTTCTTATTTCGGGGCATTCTGAAAGCTTCGAGCAATTGcc                                                                                                                                       >  6:771211/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           tttCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                                                                                                                                                      >  5:859498/1‑118 (MQ=255)

                                                                                                                                                                                                                                                           tttCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                                                                                                                                                      <  6:859498/118‑1 (MQ=255)

                                                                                                                                                                                                                                                            ttCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCCATTGCCa                                                                                                                                      <  6:963999/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             tCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGCATTAGTTTTTGGGACCCTTTCCCTGCTTCggg       >  4:280503/1‑275 (MQ=255)

                                                                                                                                                                                                                                                              ccAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGc                                                                                                                                    >  6:695266/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               cAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcgc                                                                                                                                                                                                   >  7:678677/1‑85 (MQ=255)

                                                                                                                                                                                                                                                                  cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaaa                                                                                                                               >  2:400486/1‑150 (MQ=255)

TGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTCCCTGCTTCTGGAGCAA  >  NC_000913/2560946‑2561478

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