BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,171,863	A→G	5.4%	intergenic (+30/+16)	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,863	0	A	G	5.4%	621.9 / 13.5	190	intergenic (+30/+16)	qseC/ygiZ	quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand): ref base A (85/93); new base G (6/5); total (92/98)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

Only 100 of 194 total aligned reads displayed.

AGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGA  >  NC_000913/3171592‑3172073

aGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGTATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCCCTGCCTATGAGCGCCGTTAAAACGActc                                                                                                                                                                                                                 >  2:50310/1‑275 (MQ=255)

  ggCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                                                                                                                                               >  2:689863/1‑275 (MQ=255)

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

                          aaTGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                       >  2:449588/1‑275 (MQ=255)

                                              gAGGGATAACGGCCCCGGTGTGACAACAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAAAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                                                                               <  1:69571/279‑1 (MQ=255)

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

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

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

                                                                        cAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                                                                                                                                          <  2:603203/274‑1 (MQ=255)

                                                                          gagGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAACCGACTCCCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                                                                                                    >  2:569528/1‑246 (MQ=255)

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

                                                                           agGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                                                                                      <  2:673148/275‑1 (MQ=255)

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

                                                                             gCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGACTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                                                                                     <  2:632853/274‑1 (MQ=255)

                                                                                                                           ccGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                    <  8:1521609/149‑1 (MQ=255)

                                                                                                                             gCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                                                                                                                                                  >  6:552902/1‑149 (MQ=255)

                                                                                                                              cTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                                                                                                 >  5:641591/1‑149 (MQ=255)

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

                                                                                                                                aCCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCCACAACCACAAACCAGACAAATAACGAGCGAGAATACATATACCCAAAACCCGAACATTTCCCGCACAAAGAATTTTCTta                                                                                 >  4:301850/1‑275 (MQ=255)

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

                                                                                                                                 ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                                                                                                                                              <  7:1415724/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                               <  7:1553008/149‑1 (MQ=255)

                                                                                                                                                 gCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                                                                                              <  7:922298/149‑1 (MQ=255)

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

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

                                                                                                                                                       gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                                                                                                        >  7:1468357/1‑149 (MQ=255)

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

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

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

                                                                                                                                                                       gCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTGCTAATAGCAAAATATaaa                                                                                                                                                                                   <  7:765375/138‑1 (MQ=255)

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

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

                                                                                                                                                                          aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                     <  8:1565357/149‑1 (MQ=255)

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

                                                                                                                                                                           aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATGTAAACCGTCCCCAAAAAAg                                                                                                                                                                    <  7:1378236/149‑1 (MQ=255)

                                                                                                                                                                           aaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                                                                                                    <  7:1227941/149‑1 (MQ=255)

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

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

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

                                                                                                                                                                                            gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                                                                                                   <  5:582070/149‑1 (MQ=255)

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

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

                                                                                                                                                                                               gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                                                                                                                                                >  5:193448/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                   gtaggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAAATTCTCTGTCTATGAGAGCCGTTAACCCGACTCGCATATATTTTACTCATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaag                                                                                                                                            <  6:137534/146‑2 (MQ=255)

                                                                                                                                                                                                    tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                                                                               <  1:50863/129‑1 (MQ=255)

                                                                                                                                                                                                       gAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAcc                                                                                                                                         >  7:1339048/1‑148 (MQ=255)

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

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

                                                                                                                                                                                                        aaTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                                                                                       >  6:1452931/1‑149 (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)

                                                                                                                                                                                                                   aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa                                                                                                                            <  6:33946/149‑1 (MQ=255)

                                                                                                                                                                                                                     gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                                                                                                   >  8:1449206/1‑124 (MQ=255)

                                                                                                                                                                                                                           ttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACGTATACCCa                                                                                                                    >  6:1480055/1‑149 (MQ=255)

                                                                                                                                                                                                                                gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                                                                                               >  8:338310/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                                                          gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt                                                                                                                               <  7:1054683/123‑1 (MQ=255)

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

                                                                                                                                                                                                                                            tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTGTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                                                                                                 >  5:1111483/1‑87 (MQ=255)

                                                                                                                                                                                                                                             ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCg                                                                                                  <  8:1398627/149‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                              gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGc                                                                                                 >  5:260548/1‑149 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                                                                                          tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                     >  5:737026/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                                                                                     <  8:473046/85‑1 (MQ=255)

                                                                                                                                                                                                                                                          tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                     <  6:641591/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                          tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                     >  6:39668/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTc                                                                                    <  5:1465222/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                            tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAagag                                                                                         <  6:968153/143‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                                               gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAag                                                                                           >  2:425527/1‑138 (MQ=255)

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

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

                                                                                                                                                                                                                                                                    ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                                                                                                                                          >  2:748718/1‑86 (MQ=255)

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

                                                                                                                                                                                                                                                                       ttAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaa                                                                        >  5:429682/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        tAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                                                                                      <  2:145378/86‑1 (MQ=255)

                                                                                                                                                                                                                                                                        tAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                                                                                      <  4:748712/86‑1 (MQ=255)

                                                                                                                                                                                                                                                                           aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaataat                                                                    >  6:480710/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             cGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAAtaataa                                                                  <  5:1557993/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                              gACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAg                                                                 <  7:1651286/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                               aCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGt                                                                <  6:1675726/149‑1 (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` `< 35 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`