BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,300	T→A	9.4%	intergenic (+39/+148)	cspA → / ← hokA	RNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,300	0	T	A	9.4%	456.9 / 27.6	138	intergenic (+39/+148)	cspA/hokA	RNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand): ref base T (57/67); new base A (3/10); total (60/78)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 141 total aligned reads displayed.

TAAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATCACT  >  NC_000913/3720034‑3720568

taAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAATATCCCTGCCATTTGGCgggg                                                                                                                                                                                                                                                                      >  2:524893/1‑275 (MQ=255)

taAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTATTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                                                                                                                                                                                                                                                      >  4:64057/1‑275 (MQ=255)

          aacatATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCAGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAt                                                                                                                                                                                                                                                            <  4:524891/271‑1 (MQ=255)

              tATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGACCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCCTTTGGCGGGGATTTTTTTATTTGt                                                                                                                                                                                                                                                        >  2:32557/1‑275 (MQ=255)

                tGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAACGGCTTCGCCTTCATCCCTCCTGACGATGGCTCTAAAGATGTGTTAGTACCCTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGAACGAAGGTCCGAAAGTGTCCTTCACCATCTAACGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTCACCAGCCTGTAATCTCTGCTTAAACGCGCAGAATCTAAGTCCCCTGCCTTTCGCAGGTGTTTTTTTTATtttttt                                                                                                                                                                                                                                                      >  4:271142/1‑275 (MQ=255)

                                        tAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatctaaccgt                                                                                                                                                                                                                              >  4:709512/1‑268 (MQ=255)

                                                        tGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattat                                                                                                                                                                                                              <  2:709514/275‑1 (MQ=255)

                                                         gACAAAGGCTTCGGCTTCCTCACTCCTCACGGTGGATCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGACGGCCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCTGCAGCTAGTAACGTAACCAGCCTGTAATCTATGCTTAAACGCACAGATGCTAAGATCCCTGCCAATGGGCGGGGATTTTTTTAGTTGGTTTAAGGAAATAAATAATCGATCCGTTAATAAATTCtattttt                                                                                                                                                                                                             >  4:621265/1‑275 (MQ=255)

                                                                  ttcggctttATCACGCCCGACCATGGCTCTAAAGATTCGTTCGTACACTTCTCTGCTATCCAGAACGATGGGTACAAATCTCTGGACGAAGGTCAGAAAGTGACCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgtg                                                                                                                                                                                                    <  4:32556/275‑1 (MQ=255)

                                                                  ttcggcttcATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTAGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgtg                                                                                                                                                                                                    <  2:621267/275‑1 (MQ=255)

                                                                                 ccTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTggg                                                                                                                                                                                     <  2:64058/275‑1 (MQ=255)

                                                                                          ggCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGACGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAAAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAataaataa                                                                                                                                                                                                                                       >  2:80200/1‑216 (MQ=255)

                                                                                              ctAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGACGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAAAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcga                                                                                                                                                                                                                                   <  4:80199/216‑1 (MQ=255)

                                                                                                                        gCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                                                                                                                                             >  5:646274/1‑148 (MQ=255)

                                                                                                                         cTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTgg                                                                                                                                                                                                                                                                           >  8:791207/1‑149 (MQ=255)

                                                                                                                           aTCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGc                                                                                                                                                                                                                                                                          >  6:595031/1‑148 (MQ=255)

                                                                                                                            tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtttt                                                                                                                                                                                                                                                     >  4:755963/1‑168 (MQ=255)

                                                                                                                            tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGt                                                                                                                                                                                    >  1:60034/1‑233 (MQ=255)

                                                                                                                               aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                                                                                                                                                     >  8:35610/1‑149 (MQ=255)

                                                                                                                                 aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAtt                                                                                                                                                                                                                                                                   <  7:891868/149‑1 (MQ=255)

                                                                                                                                 aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTTAAGAATAAATTTGGGTGCAATGAAAATGGGCAACGCCGTAAATCAGGCGGGAATAATTTCccg                                                                                                                                     >  4:326854/1‑275 (MQ=255)

                                                                                                                                 aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                    <  6:328268/148‑1 (MQ=255)

                                                                                                                                  aCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAAATCCCCGCCAAATGGCAGGGAtct                                                                                                                                                                                                                                                                  <  8:793942/149‑3 (MQ=255)

                                                                                                                                   cGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtct                                                                                                                                                                                                                                                                  <  5:271930/148‑3 (MQ=255)

                                                                                                                                   cGATGGTTACAAATCCCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtctt                                                                                                                                                                                                                                                                 <  8:843770/149‑1 (MQ=255)

                                                                                                                                    gATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttt                                                                                                                                                                                                                                                                <  8:206681/149‑1 (MQ=255)

                                                                                                                                     aTGGTTACAAATCTCTGGGCGAAGGTCAGAAAATGTCCTTCAACCACGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttt                                                                                                                                                                                                                                                                <  5:350372/148‑1 (MQ=255)

                                                                                                                                      tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAataaataa                                                                                                                                                                                                                                       >  2:275587/1‑172 (MQ=255)

                                                                                                                                      tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGAACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCTATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCTCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGaa                                                                                                                                >  1:142687/1‑275 (MQ=255)

                                                                                                                                          tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                                                                                                                                          <  7:361517/149‑1 (MQ=255)

                                                                                                                                          tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaa                                                                                                                                                                                                                       <  4:150137/184‑1 (MQ=255)

                                                                                                                                           aCAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtcttagattct                                                                                                                                                                                                                                                          <  5:595031/148‑8 (MQ=255)

                                                                                                                                             aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtc                                                                                                                                                                                                                                                                   <  8:841120/137‑2 (MQ=255)

                                                                                                                                             aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtc                                                                                                                                                                                                                                                                   >  7:841120/1‑136 (MQ=255)

                                                                                                                                              aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                    >  5:402187/1‑135 (MQ=255)

                                                                                                                                              aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                    <  6:402187/135‑1 (MQ=255)

                                                                                                                                                 ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                                                                                                                                    <  6:138168/148‑1 (MQ=255)

                                                                                                                                                 ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                                                                                                                                    <  6:122998/148‑1 (MQ=255)

                                                                                                                                                  tctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAg                                                                                                                                                                                                                                                  >  7:277959/1‑149 (MQ=255)

                                                                                                                                                        cGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaa                                                                                                                                                                                                                                               <  2:433333/146‑1 (MQ=255)

                                                                                                                                                         gAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaata                                                                                                                                                                                                                                            <  6:28156/148‑1 (MQ=255)

                                                                                                                                                            ggTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATcgcg                                                                                                                                                                                                                              >  4:423913/1‑159 (MQ=255)

                                                                                                                                                               cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     <  7:35610/149‑1 (MQ=255)

                                                                                                                                                               cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     >  7:140401/1‑149 (MQ=255)

                                                                                                                                                               cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     <  8:262301/149‑1 (MQ=255)

                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCg                                                                                                                                                                                                                                                                         >  6:573940/1‑112 (MQ=255)

                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                                                                                                                                              >  8:673176/1‑139 (MQ=255)

                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                                                                                                                                              <  7:673176/139‑1 (MQ=255)

                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaat                                                                                                                                                                                                                          <  2:423915/159‑1 (MQ=255)

                                                                                                                                                                    aGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatc                                                                                                                                                                                                                                 >  5:137138/1‑148 (MQ=255)

                                                                                                                                                                          cTTCACCATCGGAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaat                                                                                                                                                                                                                          >  7:528565/1‑149 (MQ=255)

                                                                                                                                                                                cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   >  8:328374/1‑118 (MQ=255)

                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                                                                                                                                                       <  7:80809/112‑1 (MQ=255)

                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCta                                                                                                                                                                                                                  <  8:528565/149‑1 (MQ=255)

                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCta                                                                                                                                                                                                                  <  7:365976/149‑1 (MQ=255)

                                                                                                                                                                                  tCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAtttt                                                                                                                                                                                                        >  2:816710/1‑159 (MQ=255)

                                                                                                                                                                                     aaaGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAAGAATCGATCGCGTAATAAAATCtatt                                                                                                                                                                                                                >  6:337392/1‑148 (MQ=255)

                                                                                                                                                                                          ggcgCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAATTCtattattat                                                                                                                                                                                                           <  6:646274/148‑1 (MQ=255)

                                                                                                                                                                                          ggcgCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattattat                                                                                                                                                                                                           >  6:872482/1‑148 (MQ=255)

                                                                                                                                                                                             gcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAtttt                                                                                                                                                                                                        <  5:872482/148‑1 (MQ=255)

                                                                                                                                                                                               tAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                                                                                                                                                     <  7:879226/149‑1 (MQ=255)

                                                                                                                                                                                               tAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtg                                                                                                                                                                                                      <  6:477455/148‑1 (MQ=255)

                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAaataaata                                                                                                                                                                                                                                        <  6:63287/111‑1 (MQ=255)

                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAaataaata                                                                                                                                                                                                                                        >  5:63287/1‑111 (MQ=255)

                                                                                                                                                                                                      ccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAatat                                                                >  2:359253/1‑275 (MQ=255)

                                                                                                                                                                                                                gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   >  6:281603/1‑86 (MQ=255)

                                                                                                                                                                                                                gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                   <  5:281603/86‑1 (MQ=255)

                                                                                                                                                                                                                gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTcaca                                                      <  4:359251/275‑1 (MQ=255)

                                                                                                                                                                                                                     gtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAAGAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                <  5:337392/148‑1 (MQ=255)

                                                                                                                                                                                                                        acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                >  6:52040/1‑145 (MQ=255)

                                                                                                                                                                                                                        acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                <  5:52040/145‑1 (MQ=255)

                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                                                                                                                                          <  6:505483/68‑1 (MQ=255)

                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                                                                                                                                          >  5:505483/1‑68 (MQ=255)

                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          >  6:344773/1‑148 (MQ=255)

                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          <  5:344773/148‑1 (MQ=255)

                                                                                                                                                                                                                           aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAAt                                                                                                                                                                         <  7:791207/149‑1 (MQ=255)

                                                                                                                                                                                                                             ccTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATg                                                                                                                                                                        <  5:754807/148‑1 (MQ=255)

                                                                                                                                                                                                                               tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                <  5:19084/138‑1 (MQ=255)

                                                                                                                                                                                                                               tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                                                                                                                                                >  6:19084/1‑138 (MQ=255)

                                                                                                                                                                                                                                      tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgat                                                                                                                                                                                                                                  >  7:233580/1‑81 (MQ=255)

                                                                                                                                                                                                                                      tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgat                                                                                                                                                                                                                                  <  8:233580/81‑1 (MQ=255)

                                                                                                                                                                                                                                       ctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          >  6:368100/1‑136 (MQ=255)

                                                                                                                                                                                                                                       ctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                                                                                                                                          <  5:368100/136‑1 (MQ=255)

                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                                                                                                                                                     >  7:742367/1‑107 (MQ=255)

                                                                                                                                                                                                                                           ttAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     <  8:55361/73‑1 (MQ=255)

                                                                                                                                                                                                                                           ttAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatc                                                                                                                                                                                                                                     >  7:55361/1‑73 (MQ=255)

                                                                                                                                                                                                                                            tAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATaaa                                                                                                                                                                                           >  5:622809/1‑114 (MQ=255)

                                                                                                                                                                                                                                             aaaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcga                                                                                                                                                                                                                                   >  4:546960/1‑73 (MQ=255)

                                                                                                                                                                                                                                               aaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                                                                                                                                                      >  5:301928/1‑148 (MQ=255)

                                                                                                                                                                                                                                                aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAgg                                                                                                                                                     >  5:258438/1‑148 (MQ=255)

                                                                                                                                                                                                                                                 gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATcgc                                                                                                                                                                                                                               <  2:546962/73‑1 (MQ=255)

                                                                                                                                                                                                                                                          cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       >  5:850609/1‑104 (MQ=255)

                                                                                                                                                                                                                                                          cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       <  6:850609/104‑1 (MQ=255)

                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       >  5:485452/1‑102 (MQ=255)

                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                                                                                                                                                       <  6:485452/102‑1 (MQ=255)

                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt                                                                                                                                         <  6:258438/148‑1 (MQ=255)

                                                                                                                                                                                                                                                            aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt                                                                                                                                         >  6:590875/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                 cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCANAACACTCCTGTTTAATATGCAATTCACACAGTGAATCCCTTATCATCCAGGGACAAAAAAAAAACGGGAaaaaaaaa     >  3:170694/1‑273 (MQ=255)

                                                                                                                                                                                                                                                                     gCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATCACt  <  1:170695/274‑1 (MQ=255)

                                                                                                                                                                                                                                                                        aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTcaca                                                      <  4:491141/219‑1 (MQ=255)

TAAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATCACT  >  NC_000913/3720034‑3720568

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