BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,171,852	0	C	T	5.7%	513.1 / 14.7	158	intergenic (+19/+27)	qseC/ygiZ	quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand): ref base C (75/74); new base T (6/3); total (81/77)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

Only 100 of 166 total aligned reads displayed.

CACAGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGA  >  NC_000913/3171589‑3172126

cccagggcagggTGGTATACGTCACGCGTAATGCTGAGAATTTCAGCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  4:19401/273‑1 (MQ=255)

  cggggCAGGGTGGGCGCCGTAACGCTGNNTGCTGATGATTTCATCGTGGGGGATAACGGNCCCGGGGGGAAACCAGAGGCCCTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGNGGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  3:43234/271‑1 (MQ=255)

   aGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTTCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCCATGTCTATGAGAGCCGTTaaa                                                                                                                                                                                                                                                                             >  2:252999/1‑268 (MQ=255)

          cGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTTGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTTAGCTGGTAACTTCTCTGTCTATTGGAGCCGTTAAAACGACTCTCAAAGa                                                                                                                                                                                                                                                               >  2:282236/1‑275 (MQ=255)

                               tGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaa                                                                                                                                                                                                                                                                             >  4:219231/1‑240 (MQ=255)

                                                                                    tggccagAATTGGCGCACGCTTCTATCGCCCACCCGGACAAATCGCTACCGGCATCGGGCGTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAAAAAGGTGGATTTGAGGCGAAGGTAATCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa                                                                                                                                                                                     <  2:46002/269‑1 (MQ=255)

                                                                                         cgAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGa                                                                                                                                                                                                                                                                         <  4:272046/186‑1 (MQ=255)

                                                                                          gAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGAct                                                                                                                                                                                                                                                                       >  2:354930/1‑187 (MQ=255)

                                                                                              tGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                                                                                                                                                                                                   <  4:354930/187‑1 (MQ=255)

                                                                                                                 gcagccGGACAAACCGGCACCTGCCGCTGGCTTGGGTGATTGGTTGTGCAGGCAATCGCCAAATTGCATGGGCTCAATGGTGAATTTGGGAATTCGGAACAAGGTGGATTTGGGGCGAAGGTAAGCTGGTAACCTCTTTGTCTATGAGAGCCGTTAATACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAATGCCACCAACCGCTAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTcc                                                                                                                                                         <  4:208107/270‑1 (MQ=255)

                                                                                                                   aCCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGc                                                                                                                                                                                                                                                                                    >  7:1050741/1‑149 (MQ=255)

                                                                                                                       ggACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGtt                                                                                                                                                                                                                                                                                <  5:281943/149‑1 (MQ=255)

                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                            >  8:1124869/1‑149 (MQ=255)

                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                            >  7:237736/1‑149 (MQ=255)

                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTACATGGCATGAATGTTGGATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                            <  5:426496/149‑1 (MQ=255)

                                                                                                                           aaaCCGCTACCGGCAGCGGGCTTGGGCTATCGATTGACCAGCGAATCGCCAAATTGCATGGAATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                            <  5:1233577/149‑1 (MQ=255)

                                                                                                                               cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                            >  7:1146896/1‑145 (MQ=255)

                                                                                                                               cGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTaaaa                                                                                                                                                                                                                                                                            <  8:1146896/145‑1 (MQ=255)

                                                                                                                                    ccGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                                                                                                                                                                                                   <  6:498170/149‑1 (MQ=255)

                                                                                                                                       gCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTACGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAACCCACCAACCCAAAACCAGACAAAGACCGAGCGAAAAAACATATACCCAAAACTCAAACATTTCCCGCAAAACAAGTTTTCTtaagac                                                                                                                                  >  4:328033/1‑274 (MQ=255)

                                                                                                                                        cggCGGGCTTGGGCTATCGATTGTGCAGCGAATCGCCAAACTGCATGGCAGTAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagata                                                                                                                                 <  2:19621/273‑1 (MQ=255)

                                                                                                                                         aGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                                                                                                                                                                                              <  6:357834/149‑2 (MQ=255)

                                                                                                                                                  gggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                                                                                                                                                                                                     <  8:121235/149‑1 (MQ=255)

                                                                                                                                                   ggCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                                                                                                                                                                                                    >  7:1310439/1‑149 (MQ=255)

                                                                                                                                                      tATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                                                                                                                                                                                                 >  6:1189521/1‑149 (MQ=255)

                                                                                                                                                          gATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                                                                                                                                                                                             >  8:1288/1‑149 (MQ=255)

                                                                                                                                                           tttGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCGGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                            <  8:630730/148‑1 (MQ=255)

                                                                                                                                                           aTTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                            >  7:1346725/1‑149 (MQ=255)

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

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

                                                                                                                                                                    cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTc                                                                                                                                                                                                                                   <  7:839430/149‑1 (MQ=255)

                                                                                                                                                                     gAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAACAGCAAAATATAAACCGtac                                                                                                                                                                                                                                  <  6:1016192/149‑3 (MQ=255)

                                                                                                                                                                      aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                                                                                                                                                                                                 <  7:439964/149‑1 (MQ=255)

                                                                                                                                                                           ccAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                                                                                                                                                            >  7:506746/1‑149 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                                                                                                                                                          <  6:1201778/117‑1 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTa                                                                                                                                                                                                                                                          >  5:1201778/1‑117 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                          <  7:72525/149‑1 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                          <  7:342451/149‑1 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                          <  7:1006299/149‑1 (MQ=255)

                                                                                                                                                                             aaaTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaaa                                                                                                                                                                                                                          >  6:105297/1‑149 (MQ=255)

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

                                                                                                                                                                               aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                        >  8:1010796/1‑149 (MQ=255)

                                                                                                                                                                               aTTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGc                                                                                                                                                                                                                        <  5:54697/149‑1 (MQ=255)

                                                                                                                                                                                 tGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                                                                                                                                                                                      <  8:237736/149‑1 (MQ=255)

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

                                                                                                                                                                                   catggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacc                                                                                                                                                                                                                    >  6:378138/1‑149 (MQ=255)

                                                                                                                                                                                    atggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccacca                                                                                                                                                                                                                   >  7:1363598/1‑149 (MQ=255)

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

                                                                                                                                                                                      ggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaac                                                                                                                                                                                                                 <  6:371422/149‑1 (MQ=255)

                                                                                                                                                                                           gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGacagag                                                                                                                                                                                                                                                          >  8:884795/1‑98 (MQ=37)

                                                                                                                                                                                           gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaa                                                                                                                                                                                                            <  7:968055/149‑1 (MQ=255)

                                                                                                                                                                                           gAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaa                                                                                                                                                                                                            >  8:1253960/1‑149 (MQ=255)

                                                                                                                                                                                            aaTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                                                                                                                                                                                           <  8:194141/149‑1 (MQ=255)

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

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

                                                                                                                                                                                               gTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                                                                                                                                                        <  8:686900/149‑1 (MQ=255)

                                                                                                                                                                                                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaaa                                                                                                                                                                                                                           <  6:105949/129‑1 (MQ=255)

                                                                                                                                                                                                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAg                                                                                                                                                                                                       >  5:1071899/1‑149 (MQ=255)

                                                                                                                                                                                                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                                                                                                                                                                                      <  7:1224600/149‑1 (MQ=255)

                                                                                                                                                                                                           aaTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                                                                                                                                                                                            <  7:1078139/149‑1 (MQ=255)

                                                                                                                                                                                                            aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                                                                                                                                                           >  5:716844/1‑149 (MQ=255)

                                                                                                                                                                                                              gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                                                                                                                                                                                         >  7:1447863/1‑149 (MQ=255)

                                                                                                                                                                                                                   aCAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAATAAAGCCACCAACCACGAACCAGACAAAGACCGAGCGAGAAt                                                                                                                                                                                    <  6:567344/149‑1 (MQ=255)

                                                                                                                                                                                                                      aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa                                                                                                                                                                                 <  7:606676/149‑1 (MQ=255)

                                                                                                                                                                                                                       ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAAt                                                                                                                                                                                                                                                     >  3:27135/1‑80 (MQ=255)

                                                                                                                                                                                                                        gTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata                                                                                                                                                                               >  5:275361/1‑149 (MQ=255)

                                                                                                                                                                                                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                                                                                                                                          >  7:374408/1‑137 (MQ=255)

                                                                                                                                                                                                                         tGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                                                                                                                                                          <  8:374408/137‑1 (MQ=255)

                                                                                                                                                                                                                           gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                                                                                                                                                                                                 <  1:27135/80‑1 (MQ=255)

                                                                                                                                                                                                                           gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAc                                                                                                                                                                            <  8:823001/149‑1 (MQ=255)

                                                                                                                                                                                                                             tttGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc                                                                                                                                                                          >  8:1398604/1‑149 (MQ=255)

                                                                                                                                                                                                                             tttGAGGCGAAGGTAAGCCGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAccc                                                                                                                                                                          >  6:164565/1‑149 (MQ=255)

                                                                                                                                                                                                                                 aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccaca                                                                                                                                                                                                                    >  1:52923/1‑102 (MQ=255)

                                                                                                                                                                                                                                 aGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCaaaa                                                                                                                                                                      >  5:187741/1‑149 (MQ=255)

                                                                                                                                                                                                                                   gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACt                                                                                                                                                                    >  7:799443/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCg                                                                                                                                                                  >  7:837425/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aGGTAAGCTGTTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                                                                                <  8:1346725/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGaa                                                                                                                                                                <  7:1253960/149‑1 (MQ=255)

                                                                                                                                                                                                                                            aGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACAttt                                                                                                                                                           <  6:593630/149‑1 (MQ=255)

                                                                                                                                                                                                                                               tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTccc                                                                                                                                                        >  7:1091612/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                 >  6:800954/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                 <  5:1189521/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                 >  8:1072079/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      ttCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATaaa                                                                                                                                                 >  6:1132466/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        taaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAaga                                                                                                                                               >  7:432898/6‑149 (MQ=39)

                                                                                                                                                                                                                                                        ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAaga                                                                                                                                               >  6:1441045/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                            <  8:1023412/52‑1 (MQ=255)

                                                                                                                                                                                                                                                            gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                                                                                                                                                                                            >  7:1023412/1‑52 (MQ=255)

                                                                                                                                                                                                                                                             tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                                                                          >  7:812668/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                                                                          <  8:432898/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             tCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGtttt                                                                                                                                          <  8:837425/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               tATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCt                                                                                                                                        >  7:745376/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTc                                                                                                                                         >  7:608024/1‑148 (MQ=255)

                                                                                                                                                                                                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTc                                                                                                                                         <  8:608024/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCtt                                                                                                                                       >  7:244688/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gagaGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa                                                                                                                                     >  8:312697/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gagaGCCGTTAAAACGTCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa                                                                                                                                     <  7:1141776/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaa                                                                                                                                     >  5:559847/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTACCGc                                                                                                                                                      >  7:232685/1‑129 (MQ=255)

                                                                                                                                                                                                                                                                       ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAg                      >  2:125327/1‑255 (MQ=255)

CACAGGGCAGCGTGGTAGACGTCACGCTGAATGCTGATAATTTCATCGTGAGGGATAACGGCCCCGGTGTGACACCAGAGGCACTGGCGCGAATTGGCGAACGCTTCTATCGCCCACCCGGACAAACCGCTACCGGCAGCGGGCTTGGGCTATCGATTGTCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGATAAGAATAATAAGTGGCGTAAGAAGAAAAAATGCTGCAATCATTGCTTGCAATATACGATCTTCGCCTGAAACGGAAGGTATACTGCAAATATTTCCTCCCGACTCAAGCATCCAGTCATACTTATTGA  >  NC_000913/3171589‑3172126

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