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,207	1	.	G	86.1%	348.6 / 47.4	129	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (14/4); new base G (58/53); total (72/57)
Fisher's exact test for biased strand distribution p-value = 7.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

AGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑G‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAA  >  NC_000913/4606059‑4606352

aGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑g‑g                                                                                                                                                    >  1:1264179/1‑151 (MQ=255)

aGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTAATCATACTTTAAATGATGTCGCCCCATTTGG‑GGG‑G‑A‑g‑g                                                                                                                                                    >  2:1931529/1‑151 (MQ=17)

    aGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑g‑g‑g‑ggg                                                                                                                                              >  1:1877420/1‑151 (MQ=255)

           aaaGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑g‑g‑g‑gggg‑ggactt                                                                                                                                      <  2:463527/151‑6 (MQ=255)

           aaaGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑g‑g‑g‑gggg‑ggactt                                                                                                                                      <  1:1498930/151‑6 (MQ=255)

           aaaGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATCTTACTTTTATTGATATCCTCCCGCTTAG‑TGC‑G‑atgctaggggactt                                                                                                                                           >  2:1738215/1‑138 (MQ=17)

            aaGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTg                                                                                                                                      <  1:2016247/151‑1 (MQ=255)

            aaGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTATTATTACTTTTAATGATATGGTACCGATGGG‑GGG‑G‑A‑G‑GCG‑GGGG‑GACtt                                                                                                                                       >  2:1767121/1‑150 (MQ=14)

              gATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑G‑G‑G‑GGGG‑GACTTGAAc                                                                                                                                   >  1:1214741/1‑151 (MQ=255)

                 aaTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAAccc                                                                                                                                 <  2:1344080/151‑1 (MQ=255)

                  aTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GTCTTGAAcc                                                                                                                                  >  1:4993/1‑149 (MQ=255)

                   ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAAcccc                                                                                                                                >  1:549701/1‑150 (MQ=255)

                     ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑GCG‑G‑G‑G‑GGGG‑GACTTGAACCCCCAc                                                                                                                             >  1:1503600/1‑151 (MQ=255)

                     ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATAATACTTTTTTTGAATTCGCCTCATTGAG‑GGC‑G‑A‑G‑GCG‑GGGG‑GACTTCAACCCCCa                                                                                                                              >  2:1708606/1‑150 (MQ=17)

                       tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACTCCCa                                                                                                                              >  1:1019832/1‑148 (MQ=255)

                          gAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCg                                                                                                                        >  2:744502/1‑151 (MQ=255)

                           aggaggGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTc                                                                                                                          <  1:1417464/148‑1 (MQ=255)

                            ggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGt                                                                                                                       <  1:142684/150‑1 (MQ=255)

                            ggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTa                                                                                                                      <  2:212688/151‑1 (MQ=255)

                             gggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTaa                                                                                                                     <  2:1599997/151‑1 (MQ=255)

                               gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGa                                                                                                                   <  2:75641/151‑1 (MQ=255)

                                 gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAac                                                                                                                 <  2:287380/151‑1 (MQ=255)

                                  ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCACCTTACTTTTTTTGATATCCTCCCGTTTGG‑TGC‑G‑A‑G‑G‑C‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAacac                                                                                                               >  1:998188/1‑151 (MQ=255)

                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCCTCTCGTCTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAAGCCCAACGTCCGTAAGAACACt                                                                                                              >  2:1705326/1‑151 (MQ=25)

                                        gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAac                                                                                                           >  1:1531985/1‑150 (MQ=255)

                                        gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAac                                                                                                           <  1:56740/150‑1 (MQ=255)

                                        gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGCTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGCCCGTAAGAACACt                                                                                                              >  1:814997/1‑147 (MQ=255)

                                         aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTATGATATCGTATCGTTTAG‑GGC‑G‑A‑G‑G‑GTGGGG‑GATTTCAACTCTCAATTGCGAAAGAACACt                                                                                                              >  2:1882688/1‑146 (MQ=17)

                                           cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                                        <  2:1981501/150‑1 (MQ=255)

                                           cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                                       <  1:441593/151‑1 (MQ=255)

                                           cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATAACGCCTCGGTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACAccg                                                                                                       <  1:1092538/151‑2 (MQ=255)

                                                    ccGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                               >  1:1788532/1‑150 (MQ=255)

                                                        aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                                                                          >  1:1500414/1‑151 (MQ=255)

                                                        aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                                                                          >  1:122297/1‑151 (MQ=255)

                                                         aGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                                                                         <  1:1712204/151‑1 (MQ=255)

                                                          tatcgcTAACACCTGAAGCTAGCGAATCAACCAATTGCGCCACCTTCGCACAATCAACTTACTTTTTTTGATATCGACTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                                                                              <  2:245169/140‑1 (MQ=25)

                                                          ggccacTAACACCTTATGCTTGCTAGTCTTCCACTTCCGCCTCCGTAGCACAGTCATCTTACTTTTTTTGTTATTGACTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                                                                            <  2:83081/144‑1 (MQ=18)

                                                              acTAACACCTTAAGGTGGATCGTCTACCAAATCCGCCGCCGTCGCACAGTGAACTAACTTATATAGATATCGACTAGATTGGTTCC‑C‑G‑G‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                             <  2:526533/142‑1 (MQ=12)

                                                                tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTAG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTCAACTCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCTTCTACCa                                                                                    >  2:2070504/1‑149 (MQ=37)

                                                                 aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTCAACTCCCACGTCAGAAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                                 >  2:1604007/1‑151 (MQ=37)

                                                                   caccataaGCTAGCTAATCAACCAATTCCGCCACCGTCGAACAATGAACTAACTATTTTAGAAATCAAATCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                                <  2:886012/144‑1 (MQ=9)

                                                                   cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                               >  2:1692009/1‑151 (MQ=255)

                                                                   cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTTTGATATCGTCTCGTTTAG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACTCCCACGTCCGTATGAACACTAAAACCTGAAGCTATCGCGTCTACCAATTcc                                                                               >  2:969453/1‑151 (MQ=25)

                                                                    acCTGAAGCTAGCGCGTCTACCAATTCCGCCGCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                              <  2:1535120/151‑1 (MQ=255)

                                                                         aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAGCACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                                         >  2:1976221/1‑151 (MQ=255)

                                                                           gCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                                       <  1:784252/151‑1 (MQ=255)

                                                                              aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                    >  1:406892/1‑151 (MQ=255)

                                                                               gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                                   >  1:1000032/1‑151 (MQ=255)

                                                                                cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                    >  2:1516521/1‑149 (MQ=255)

                                                                                cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGAAATAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                      >  2:840000/1‑147 (MQ=255)

                                                                                 gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                                  >  1:637765/1‑150 (MQ=255)

                                                                                       aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                                           <  1:1711047/151‑1 (MQ=255)

                                                                                         cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            >  2:53389/1‑148 (MQ=255)

                                                                                          aaTTCCGCCGCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCa                                                        <  2:863775/151‑1 (MQ=255)

                                                                                          aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc                                                         >  2:1446635/1‑150 (MQ=255)

                                                                                              ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTCCGCCCAGTGCGAGCAata                                                     >  2:755005/1‑150 (MQ=255)

                                                                                               cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                                   <  2:1997349/151‑1 (MQ=255)

                                                                                                 ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                                                     >  2:925600/1‑147 (MQ=255)

                                                                                                  cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                                >  1:803826/1‑151 (MQ=255)

                                                                                                   aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                               <  1:1354626/151‑1 (MQ=255)

                                                                                                     cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                             >  2:1204098/1‑151 (MQ=255)

                                                                                                        cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                          <  1:2044172/151‑1 (MQ=255)

                                                                                                        cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                          >  2:640179/1‑151 (MQ=255)

                                                                                                           ataGTCATATTACCATGTTTGATATCGACTAGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                             <  2:1295895/143‑1 (MQ=18)

                                                                                                           acaGTGATATTCCTTTGTATGATATCGACTCGTTTTG‑TGG‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATGCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                       <  2:1430666/151‑1 (MQ=18)

                                                                                                           acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                       >  1:1650302/1‑151 (MQ=255)

                                                                                                           acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                        <  2:340637/150‑1 (MQ=255)

                                                                                                              gTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                      <  2:1090564/149‑1 (MQ=255)

                                                                                                               tCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGt                                    >  2:149117/1‑150 (MQ=255)

                                                                                                                cATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                       >  1:330043/1‑146 (MQ=255)

                                                                                                                cATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                                   >  1:639068/1‑150 (MQ=255)

                                                                                                                cATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                                  >  1:1402220/1‑151 (MQ=255)

                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATAAACGTGGTTTTTGGAGCg                                 >  2:1378678/1‑151 (MQ=255)

                                                                                                                  tCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                 >  1:1650638/1‑150 (MQ=255)

                                                                                                                    ttACCTTTTTTGAAAGCGCCTCGTTTGG‑TTC‑C‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                <  2:327622/149‑1 (MQ=21)

                                                                                                                     tACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                 <  2:1909381/147‑1 (MQ=255)

                                                                                                                     tACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                <  1:552847/148‑1 (MQ=255)

                                                                                                                      aCTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTACGAGCACTATCAACGTGGTTTTTGCAGCGAgggg                            >  2:2063733/1‑151 (MQ=21)

                                                                                                                        tttttttGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg                          >  2:891911/1‑151 (MQ=255)

                                                                                                                         ttttttGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGAGCGAgggggg                          >  2:940596/1‑150 (MQ=255)

                                                                                                                          tttttGATATCGCCTCGTTTTG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg                         <  1:1540714/150‑1 (MQ=255)

                                                                                                                           ttttGATATCTCCTCGTTTGG‑TGC‑G‑A‑GTG‑G‑GGGG‑GACTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggg                        <  2:413369/150‑1 (MQ=255)

                                                                                                                           ttttGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggtggggg                        <  2:1620016/150‑1 (MQ=255)

                                                                                                                            tttGATATCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGa                       <  1:1805659/150‑1 (MQ=255)

                                                                                                                               gaaatCGACTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACtt                    <  2:1567061/147‑1 (MQ=255)

                                                                                                                                 tatCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGa                  >  1:135116/1‑150 (MQ=255)

                                                                                                                                 tatCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAAGAGGTTTTTGGTGCGAGGGGGGGGAc                      >  2:1225718/1‑146 (MQ=255)

                                                                                                                                  atCGCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCGAGCAATATCAACGCGGTTTTTGGTGCGAGGGGGGGGACTTGa                  >  1:872361/1‑149 (MQ=255)

                                                                                                                                  atCGCCTCGTTTGG‑TGC‑C‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCGCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAc                <  2:648108/151‑1 (MQ=255)

                                                                                                                                     gCCTCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGGAcccc             >  2:1273171/1‑151 (MQ=255)

                                                                                                                                     gCCTCGTTTGG‑TGC‑G‑A‑G‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAccccc            >  1:1653286/1‑151 (MQ=255)

                                                                                                                                        tCGTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATTTCAACGTGGTTTTAGGTGCGAGGGGGGGGACTTGAACCCCCa           >  2:1546974/1‑150 (MQ=255)

                                                                                                                                          gTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGt        <  1:956721/151‑1 (MQ=255)

                                                                                                                                          gTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGt        >  2:1148013/1‑151 (MQ=255)

                                                                                                                                          gTTTGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTCGTGCGAGGGGGGGGACTTGAACCCCCAcct        >  2:873217/1‑149 (MQ=255)

                                                                                                                                           tttGG‑TGCGG‑G‑G‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCa           <  1:336728/147‑1 (MQ=255)

                                                                                                                                           tttGG‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc          <  1:1976959/148‑1 (MQ=255)

                                                                                                                                              gg‑TGC‑G‑AGG‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGt    <  2:233292/151‑1 (MQ=255)

                                                                                                                                               g‑TGC‑G‑A‑G‑G‑G‑GGGGACACTTGCACCCCCACCTCCTTAAGAACACCAACACCTGAAGATAGAGCGTATACCAATTCCGCCAACTTCGCCCAGAGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc      <  2:28288/148‑1 (MQ=25)

                                                                                                                                                     g‑A‑G‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTGCGAGCAATATCAACGAGGTTTTTGGTGCGAGGGGGGGGCCTTGAACCCCCACGTCCGTaa  >  2:957697/1‑147 (MQ=255)

AGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGG‑TGC‑G‑A‑G‑G‑G‑GGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAA  >  NC_000913/4606059‑4606352

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