BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,171,863	A→G	5.0%	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.0%	466.2 / 7.7	141	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 (67/64); new base G (5/3); total (73/68)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

Only 100 of 147 total aligned reads displayed.

CACAGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGA  >  NC_000913/3171589‑3172126

ccccgggTAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGGGGGATAACGGCCACGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  2:148974/271‑1 (MQ=255)

 acagggcaagtgggTAGACGTAAAGCTTAGTGGTGAGAAATTCGTCGTTAGGGAGAACGGCCCCGGTGTGACACCAGCGGCACTGGCGCGAATTGGCGAAAGCTTCTAACGCCCACCCAGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTTGATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  4:8864/262‑1 (MQ=255)

      gtaGCGTGGTAGCCGGCACGGCTAATCCTGATGATTTATTCGTGAGTGGTAACGGCCCCGGGGTGACACCAGAGGCACTGGCGCGCATGGGCGAACGCTTCTCTTGCCCACCCGGGCAAGCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAGTCGGCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGGACAAGGTGGATTTGAGACGAAGGTAGGCTGGTAATTTCTCTGTCTGTGAGAGCCGTTAAAACGGCTCTCa                                                                                                                                                                                                                                                                   <  4:255422/273‑1 (MQ=255)

                                                  aGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGGTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAAACCATCCCCAAAAAAGc                                                                                                                                                                                                                        >  1:164986/1‑274 (MQ=255)

                                                      ttAACGGCCCCGGTGTGACACCAGGGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                                                                                                                                    <  2:528839/273‑1 (MQ=255)

                                                              cccGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGATCTCGATGTTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTTGATTTGGGAATGCGGAACCAGGTGGATTTTAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGGCTCCCAAAGACTTTACTACTTGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCAaaaa                                                                                                                                                                                                           >  4:528839/1‑275 (MQ=255)

                                                               ccGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                                                                                                                                                      >  4:487156/1‑215 (MQ=255)

                                                                   tgtgACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCAt                                                                                                                                                                                                                                                                  <  2:487156/215‑1 (MQ=255)

                                                                                     ggCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt                                                                                                                                                                                    >  4:231529/1‑275 (MQ=255)

                                                                                                                cccAACCAGTCAAACCGCTTGCGGCTGCGGGTGTGGGTTTTCTGTTGTCCAGGAAATCGCCAAATTGAATGGCATGANTGTTGAATTTGGTAATGCGGAACAAGGTTGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTTTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAGAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                                                                                                                                                         <  3:164986/275‑1 (MQ=255)

                                                                                                                    cccGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  1:62531/159‑1 (MQ=255)

                                                                                                                      cGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  4:381771/157‑1 (MQ=255)

                                                                                                                                 cTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                                                                                                                                                                                      <  5:690722/149‑1 (MQ=255)

                                                                                                                                 cTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                                                                        <  2:231529/275‑1 (MQ=255)

                                                                                                                                    ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                                                                                                                                                                                                   <  5:391529/149‑1 (MQ=255)

                                                                                                                                     cGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCAt                                                                                                                                                                                                                                                                  <  6:98305/149‑1 (MQ=255)

                                                                                                                                             ggCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCAAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCAGTTAAGACGACTCTCATAGATTTTa                                                                                                                                                                                                                                                          <  6:476931/149‑1 (MQ=255)

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

                                                                                                                                                   ggCTATCGATTGTCCAGCGAATGGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTACGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                    <  8:890166/149‑1 (MQ=255)

                                                                                                                                                   ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                    <  8:80882/149‑1 (MQ=255)

                                                                                                                                                       aTCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaa                                                                                                                                                                                                                                                <  7:338728/149‑1 (MQ=255)

                                                                                                                                                        tCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGCGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                                                                                                                                                                                               <  5:600726/149‑1 (MQ=255)

                                                                                                                                                           aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                    <  5:104916/141‑1 (MQ=255)

                                                                                                                                                           aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                    >  6:104916/1‑141 (MQ=255)

                                                                                                                                                             tGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                        >  3:34157/1‑167 (MQ=255)

                                                                                                                                                               tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                                                                        <  7:49466/149‑1 (MQ=255)

                                                                                                                                                               tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                                                                                                                                                                                        >  7:960434/1‑149 (MQ=255)

                                                                                                                                                                 cAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAcc                                                                                                                                                                                                                                      >  8:160335/1‑149 (MQ=255)

                                                                                                                                                                  aGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCg                                                                                                                                                                                                                                     >  8:712290/1‑149 (MQ=255)

                                                                                                                                                                    cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                                                                                                                                                                                                   >  6:666512/1‑149 (MQ=255)

                                                                                                                                                                     gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTcc                                                                                                                                                                                                                                  >  5:334299/1‑149 (MQ=255)

                                                                                                                                                                      aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                                                                                                                                                                                                 <  6:279473/149‑1 (MQ=255)

                                                                                                                                                                         tgCCCAATTGCATGGGATGGAAGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTAGGAGAGCCGTTAAAACGACTCTCATAGATTGTACTAATAGCAAAAGATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAAtgc                                                                                                 <  4:127775/273‑1 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                          >  8:752284/1‑149 (MQ=255)

                                                                                                                                                                               aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                        >  7:804162/1‑149 (MQ=255)

                                                                                                                                                                                 tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                                                                                                                                                      <  7:310382/149‑1 (MQ=255)

                                                                                                                                                                                 tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                                                                                                                                                      <  5:843722/149‑1 (MQ=255)

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

                                                                                                                                                                                    atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCCCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                                                                                                                                                                                                   >  6:665785/1‑149 (MQ=255)

                                                                                                                                                                                     tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                                                                                                                                                                                                  <  5:836714/149‑1 (MQ=255)

                                                                                                                                                                                          tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                                                                                                                                                                                             >  5:694186/1‑149 (MQ=255)

                                                                                                                                                                                            aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                                                                                                                                                                                           >  8:485337/1‑149 (MQ=255)

                                                                                                                                                                                             aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                                                                                                                                                          <  6:192461/149‑1 (MQ=255)

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

                                                                                                                                                                                                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                                                                                                                                                                                              <  5:608352/110‑1 (MQ=255)

                                                                                                                                                                                                     tttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaaa                                                                                                                                                                                                  >  7:858310/1‑149 (MQ=255)

                                                                                                                                                                                                       tGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATa                                                                                                                                                                                                <  5:448593/149‑1 (MQ=255)

                                                                                                                                                                                                             tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                                                                                                                                                                                                     >  2:465237/1‑90 (MQ=255)

                                                                                                                                                                                                             tGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt                                                                                                                                       <  2:200035/200‑1 (MQ=255)

                                                                                                                                                                                                                ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                    <  8:628481/88‑1 (MQ=255)

                                                                                                                                                                                                                ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                                                                                                                                                       >  5:530585/1‑149 (MQ=255)

                                                                                                                                                                                                                 gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                                                                                                                                                                                      >  8:350545/1‑149 (MQ=255)

                                                                                                                                                                                                                      aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa                                                                                                                                                                                 <  8:214418/149‑1 (MQ=255)

                                                                                                                                                                                                                       ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat                                                                                                                                                                                >  7:625323/1‑149 (MQ=255)

                                                                                                                                                                                                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                                                                                                                                                                                                     >  7:592035/1‑63 (MQ=25)

                                                                                                                                                                                                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGctct                                                                                                                                                                                                                                                                     <  8:592035/63‑1 (MQ=25)

                                                                                                                                                                                                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata                                                                                                                                                                               >  8:8341/1‑149 (MQ=255)

                                                                                                                                                                                                                                gAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaa                                                                                                                                                                       >  6:699982/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaaa                                                                                                                                                                      <  6:800308/149‑1 (MQ=255)

                                                                                                                                                                                                                                 aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAt                                                                                                                                                             >  4:676983/1‑158 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                                                                                                                                                   <  5:768317/68‑1 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAg                                                                                                                                                                                                                                                   >  6:768317/1‑68 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaattcaa                                                                                                                                                                  <  8:907436/149‑6 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                                                                                                                                                                  >  5:698765/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                                                                                                                                                         <  2:676983/158‑1 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGaaa                                                           >  2:249120/1‑252 (MQ=255)

                                                                                                                                                                                                                                      aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGa                                                                                                                                                                 <  6:698765/149‑1 (MQ=255)

                                                                                                                                                                                                                                      aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAag                                                                                                                                                >  1:181299/1‑166 (MQ=255)

                                                                                                                                                                                                                                       aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa                                                                                                                                                                                   >  2:739639/1‑130 (MQ=255)

                                                                                                                                                                                                                                         gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACa                                                                                                                                                              >  8:305751/1‑149 (MQ=255)

                                                                                                                                                                                                                                         gTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGa                                                       <  4:249120/252‑1 (MQ=255)

                                                                                                                                                                                                                                          tAAGCTGGTAACTTCTCTGTATATGAGNGCAGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCAACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtt                                                                                                                                            <  3:181299/166‑1 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAt                                                                                                                                                             >  4:75585/1‑143 (MQ=255)

                                                                                                                                                                                                                                                 gTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                                                                                                                                                         >  4:354507/1‑82 (MQ=255)

                                                                                                                                                                                                                                                   aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATACCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCAt                                                                                                                                                    >  6:10628/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                                                                                                                                                                                                     <  2:354507/82‑1 (MQ=255)

                                                                                                                                                                                                                                                     cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCa                    <  2:561182/275‑1 (MQ=255)

                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                 >  5:312626/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                 >  7:237144/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tctcTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATACCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAag                                                                                                                                                <  5:10628/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACTTATACCCAAAACTCGAACATTTCCCGCATAAagag                                                                                                                                              <  7:350545/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         tctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAagag                                                                                                                                              <  7:804587/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                                                                                                                                          <  1:32396/102‑1 (MQ=255)

                                                                                                                                                                                                                                                              cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCACGCATAAAGAGTTTTc                                                                                                                                         <  8:102143/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              cTATGAGAACCGTTTAAACGACTCTCATAGACTTTACTAATTGCAAAATATAAACCGTCCCCAAAAAAGCAAACAACCACAAACCACACAAAAACCGAGCGAGAACACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTc                                                                                                                                         <  7:668037/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                                                                        >  5:917770/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATAc         >  2:80229/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                 tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTta                                                                                                                                      >  7:486717/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa                                                                                                                                     <  6:530585/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   agagTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaag                                                                                                                                    >  7:992746/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                                                                        >  8:1043932/1‑144 (MQ=255)

                                                                                                                                                                                                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                                                                        <  7:1043932/144‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata                                                                                                                                 >  8:449797/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCAttg                                                                                      >  4:666003/1‑191 (MQ=255)

                                                                                                                                                                                                                                                                         gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagataagg                                                                                                                              >  5:900113/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                          ttAAAACGCCTCTCATAGATTTTACTAATAGCAAAATATAATCCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGACTACATATACCCAAACCTCTAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaa                                                                                                                             <  7:979449/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          ttAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaa                                                                                                                             >  7:861649/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           tAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAACCCCTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaat                                                                                                                            <  5:110956/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                           tAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                                                                                                                                                                                      <  4:241143/75‑1 (MQ=255)

                                                                                                                                                                                                                                                                              aaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGaataat                                                                                                                         <  7:752284/149‑1 (MQ=255)

CACAGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGA  >  NC_000913/3171589‑3172126

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