BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,294	T→C	9.9%	intergenic (+33/+154)	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,294	0	T	C	9.9%	430.6 / 29.3	141	intergenic (+33/+154)	cspA/hokA	RNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand): ref base T (60/67); new base C (4/10); total (64/77)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

Only 100 of 144 total aligned reads displayed.

TAAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATC  >  NC_000913/3720034‑3720565

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)

               aTGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCCAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                                                                                                                                                    >  2:598615/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)

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

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

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

                                                                                                                  ttCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGc                                                                                                                                                                                                                                                                               >  8:365976/1‑149 (MQ=255)

                                                                                                                        gCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                                                                                                                                          >  6:499802/1‑148 (MQ=255)

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

                                                                                                                          tATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGcc                                                                                                                                                                                                                                                                              >  5:111048/1‑142 (MQ=255)

                                                                                                                          tATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGcc                                                                                                                                                                                                                                                                              <  6:111048/142‑1 (MQ=255)

                                                                                                                           aTCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCaaatggc                                                                                                                                                                                                                                                                       >  6:271930/1‑142 (MQ=255)

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

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

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

                                                                                                                               aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                                                                                                                                                                                                                                                   >  5:339894/1‑148 (MQ=255)

                                                                                                                               aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                                                                                                                                                  <  7:876130/149‑1 (MQ=255)

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

                                                                                                                               aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCaggg                                                                                                                                                                                                                                                                   <  5:838679/148‑1 (MQ=255)

                                                                                                                                gAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTACACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                                                                                                                                                 <  7:773789/149‑1 (MQ=255)

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

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

                                                                                                                                   cGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtct                                                                                                                                                                                                                                                               <  6:549260/148‑3 (MQ=255)

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

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

                                                                                                                                      tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaa                                                                                                                                                                                                                        >  2:150138/1‑184 (MQ=255)

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

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

                                                                                                                                          tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcga                                                                                                                                                                                                                                <  4:275586/172‑1 (MQ=255)

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

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

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

                                                                                                                                              aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGttt                                                                                                                                                                                                                                                   <  7:637250/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                                                                                aGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCg                                                                                                                                                                                                                                                                      <  5:573940/112‑1 (MQ=255)

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

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

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

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

                                                                                                                                                                         ccTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGta                                                                                                                                                                                                                         >  6:369718/1‑148 (MQ=255)

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

                                                                                                                                                                              aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaaa                                                                                                                                                                                                                   >  7:376590/1‑149 (MQ=255)

                                                                                                                                                                                cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                                                                                                                                <  7:328374/118‑1 (MQ=255)

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

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

                                                                                                                                                                                   cGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCta                                                                                                                                                                                                               <  5:499802/148‑1 (MQ=255)

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

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

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

                                                                                                                                                                                           gcgcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                                                                                                                                 <  5:231677/106‑1 (MQ=255)

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

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

                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTTCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGTAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTga                                                                                                                                                                                                <  6:286961/148‑1 (MQ=255)

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

                                                                                                                                                                                                  aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTga                                                                                                                                                                                                >  5:286961/1‑148 (MQ=255)

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

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

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

                                                                                                                                                                                                                     gtaacCAGCCTTTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTg                                              <  2:326856/275‑1 (MQ=255)

                                                                                                                                                                                                                      taacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAAt                                                                                                                                                                            >  6:754807/1‑148 (MQ=255)

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

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

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

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

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

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

                                                                                                                                                                                                                             ccTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATg                                                                                                                                                                     <  6:339894/148‑1 (MQ=255)

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

                                                                                                                                                                                                                                 tAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTACTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCaa                                                                                                                                                                <  8:140401/149‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaa                                                                                                                                                        <  8:376590/149‑1 (MQ=255)

                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgta                                                                                                                                                         <  5:369718/148‑1 (MQ=255)

                                                                                                                                                                                                                                         gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgta                                                                                                                                                         <  6:137138/148‑1 (MQ=255)

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

                                                                                                                                                                                                                                            tAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATaaa                                                                                                                                                                                        <  6:622809/114‑1 (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)

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

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

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

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

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

                                                                                                                                                                                                                                                                    tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGa                                                                                                                              >  5:634353/1‑148 (MQ=255)

TAAAGGTAATACACTATGTCCGGTAAAATGACTGGTATCGTAAAATGGTTCAACGCTGACAAAGGCTTCGGCTTCATCACTCCTGACGATGGCTCTAAAGATGTGTTCGTACACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCTTTCAATTTGCAGGCTAAAAACGCCGCCAGCTCATAACTCTCCTGTTTAATATGCAATTCACACAGTGAATCTCTTATCATCCAGGTGAAAAATAAAAGCGTGAAACAAATC  >  NC_000913/3720034‑3720565

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