BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,215	(G)_8→9	noncoding (72/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,215	1	.	G	84.2%	144.0 / 18.5	57	noncoding (72/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (9/0); new base G (26/22); total (35/22)
Fisher's exact test for biased strand distribution p-value = 9.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑G‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGT  >  NC_000913/4605965‑4606453

ttACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg‑g‑g‑gg‑gg‑g                                                                                                                                                                                                                                                    >  2:123818/1‑251 (MQ=255)

   cATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAAGTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTGGATATCGCCTCGTTTGGTGCGAg‑g‑g‑gg‑gg‑gga                                                                                                                                                                                                                                                  <  1:607634/251‑3 (MQ=255)

    aTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg‑g‑g‑gg‑gg‑ggact                                                                                                                                                                                                                                                <  1:115609/251‑5 (MQ=255)

            aCCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTA‑ATCTTACTTTTTTTGATATGGCCTCGTGTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑ATTGAAcccc                                                                                                                                                                                                                                      >  1:12063/1‑251 (MQ=255)

              cATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCa                                                                                                                                                                                                                                    <  1:123818/251‑1 (MQ=255)

               attgttccgggcatttatgggtgtgatcatGCTGATGATAAGGAATTTGAAATAGGTATGAAAAAACAAGTAAGAAAGAAGAAGGGAAGTAAAGAGAATTGGTGCGAGGGGGGTG‑ACTTGACCCCCCACGTCCGTAATGTCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCGCCTTCGCACAGGC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCAc                                                                                                                                                                                                                                   <  2:190743/226‑1 (MQ=21)

               aTCGTGCCTGGCATTTTTGGCTGGCAGCATGCTGATGATAAGGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTTGCGCGTCTACCAATGCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCAc                                                                                                                                                                                                                                   <  2:647953/251‑1 (MQ=37)

                   tGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑GCGG‑GG‑G‑‑G‑CTTGACCCCCCACGTcc                                                                                                                                                                                                                               >  2:237659/1‑251 (MQ=255)

                    gCCTGGCGATTATTGCTGGAATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAAATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCg                                                                                                                                                                                                                              <  2:119032/251‑1 (MQ=38)

                           aCTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCACCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑G‑‑A‑CTTGAAccccccagttctgaaggaac                                                                                                                                                                                                                      >  2:236473/1‑236 (MQ=255)

                           aCTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAAACAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAa‑c                                                                                                                                                                                                                      <  1:694320/251‑1 (MQ=255)

                            cTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAa‑c                                                                                                                                                                                                                      >  1:465546/1‑250 (MQ=255)

                              tATTGCTTGCATCATGCGGATGAGAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAGGGAAGTAAAAATAATTGGTGGCAGGGGGGGG‑ACTTGAACCCCCCCGTCCGTAAGGACACTAACACCTGAAGCTAGAGCATCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATAGCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAa‑cac                                                                                                                                                                                                                    <  2:208653/250‑1 (MQ=37)

                                     ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTGTGTGGTCGCTGCGTTTGTCTCGAG‑G‑G‑GG‑GGCG‑‑A‑CTTAAACCCCCAAGTACGTAAGAA‑CACTAacactt                                                                                                                                                                                                            >  2:151913/1‑249 (MQ=25)

                                                    atAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCt                                                                                                                                                                                             >  2:253521/1‑251 (MQ=255)

                                                      aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTATTGATATCGACTCGTTTGGTGCGAGTG‑G‑GG‑GG‑G‑‑A‑CTTGAACCCCCCCGTCCGTAAGAA‑CACTAACACCTGAAGCTAGcgcggcgac                                                                                                                                                                                           >  1:202166/1‑246 (MQ=255)

                                                            tttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACCTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                                                                                                                                     >  1:295462/1‑251 (MQ=255)

                                                                          gAAAAATCAAGGAAGAAACAAGAAAGGAATTAAAGATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACCCCTGAAGCTAGCGCGTCTACCAATTCCGCCGCCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GGGA‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                                                                                                                                        <  1:403016/251‑1 (MQ=37)

                                                                            naaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTGCGCACAGTC‑ATCTTACTTTTTTAGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GGTGG‑G‑‑A‑CTTGAACCCCTACGTCGGTAAGAA‑CACTAACACCTGAGGGTAGCGCGTCTATGACTGCCGCCACCTTCGccc                                                                                                                                                                       >  2:304022/2‑249 (MQ=37)

                                                                             aaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                                                                                                                                    >  1:444762/1‑251 (MQ=255)

                                                                                        gaaACAAGAAATGAAGTAAAGATAATTGGTGCGGGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCCTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGACTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat                                                                                                                                                         <  2:444762/251‑1 (MQ=255)

                                                                                              aGAAAGGAAGTACAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATAAACGt                                                                                                                                                   >  1:492407/1‑251 (MQ=255)

                                                                                               gAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg‑g‑g‑gg‑gg‑ggactt                                                                                                                                                                                                                                               >  1:283503/1‑156 (MQ=255)

                                                                                               gAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg‑g‑g‑gg‑gg‑ggactt                                                                                                                                                                                                                                               <  2:283503/161‑6 (MQ=255)

                                                                                                aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                                                                                                                                 <  2:671002/251‑1 (MQ=255)

                                                                                                          aaGACAATTGGTGCGATGTGGGGG‑TCTTGAACCCCCACGTCCGTAAGGTCACTAACACCTGAAGCTAGCGCGGCTACCAATTCCGCCACCTTCCCACAGTC‑ATCTTACTTTTTATGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGttc                                                                                                                                       <  1:72234/251‑3 (MQ=37)

                                                                                                             aTAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGAAACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTAGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGATTGGTGCGAG‑GAG‑GG‑GG‑G‑‑A‑CCTGAACCATCACGTCGGTAAAAAGCA‑TAGCACCCGATGGTAGCGGGTACACCAAATCCGACAAATTCGCCCAGTGAGCGAAATATTAACGGGGTTTTTgtggca                                                                                                                                      >  2:206146/1‑244 (MQ=12)

                                                                                                                       cGAGGGGGGGG‑AATTGAACCCCCACGTCCGTAAGGACACTAACTCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTGCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCTTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACt                                                                                                                          <  1:206146/251‑1 (MQ=255)

                                                                                                                         aGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGCCCGAAAGAA‑CACTAAAACCTGAAGCTAGTCCGTCTACCACTCCCGCCACCTCCACCCACTGTTACGAATATCAACGGGGGTTTTGGGGCGAGGGGGGGGATTTg                                                                                                                        >  1:21230/1‑251 (MQ=21)

                                                                                                                            gggggg‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGcgcg                                                                                                                                                                                                >  2:648331/1‑176 (MQ=255)

                                                                                                                            gggggg‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGcgcg                                                                                                                                                                                                <  1:648331/176‑1 (MQ=255)

                                                                                                                                gg‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCTATATCAACGTGGTTTTTGGGGCGAGGGGGGGGACTTGGAccccc                                                                                                                 >  1:170149/1‑251 (MQ=255)

                                                                                                                                 g‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTA‑ACCTTACTTTAGTTGATATCGCCTCGTTTGGGGCGAG‑G‑G‑GG‑GG‑G‑‑GGCTCGAACCCCCACGTCCGTAAGAA‑CACTAACACATGAAGATAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGCGAGCAATAGCAGCGTGGAGTTTGGTGCGAggggggggggatgacaaccc                                                                                                                 >  1:553482/1‑238 (MQ=25)

                                                                                                                                      tgacccccccCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATGCCGCCACCGTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGt                                                                                                             <  2:465546/246‑1 (MQ=38)

                                                                                                                                       gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCGTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                                                                                           <  1:237659/251‑1 (MQ=255)

                                                                                                                                       aagcccacACGTGCGTAAGGACACTATCAGCTGAAGGTAGCGTGTTAACCAATATCGCAGCCTTAGCACAGTCGAT‑TAACTTGTTTTGAAATCGCTTGGTTTGGTGGGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCACCGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                                                                                           <  2:626719/244‑1 (MQ=16)

                                                                                                                                            cccACGTCCGTAAGGCCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCC‑ATCTTACTTTTCTTGATATCGCCTCGTTTGGGGCGAG‑GTG‑AG‑GG‑G‑‑A‑CTTGAACTCCAACGTCCGTATGAA‑CACTAAACCCTGAAGATAGTGCGTCTACCAATTCCGCCACCTTCGCCAAGTGCGACAAATATCAACGTGGTATTTGGTGCGACGGGCAGGGCTTGAACCCACACGACCgaagg                                                                                                      >  1:198838/1‑247 (MQ=19)

                                                                                                                                            cccACGTCCGTAAGGACACTAACCCCTGAAGCTAGCGCGTCTACCCATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATAGCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGACCCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                      <  1:304022/251‑1 (MQ=255)

                                                                                                                                             ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GGGA‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                        >  1:670106/1‑249 (MQ=255)

                                                                                                                                               aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATATTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCCCGTCCGTAAGGac                                                                                                   >  1:602796/1‑251 (MQ=255)

                                                                                                                                                     gTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATGTCAACGTGGTTTTTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCAGAAAGGACACTAac                                                                                             >  1:115929/1‑251 (MQ=255)

                                                                                                                                                      tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTACGGACACTAac                                                                                             >  1:171107/1‑250 (MQ=255)

                                                                                                                                                                               cgcgTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                   <  1:253521/251‑1 (MQ=255)

                                                                                                                                                                                  gTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                >  1:110855/1‑251 (MQ=255)

                                                                                                                                                                                     tACCAATTCCGCCAACTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGccc                                                             >  1:183592/1‑250 (MQ=255)

                                                                                                                                                                                      aCCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAGCCCCACGTCCGTAAGGACACTAAACCCTGAAGATAGCGCGTCTACCAATTACGCCAc                                                            >  2:197184/1‑251 (MQ=38)

                                                                                                                                                                                             ccGCCACCTTCGCACAGTC‑ATCTTACTTAGTTTGATATCGACTCGTTTGGTGGGGG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                      <  1:151913/250‑1 (MQ=37)

                                                                                                                                                                                               gCCACATTCGCACAGTC‑ATCTTAATTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACTAATTCCGCCACCTTCGCCCAGCGCGAGGAATAGAAACGTGGTCTTTGGAGCGAGGGCGGTGAGTTGAACCCCCACGCCCATAAGGTAACAAACAACAGTAGCTAGCGAGCCTACTAATTCCGccgtcgtggtac                                                    >  2:299282/1‑240 (MQ=16)

                                                                                                                                                                                                       cGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt                                           <  1:197184/251‑1 (MQ=255)

                                                                                                                                                                                                        cgcccagccatattatgattgttGATATCGCCTCTGTTGGTCCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTACGAA‑CACTACCCCCTGATGCTAGCGCTTCTCCCAATTCCTCCACCATCGCCCAATGCGAGCGATATTAACGTGGTTTTTGTTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCGGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                            <  1:299282/229‑1 (MQ=18)

                                                                                                                                                                                                          acaGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCGAGCGCGTCTACCAATTACGCCACCTACGCATACCATCAATTCt                                        >  1:109827/1‑251 (MQ=255)

                                                                                                                                                                                                             gTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGCACACTAACACCTTCAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaa                                     >  2:443488/1‑250 (MQ=255)

                                                                                                                                                                                                                           ttttGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAc                        >  1:229351/1‑251 (MQ=255)

                                                                                                                                                                                                                                   gggggtttttggcgggggg‑g‑g‑gg‑gg‑g‑‑C‑CTTGACCCCCCCCCTCAGTAGGAA‑CCCTAAAACCTGAGGTTAGCGCATCTCCCAATGCCGCGCGCTTGGCGCATTCCGAGCTATATCACGGGGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg                <  2:472704/232‑1 (MQ=12)

                                                                                                                                                                                                                                     gCCTCGTTTGGTACGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg              >  1:494590/1‑251 (MQ=255)

                                                                                                                                                                                                                                      ctgcGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGGGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                          <  2:565243/155‑1 (MQ=255)

                                                                                                                                                                                                                                      ccTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcgc             >  1:175576/1‑251 (MQ=255)

                                                                                                                                                                                                                                      ccTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCACACGTCCGTAAGCA‑CAATAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                          >  1:565243/1‑158 (MQ=25)

                                                                                                                                                                                                                                       cTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑G‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAt           >  1:678895/1‑251 (MQ=255)

                                                                                                                                                                                                                                         agTTTGGTGCGGG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCGCGTGCGTAAGAA‑CCCTAACACCTGAAGCTAGCGCGTCTTCCAATTCCGCCAACTACGCCCAGTGCGAGCAATATAAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAtt          <  1:86080/250‑1 (MQ=37)

                                                                                                                                                                                                                                           tttGGTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTGCCAGCACTATCAACGTGGTTTTTGGTGCGAGGGGGGGGCTTTGGCCCCCCACGTCCGTAAGGATAGTAACACCTGGAGCTAGCGCGTCTACCAATTCCGACACCTTTGCATACGAACAATTCTTGAGAAGAATAGCTACCACGGACgcgcttttt        >  2:240614/1‑246 (MQ=21)

                                                                                                                                                                                                                                             ggggcgagg‑g‑g‑gg‑gg‑g‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt     <  1:57101/243‑1 (MQ=255)

                                                                                                                                                                                                                                               gTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGAGCGAGGCGCGGTACTTGGACCCCCACCTCCGCAAGGACACTAACACCTGAAGATAGCGCGTATACCAATTACACCACCTACGCGTAACATCAATTCTTAAAAAGAATTGCTACCACggccagtcgtcatattg    >  2:372145/1‑236 (MQ=21)

                                                                                                                                                                                                                                               gTGCGAG‑G‑G‑GG‑GG‑GG‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTACGAGCAATATCAACGAGGTTTTTGGTGCGAGGGGGGGGACTGGAACCCCCACGGCCGCAAGAAAAATAAAACATAAAGCGAACGCGTCTACCAATTCCGCCACCTTAGCAGACAATCAATTATTCAAAAGGATTGCTACCAcaaggacgcattatattg    >  1:204853/1‑234 (MQ=16)

                                                                                                                                                                                                                                                gcgagg‑g‑g‑gg‑gg‑g‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                 >  2:633788/6‑188 (MQ=255)

                                                                                                                                                                                                                                                gcgagg‑g‑g‑gg‑gg‑g‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                 <  1:633788/183‑1 (MQ=255)

                                                                                                                                                                                                                                                gcgagg‑g‑g‑gg‑gg‑g‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGt  >  1:448602/6‑251 (MQ=255)

TTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTC‑ATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑G‑G‑GG‑GG‑G‑‑A‑CTTGAACCCCCACGTCCGTAAGAA‑CACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGT  >  NC_000913/4605965‑4606453

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