BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,171,852	C→T	5.1%	intergenic (+19/+27)	qseC → / ← ygiZ	quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,171,852	0	C	T	5.1%	616.7 / 12.8	191	intergenic (+19/+27)	qseC/ygiZ	quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand): ref base C (84/95); new base T (7/3); total (92/99)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 199 total aligned reads displayed.

AGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGA  >  NC_000913/3171592‑3172073

                   cACGATGAATGACGATAATGTCCTCGTGAGGGAGAACGGCCCCGGTGGTACACCAGATGCACTGGCGCGAATGTGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTGGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                                                                                                                                             <  1:6949/276‑1 (MQ=255)

                     cGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                                                                                                                                             <  2:360333/274‑1 (MQ=255)

                          aaTGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATTAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATCGCCAAata                                                                                                                                                                                       >  4:603199/1‑275 (MQ=255)

                                    aTTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                                                                                                                                                  >  2:599026/1‑238 (MQ=255)

                                                  gtataagagacaggTGGGACAACCGGGGAACTGGGGGAAATGGGCGTAAGCTTCAAGCGACCACCCGGACAAAGCGCTACCGGCAGCCTGATTGGGGTATTGTTTGTGCAGCGGAACACAAAATTGAATGACATGGATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGGCTATGAGAGCCGTTAAAACGACTCTCCTAGATTTTGCTGATAGCAAAATATAAACCGTCCCCAACAAAGccac                                                                                                                                                                <  4:569524/262‑1 (MQ=255)

                                                  ggtAACGGCCCGGGTTGGACACCCGCGGCACTGGCGCGAATGGGCGGACGCTTCCACCGCCCAGCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCTCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                                                                               <  4:419781/273‑1 (MQ=255)

                                                       cGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacca                                                                                                                                                          <  2:629582/275‑1 (MQ=255)

                                                                     cacCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCCTGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGGCACCAAAAAAGCCACCAACCCCAAACCAAAAAAAt                                                                                                                                            >  4:497269/1‑275 (MQ=255)

                                                                        cagggggaatGGCGCGCATTGGCGAACGTCTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                                                                                                                                          <  4:308843/265‑1 (MQ=255)

                                                                          gagGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                                                                                     >  1:184285/1‑261 (MQ=255)

                                                                            ggCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGATTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                                                                                     <  2:245247/275‑1 (MQ=255)

                                                                             gCACTGGTGCTAATTGGGTAACGCTTCTATCTGCCACCCGGACAAACCGCTACCGGCAGAGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                                                                                <  1:86283/279‑1 (MQ=255)

                                                                                   gcgcGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                                                                                                                                             <  2:567815/212‑1 (MQ=255)

                                                                                                                aCCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGc                                                                                                                                                                                                                               >  8:1378236/1‑149 (MQ=255)

                                                                                                                aCCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGc                                                                                                                                                                                                                               <  8:520711/149‑1 (MQ=255)

                                                                                                                   cGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGt                                                                                                                                                                                                                            <  5:103268/149‑1 (MQ=255)

                                                                                                                    ggACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                                                                                                                                           >  5:279295/1‑149 (MQ=255)

                                                                                                                      aCAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaa                                                                                                                                                                                                                         <  5:1395645/149‑1 (MQ=255)

                                                                                                                          aCCGCTACCGGCAGCGGGCTTGGGAAATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGGTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACg                                                                                                                                                                                                                     >  6:167913/1‑149 (MQ=255)

                                                                                                                                aCCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                                                                                                                                               <  5:934544/149‑1 (MQ=255)

                                                                                                                                 ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCa                                                                                                                                                                                                              <  8:958739/149‑1 (MQ=255)

                                                                                                                                  cGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                                                                                                                                             <  7:1070997/149‑1 (MQ=255)

                                                                                                                                  cGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                                                                                                                                             >  6:990772/1‑149 (MQ=255)

                                                                                                                                      acggggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATGTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCGGGTAACTTCACTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata                                                                            <  2:444593/271‑1 (MQ=255)

                                                                                                                                      aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCg                                                                                                                                                                                                                             <  1:267073/129‑1 (MQ=255)

                                                                                                                                      aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCg                                                                                                                                                                                                                             >  3:267073/1‑129 (MQ=255)

                                                                                                                                         gggCTTGGGCTATCGATTGTGCAGCGAATCGCCAAATTGCATGGCATGACGGTTGCATTTGGGAATTCGGAACAAGGTGGGTTTGAGGCGAAGGTGAGCTGGAAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtgtt                                                                                                                                                                                                      >  6:184966/1‑149 (MQ=255)

                                                                                                                                         gggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAtttt                                                                                                                                                                                                      >  6:943004/1‑149 (MQ=255)

                                                                                                                                          ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                                                                                                     <  7:814101/149‑1 (MQ=255)

                                                                                                                                           gCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTAc                                                                                                                                                                                                    <  6:458392/149‑1 (MQ=255)

                                                                                                                                            cTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                                                                                                                                   >  7:290072/1‑149 (MQ=255)

                                                                                                                                             ttGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                                                                                                                                  >  7:1557946/1‑149 (MQ=255)

                                                                                                                                               gggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                                                                                                      >  6:1282851/1‑127 (MQ=255)

                                                                                                                                               gggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAAc                                                                                                                                                                                                                      <  5:1282851/127‑1 (MQ=255)

                                                                                                                                                ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                               <  6:757400/149‑1 (MQ=255)

                                                                                                                                                 gCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                                                                                              >  5:598381/1‑149 (MQ=255)

                                                                                                                                                     tCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGACGTTAAAACGACTCTCATAGATTTTACTAATAGTAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATATCGAGCGCGACTCAATATACACAAAACTCGAACATTTCCCCCATCAAAAGTCTTATTAATATACGACTAATAAgctgtc                                                            >  4:444589/1‑270 (MQ=255)

                                                                                                                                                       gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                                                                                                        <  5:1425575/149‑1 (MQ=255)

                                                                                                                                                          tGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcac                                                                                                                                                         >  2:639764/1‑177 (MQ=255)

                                                                                                                                                             ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                                                                                                                                  >  7:98056/1‑149 (MQ=255)

                                                                                                                                                              cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                                                                                     <  4:639759/177‑1 (MQ=255)

                                                                                                                                                               aGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGcc                                                                                                                                                                                                                              >  5:756113/1‑103 (MQ=255)

                                                                                                                                                               aGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                                                                                                                                                >  7:1342483/1‑149 (MQ=255)

                                                                                                                                                               aGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGAGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGcc                                                                                                                                                                                                                              <  6:756113/103‑1 (MQ=255)

                                                                                                                                                                gCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGt                                                                                                                                                                               >  6:279629/1‑149 (MQ=255)

                                                                                                                                                                gCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGaaa                                                  <  2:497273/274‑1 (MQ=255)

                                                                                                                                                                 cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                                                                                                                                              <  6:1468125/149‑1 (MQ=255)

                                                                                                                                                                       gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                   >  8:765375/1‑138 (MQ=255)

                                                                                                                                                                        ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                                                                                                       <  6:598381/149‑1 (MQ=255)

                                                                                                                                                                        ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                                                                                                       >  6:1076664/1‑149 (MQ=255)

                                                                                                                                                                          aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                     <  5:137353/149‑1 (MQ=255)

                                                                                                                                                                            aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTAAGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                   <  8:729757/149‑1 (MQ=255)

                                                                                                                                                                                 atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                                                                                                                                              >  7:54381/1‑149 (MQ=255)

                                                                                                                                                                                 atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                                                                                                                                              >  5:519172/1‑149 (MQ=255)

                                                                                                                                                                                       tgaatgttgaatgtGGGAATGCGGAATAAGGTGGATTTGAGGCTAAGGTCATAAGGTAACTTCTCGGTCGATGAGAGCCGTTAAAACGGATCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCAACAACCAAAAACCAGGGAAAGACCGAGCGGGGATACATATACCAAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGGGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCaa                          <  4:689858/275‑1 (MQ=255)

                                                                                                                                                                                         aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                                                                                                                                      <  6:1390232/149‑1 (MQ=255)

                                                                                                                                                                                          aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                                                                                     <  5:71706/149‑1 (MQ=255)

                                                                                                                                                                                             ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                                                                                                                                  >  6:1628745/1‑149 (MQ=255)

                                                                                                                                                                                              tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGAt                                                                                                                                                                                                         >  5:1706123/1‑93 (MQ=255)

                                                                                                                                                                                                 aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                                                                                                                                              >  7:177089/1‑149 (MQ=255)

                                                                                                                                                                                                  tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                                                                                                                                             >  5:333554/1‑149 (MQ=255)

                                                                                                                                                                                                     gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                    >  5:220815/1‑75 (MQ=255)

                                                                                                                                                                                                       gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                                                                                                                                         <  8:1339048/148‑1 (MQ=255)

                                                                                                                                                                                                        aaTTCGGATCAAAGTGGGTTTGAGGCGAAAGTAAGCTGGTAGCTGCTCTGTCTATGAGAGACCTGAAAACGGATATCAGAGATTTTACTAATAGAAAAATAGAAACCGGCCCCAAAAAAGCCACCAAACACAAACCAGACAAAAACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAACAATAATAAGAGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCt          <  4:188316/274‑1 (MQ=255)

                                                                                                                                                                                                            cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                                                                                                   <  7:1517371/149‑1 (MQ=255)

                                                                                                                                                                                                            cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                                                                                                                                   <  6:150688/149‑1 (MQ=255)

                                                                                                                                                                                                            cGGAACAAGGTGGATGTGAGGTGATGGTAAACTGGGAACTTCTCTGTCTATGAGAGCCGTTAAAACGAATCTCACAGATGTTACTAATATCAAAATATAAACCGTCCCCAAAGAAGCCAGCATCCACAAACCAGACAAATACcgagcga                                                                                                                                   >  8:1045898/1‑149 (MQ=255)

                                                                                                                                                                                                               aaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                                                                                <  5:1678686/149‑1 (MQ=255)

                                                                                                                                                                                                               aaCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGa  <  2:301854/275‑1 (MQ=255)

                                                                                                                                                                                                                 cAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa                                                                                                                              >  5:1509798/1‑149 (MQ=255)

                                                                                                                                                                                                                     gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                                                                                                   <  7:1449206/124‑1 (MQ=255)

                                                                                                                                                                                                                                   aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGa                                                                                                            <  8:1468357/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                           >  7:438829/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                           <  8:1557946/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                           <  8:1432043/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                           <  8:589893/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaaa                                                                                                                 >  5:865465/1‑135 (MQ=255)

                                                                                                                                                                                                                                            tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTccc                                                                                                   >  8:1651286/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCg                                                                                                  >  7:1398627/1‑149 (MQ=255)

                                                                                                                                                                                                                                              gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                   >  6:377381/1‑67 (MQ=255)

                                                                                                                                                                                                                                              gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                   <  5:377381/67‑1 (MQ=255)

                                                                                                                                                                                                                                              gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                     >  6:1598443/1‑81 (MQ=255)

                                                                                                                                                                                                                                                aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                                                                                                                                         >  2:674065/1‑107 (MQ=255)

                                                                                                                                                                                                                                                aaCTTCTCTGTCAATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                                                                                                                                             <  8:375074/87‑1 (MQ=255)

                                                                                                                                                                                                                                                 aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATa                                                                                              >  6:152154/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAaga                                                                                          <  7:756273/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAaga                                                                                          >  8:878884/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtt                                                                                       >  7:63266/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtt                                                                                       >  7:1333701/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        tGTCTATGAGAGCCGTTAAAACGACTCTAATAGTTTGTACTAATAGCAAAATATAAACCGTCCGCAAAAAAGTCACCAACGAAAACCCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtc                                                                                       >  6:1229407/1‑148 (MQ=255)

                                                                                                                                                                                                                                                         gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGttt                                                                                      >  7:340076/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                     <  5:1655416/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAagag                                                                                         >  5:968153/1‑143 (MQ=255)

                                                                                                                                                                                                                                                            tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                   >  8:1154698/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                   >  8:1027512/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta                                                                                 >  6:248308/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaag                                                                               >  6:800343/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                                                                                                                                          >  4:145376/1‑86 (MQ=255)

                                                                                                                                                                                                                                                                    ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataa                                                                           <  5:780090/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta                                                                                 >  2:718472/1‑143 (MQ=255)

AGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGA  >  NC_000913/3171592‑3172073

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