BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,307,294	A→T	15.6%	L169Q (CTG→CAG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,294	0	A	T	15.6%	418.3 / 59.1	149	L169Q (CTG→CAG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand): ref base A (64/61); new base T (12/11); total (77/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 149 total aligned reads displayed.

CGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTG  >  NC_000913/4307022‑4307560

cGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTAACAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCCCCCAGGATTGCCACCCGCCCAGTGTTGGGGCATATCCCCCAg                                                                                                                                                                                                                                                                           >  4:422763/1‑274 (MQ=255)

                                            cccAGGACCACCGCATCGGGATTGAACAGATTAATGCTGGTGGCAATGGCCCTTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGTCTGTTTTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCa                                                                                                                                                                                                                              <  1:831/275‑1 (MQ=255)

                                                  aTCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTAAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATACCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                                                                                                                                                                                                                        <  2:283342/275‑1 (MQ=255)

                                                           tCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCNTACCAGCGTCTTAGCGCTATTCCCGAGCAATTGGTTTCCAGGCCCCCACGATTGCCCCACCCGCAGGTCNTGGCCATATCTCCCGGGGGCAGACACCCCAGTTCGCCCGCCACCCNCNCTGCACCCCGCCAcccgccccccccct                                                                                                                                                                                                                >  3:44328/1‑261 (MQ=255)

                                                                    aaCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAGTGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCCTTCCAGAGCAATTGGTTTCCAGGAACCCAGGATTGCCACACGCGCAGTGTTGGGCCATATCTCCCAGGGGGCTGTGACCCCGTTCGCCTGCCACACCGTGTCCACCCGTCCAAGGCGCACCGTTCATCcacccct                                                                                                                                                                                                       >  1:64439/1‑270 (MQ=255)

                                                                                   gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTACCACACGCCCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCAAGTTGGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTAATCCCCACTGCGAACCCCAACCCCg                                                                                                                                                                                        >  1:4157/1‑274 (MQ=255)

                                                                                                gTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCa                                                                                                                                                                                                                                                             >  2:278723/1‑192 (MQ=255)

                                                                                                    cgcgTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTc                                                                                                                                                                                                                                                         <  4:278722/192‑1 (MQ=255)

                                                                                                        ttttCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCCCACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                               >  2:473333/1‑230 (MQ=255)

                                                                                                                       gggACGAAAGGGGCG‑TTTCCGCATGGACGAAAGGATCGCGCAATGGGTAATTTCGGGGCTTTTGTTCGTACCAGCTCCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                                                                                                                                  <  2:430957/274‑1 (MQ=255)

                                                                                                                              aaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                                                                                                                                                                                                           >  4:158685/1‑212 (MQ=255)

                                                                                                                                  ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTg                                                                                                                                                                                                       <  2:158686/212‑1 (MQ=255)

                                                                                                                                     gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta                                                                                                                                     <  2:417239/275‑1 (MQ=255)

                                                                                                                                      ttgtCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGatat                                                                                                                                                                                                                                                                   >  5:969776/4‑148 (MQ=255)

                                                                                                                                       tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATg                                                                                                                                                                                                                                                                  <  6:870244/148‑1 (MQ=255)

                                                                                                                                        ttCCGCATAGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                 <  6:969776/148‑1 (MQ=255)

                                                                                                                                                 gACGAAAAGATCGCTCAATGGGTCATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgg                                                                                                                                                                                                                                                                          <  8:933853/130‑1 (MQ=255)

                                                                                                                                                 gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTAGGTCATATCTCCCAgg                                                                                                                                                                                                                                                                          >  7:933853/1‑130 (MQ=255)

                                                                                                                                                            cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTa                                                                                                              >  2:129308/1‑275 (MQ=255)

                                                                                                                                                                      gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                      <  6:1073161/113‑1 (MQ=255)

                                                                                                                                                                               gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                                                                                                                                         >  8:192793/1‑101 (MQ=255)

                                                                                                                                                                               gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                                                                                                                                         <  7:192793/101‑1 (MQ=255)

                                                                                                                                                                                  cTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCa                                                                                                                                                                                                                       <  6:27812/148‑1 (MQ=255)

                                                                                                                                                                                   tgttgtTCGTCCAAGCGTCTTAGCGCCATTCCAGGGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCg                                                                                                                                                                                                     <  4:185339/165‑1 (MQ=255)

                                                                                                                                                                                     ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                                                                                                                                                    >  7:177610/1‑148 (MQ=255)

                                                                                                                                                                                        ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                                                                                                                                                 <  8:486633/148‑1 (MQ=255)

                                                                                                                                                                                               cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                       <  6:1020314/87‑1 (MQ=255)

                                                                                                                                                                                               cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                                                                                                                                       >  5:1020314/1‑87 (MQ=255)

                                                                                                                                                                                                aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                >  5:473432/1‑141 (MQ=255)

                                                                                                                                                                                                aGCGTCTTAGCGCCATGCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                <  6:473432/141‑1 (MQ=255)

                                                                                                                                                                                                 gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                                                                                                                                                >  7:745142/1‑92 (MQ=255)

                                                                                                                                                                                                 gCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                                                                                                                                                <  8:745142/92‑1 (MQ=255)

                                                                                                                                                                                                    tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                                                                                                                                            <  7:708041/125‑1 (MQ=255)

                                                                                                                                                                                                       tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                                                                                                                                                                                                                                 >  7:437624/1‑85 (MQ=255)

                                                                                                                                                                                                       tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                                                                                                                                  >  5:562660/1‑148 (MQ=255)

                                                                                                                                                                                                           gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCg                                                                                                                                                                                                                                            >  8:448999/1‑102 (MQ=255)

                                                                                                                                                                                                           gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCg                                                                                                                                                                                                                                            <  7:448999/102‑1 (MQ=255)

                                                                                                                                                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                      >  7:557774/1‑75 (MQ=255)

                                                                                                                                                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                                                                                                                                      <  8:557774/75‑1 (MQ=255)

                                                                                                                                                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                                                                                                                                                 >  7:909748/1‑80 (MQ=38)

                                                                                                                                                                                                              aTTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                          >  6:913263/1‑101 (MQ=255)

                                                                                                                                                                                                              aTTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                          <  5:913263/101‑1 (MQ=255)

                                                                                                                                                                                                              aTTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                                                                                                                                           >  7:413636/1‑148 (MQ=255)

                                                                                                                                                                                                               ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTCCACCCGTCCACGGCGCACCGTTCATCCAAACTGCGAACCCCTTCCCCGTACCGAGATAGGCCACCAGGACCAGTTGTTGCGTAAGGCGGTTTCACAATGCGTCCCAGAGCAGTTACCGGCAAATCGCGCGGGAAAGCACCACCGGACTCTTCAGGCTATTTCCAAGCTTacct                                                           >  4:260099/1‑272 (MQ=255)

                                                                                                                                                                                                                    gagCAATTGGTTTCCAGGCACCCAGGATTGCCAAACGCGAAGTGTTGGGTCATATCTCCAAGGGGGATATGGCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCgag                                                      <  2:7458/275‑1 (MQ=255)

                                                                                                                                                                                                                     agCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc                                                                                                                                                                                    >  8:1054041/1‑148 (MQ=255)

                                                                                                                                                                                                                         aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                                                                                                                                                <  7:1054041/148‑1 (MQ=255)

                                                                                                                                                                                                                           tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCa                                                                                                                                                                                                                                      >  8:882908/1‑92 (MQ=255)

                                                                                                                                                                                                                           tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCg                                                                                                                                                                                   >  8:1106480/1‑143 (MQ=255)

                                                                                                                                                                                                                           tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCg                                                                                                                                                                                   <  7:1106480/143‑1 (MQ=255)

                                                                                                                                                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                          >  6:124431/1‑85 (MQ=255)

                                                                                                                                                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                                                                                                                                          <  5:124431/85‑1 (MQ=255)

                                                                                                                                                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgc                                                                                                                                                                           <  8:413636/148‑1 (MQ=255)

                                                                                                                                                                                                                                tCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                                                                                                                                           <  8:1042931/82‑1 (MQ=255)

                                                                                                                                                                                                                                tCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                                                                                                                                           >  7:1042931/1‑82 (MQ=255)

                                                                                                                                                                                                                                    ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggg                                                                                                                                                                                                                                                                         >  5:416072/1‑48 (MQ=255)

                                                                                                                                                                                                                                    ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAg                                                                                                                                                                     >  5:802576/1‑148 (MQ=255)

                                                                                                                                                                                                                                         ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                                                                                                                                                <  7:989088/52‑1 (MQ=18)

                                                                                                                                                                                                                                          cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                                                                                                                                      >  7:525671/1‑61 (MQ=255)

                                                                                                                                                                                                                                          cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                                                                                                                                               >  7:875706/1‑148 (MQ=255)

                                                                                                                                                                                                                                           aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                                                                                                                                                                                             >  5:270207/1‑101 (MQ=255)

                                                                                                                                                                                                                                            ggATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTTCCAGGAGAGTTGCAGGTTACAGTCGCGGGAAAACTCAACCGGACAATTTAGCGTATTTTCGCGCTTATCGGCGCGACCATATAAATCCGCCGCTGTTaa                              >  2:394992/1‑275 (MQ=255)

                                                                                                                                                                                                                                             gATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                            <  6:562660/148‑1 (MQ=255)

                                                                                                                                                                                                                                              aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                           <  6:802576/148‑1 (MQ=255)

                                                                                                                                                                                                                                               ttGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                                                                                                                                          <  8:177610/148‑1 (MQ=255)

                                                                                                                                                                                                                                               ttGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                                                                                                                                              <  6:1070305/144‑1 (MQ=255)

                                                                                                                                                                                                                                               ttGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                                                                                                                                              >  5:1070305/1‑144 (MQ=255)

                                                                                                                                                                                                                                                   cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                                                                                                                                           <  7:17300/47‑1 (MQ=255)

                                                                                                                                                                                                                                                   cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                                                                                                                                           >  8:17300/1‑47 (MQ=255)

                                                                                                                                                                                                                                                     cacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                                                                                                                                                   >  7:1076045/1‑133 (MQ=255)

                                                                                                                                                                                                                                                     cacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                                                                                                                                                   <  8:1076045/133‑1 (MQ=255)

                                                                                                                                                                                                                                                         cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg                                                                                                                                                                            >  6:28711/1‑120 (MQ=255)

                                                                                                                                                                                                                                                         cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg                                                                                                                                                                            <  5:28711/120‑1 (MQ=255)

                                                                                                                                                                                                                                                         cgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAg                                                                                                                                                                             <  6:1020246/119‑1 (MQ=255)

                                                                                                                                                                                                                                                           cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                                                                                                                                                >  7:827621/1‑114 (MQ=255)

                                                                                                                                                                                                                                                           cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                           >  5:766839/1‑135 (MQ=255)

                                                                                                                                                                                                                                                           cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                                                                                                                                           <  6:766839/135‑1 (MQ=255)

                                                                                                                                                                                                                                                             gtgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                                                                                                                                                               <  5:823654/129‑1 (MQ=255)

                                                                                                                                                                                                                                                                    gTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                                                                                                                                                         >  8:95872/1‑112 (MQ=255)

                                                                                                                                                                                                                                                                     tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                                                                                                                                            <  5:1003412/92‑1 (MQ=255)

                                                                                                                                                                                                                                                                     tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                              >  5:1100344/1‑106 (MQ=39)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                                                                                                                                           <  8:684445/92‑1 (MQ=255)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc                                                                                                                                          >  6:577765/1‑141 (MQ=255)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc                                                                                                                                          <  5:577765/141‑1 (MQ=255)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                                                                                                                                   >  6:583110/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                      cATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                                                                                                                                   >  5:8767/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCTTGCCACACCGTGTGTACCCCTCCCCGGCGCACCGTTCATCCACACTGCGCACCCCATCCCCGTACCTAGACAAGCCGCCTCAACCGGTTGTTTCGCAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCCGTTTAACGTCGCGGGAACACTCTTCAGGTCACTTCCCTGTACTTGCGAGCTTATCTGCGAGAGCACATAAGTCCGCCGCTGTTCACGGCAGCGTACGCGTaataaatatgtga  >  4:524449/1‑263 (MQ=255)

                                                                                                                                                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGACCTTATCGGCGAGATCATATAAATCCGCCCATGTTAACGGCAGGTTAGGCGTAGAcatacttggt  >  4:460136/1‑269 (MQ=255)

                                                                                                                                                                                                                                                                        tatCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                                                                                                                                                <  5:483764/69‑1 (MQ=38)

                                                                                                                                                                                                                                                                          tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                            >  8:201585/1‑119 (MQ=255)

                                                                                                                                                                                                                                                                          tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                            <  7:201585/119‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tcccccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc                                                                                                                               >  7:982698/4‑148 (MQ=255)

                                                                                                                                                                                                                                                                             cccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg                                                                                                                            <  6:8767/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                             cccAGGGGGATATGACCCAGTTCTCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg                                                                                                                            <  8:793187/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                              ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                              <  8:392137/97‑1 (MQ=255)

                                                                                                                                                                                                                                                                              ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                                                                                                                                              >  7:392137/1‑97 (MQ=255)

                                                                                                                                                                                                                                                                               cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                               <  6:1121500/63‑1 (MQ=255)

                                                                                                                                                                                                                                                                               cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                                                                                                                                               >  5:1121500/1‑63 (MQ=255)

                                                                                                                                                                                                                                                                               cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGt                                                                                                                                              >  6:455172/1‑128 (MQ=255)

                                                                                                                                                                                                                                                                               cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgag                                                                                                                          <  8:875706/148‑1 (MQ=255)

CGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTG  >  NC_000913/4307022‑4307560

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