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	88.5%	278.7 / 28.6	96	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (7/4); new base G (44/41); total (51/45)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

AGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCC  >  NC_000913/4606059‑4606348

aGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg                                                                                                                                               <  2:1103218/151‑1 (MQ=255)

 gAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g                                                                                                                                             <  4:1507797/151‑1 (MQ=255)

   aaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑gg                                                                                                                                          <  3:283742/151‑1 (MQ=255)

    aGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑ggg                                                                                                                                         <  4:492498/151‑1 (MQ=255)

     ggAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCTTCGCTTGGTGC‑GA‑‑gg‑g‑gggg                                                                                                                                        >  2:482607/1‑151 (MQ=255)

      gAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACGGAACATACTTTTGTAGGGATTGCTTCTTTTGTGGC‑GA‑‑gg‑g‑gggg                                                                                                                                        >  2:588644/1‑150 (MQ=14)

       aaGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑gggg                                                                                                                                        >  4:230737/1‑149 (MQ=255)

        aGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑gggggga                                                                                                                                     <  2:1103173/151‑3 (MQ=255)

        aGTAAAGATAATTGGGGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑gggggga                                                                                                                                     <  2:120467/151‑3 (MQ=255)

         gTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTCACACCTGAAGCTAGCGCGTCTACCAATTCCGCAACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTGTAGTGC‑GA‑‑gg‑g‑ggggggac                                                                                                                                    >  2:1424329/1‑148 (MQ=255)

           aaaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑ggggggactt                                                                                                                                  <  4:624763/151‑6 (MQ=255)

           aaaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑ggggggact                                                                                                                                   <  3:271327/150‑5 (MQ=255)

           aaaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑ggggggact                                                                                                                                   <  2:1195402/150‑5 (MQ=255)

           aaaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTATTGATATCGCATCGTGTGGGGG‑GA‑‑gg‑g‑tgggggtctt                                                                                                                                  >  2:1417239/1‑146 (MQ=18)

           aaaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATATTACTTTATTTGATATCCTCTCGCTGAGGGC‑GA‑‑GG‑GTGGGGGACtt                                                                                                                                   >  1:473841/1‑151 (MQ=17)

            aaGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑gg‑g‑ggggggactt                                                                                                                                  <  4:595736/151‑6 (MQ=255)

             aGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGa                                                                                                                                 <  3:237691/151‑1 (MQ=255)

             aGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGa                                                                                                                                 <  3:440402/151‑1 (MQ=255)

               aTAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGaa                                                                                                                                <  4:780576/150‑1 (MQ=255)

                 aaTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGaa                                                                                                                                >  4:574477/1‑148 (MQ=255)

                 aaTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GG‑G‑GGGGGACTTGAAcccc                                                                                                                            >  4:1265152/1‑151 (MQ=255)

                      gTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATTTCGTCTCGCTTGGGGC‑GA‑‑GGCT‑GGGGGACTTCAAc                                                                                                                               >  2:950699/1‑144 (MQ=18)

                      gTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACTCACATTACTTTTTTTGATTTCGTCCCGTTTGGTGG‑GA‑‑GGCT‑AGGGGATTTGAAc                                                                                                                               >  1:284376/1‑144 (MQ=14)

                        gCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTc                                                                                                                      <  2:730934/151‑1 (MQ=255)

                         cGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGt                                                                                                                       >  2:948334/1‑149 (MQ=255)

                         cGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTcc                                                                                                                     >  4:133753/1‑151 (MQ=255)

                            cggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCg                                                                                                                    >  4:673674/2‑150 (MQ=255)

                             ggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCGCCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTa                                                                                                                  <  2:1215154/151‑1 (MQ=255)

                              gggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCACATTACTTTTTTTGAAATCGTCTCGCTGAGGGC‑GA‑‑GG‑G‑GGGGGACtttcacc                                                                                                                              >  2:1672406/1‑132 (MQ=17)

                               ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                                <  2:1393295/151‑1 (MQ=255)

                                gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                               >  4:1419096/1‑150 (MQ=255)

                                  gagACTTGAACCCCCACGTCCGTAAGGACTCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                                             <  2:859107/149‑1 (MQ=255)

                                     aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                           <  4:416952/150‑1 (MQ=255)

                                      cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCCCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                          >  3:442733/1‑150 (MQ=255)

                                        tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGAGC‑GAG‑GG‑G‑GGGGGCCTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                        >  2:1623946/1‑150 (MQ=255)

                                            cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                                   <  4:1296619/151‑1 (MQ=255)

                                              cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                                                                 <  4:1032340/151‑1 (MQ=255)

                                                     ccGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                                                                          <  4:132622/151‑1 (MQ=255)

                                                      cGTAAGGACACTAACACCTGATGCTAGCTCGTCTACCAATTCCGCCACCGTAGCACAGGCAACTTACTTTTTTTGATATCGACTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                         <  1:437759/151‑1 (MQ=25)

                                                       gTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCACATTACTTTTTTTGATATCGCCCCGTTTAGTGC‑CA‑‑GG‑GTGGGGGACTTCAACCCCCACGCCCTTAAAACCACTAACACCTCAAGCTAacgc                                                                                        >  4:563954/1‑151 (MQ=21)

                                                             acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCGCCTTCGCACAGTCATCTTTCTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                  <  4:582906/151‑1 (MQ=255)

                                                             acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATAGCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                  <  4:260079/151‑1 (MQ=255)

                                                              cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                                 >  4:787301/1‑151 (MQ=255)

                                                               acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGGCCGTAAGa                                                                                                               >  3:304083/1‑120 (MQ=255)

                                                                cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCACCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                                <  2:946326/150‑1 (MQ=255)

                                                                 tAACCCCTGAAGCCAGCGCGTCTACCAATTCCGCCTCCTTCGCCCAGTCATCTTCCTTTCTTTGATCCCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                            <  1:304066/121‑1 (MQ=25)

                                                                    cacCTGAAGCTAGCGAGTCTACCAATTCCGCCACCTTCGCACAGTGATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                           <  2:1523639/151‑1 (MQ=255)

                                                                              tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                                 >  3:372231/1‑151 (MQ=255)

                                                                               aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                >  4:592904/1‑151 (MQ=255)

                                                                                 cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                              >  4:283708/1‑151 (MQ=255)

                                                                                  gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                              >  2:1405774/1‑150 (MQ=255)

                                                                                   cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGTTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                              >  2:835172/1‑149 (MQ=255)

                                                                                    gTCTACCAATACCGCCACCTAAACACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                                                             <  4:1102538/149‑1 (MQ=21)

                                                                                     tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                            <  4:567653/149‑1 (MQ=255)

                                                                                      cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                         >  4:1196623/1‑151 (MQ=255)

                                                                                         ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGACCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                        >  1:32575/1‑149 (MQ=255)

                                                                                          cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc                                                     >  1:307170/1‑151 (MQ=255)

                                                                                            aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa                                                   <  4:123830/151‑1 (MQ=255)

                                                                                              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  3:33172/1‑151 (MQ=255)

                                                                                              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  1:345058/1‑150 (MQ=255)

                                                                                              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  2:1504343/1‑150 (MQ=255)

                                                                                              tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  4:1169235/1‑150 (MQ=255)

                                                                                                cgcggccTACGCACAGTCATCTTACTTATTTTGATATCACCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                                                 <  2:1425778/144‑1 (MQ=21)

                                                                                                cGCCACTTTCGAACAGTCATCTTACTTTTTTTGATATTGCCTCGTTTGGTGCGGG‑‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                                                 <  2:1569671/149‑1 (MQ=21)

                                                                                                    aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                           >  4:744298/1‑151 (MQ=255)

                                                                                                      cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                           >  1:13628/1‑150 (MQ=255)

                                                                                                       ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                        >  4:226866/1‑151 (MQ=255)

                                                                                                         cGCACAGACATATTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa                                             <  2:406612/144‑1 (MQ=255)

                                                                                                           cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCCACCAATTCCGCCACCTTCCCCCAGCGCGAGCATTATCAACGttgatt                                     >  2:150503/1‑146 (MQ=21)

                                                                                                           cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                    >  2:445579/1‑151 (MQ=255)

                                                                                                            acagagaactaaCTATTTATGATATCGCATAGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                    <  2:164769/139‑1 (MQ=18)

                                                                                                              aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                 <  2:896507/151‑1 (MQ=255)

                                                                                                               gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                <  3:180379/151‑1 (MQ=255)

                                                                                                                tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                               >  2:1463534/1‑151 (MQ=255)

                                                                                                                 cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                 <  2:40435/148‑1 (MQ=255)

                                                                                                                   tCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                             >  4:1201128/1‑150 (MQ=255)

                                                                                                                    cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                             >  4:1504952/1‑149 (MQ=255)

                                                                                                                     ttACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgg                          >  4:750283/1‑151 (MQ=255)

                                                                                                                       aCTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                        >  2:1520601/1‑151 (MQ=255)

                                                                                                                           tttttGATATCGCCTCGTTTGGTGC‑GA‑‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggg                    <  3:414840/150‑1 (MQ=255)

                                                                                                                            ttttGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggg                    >  4:1125560/1‑150 (MQ=255)

                                                                                                                            ttttGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGAGCGAgggg                        >  2:64869/1‑146 (MQ=255)

                                                                                                                             tttGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACCCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGa                   <  2:1234125/150‑1 (MQ=255)

                                                                                                                             tttGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGa                   <  4:457359/150‑1 (MQ=255)

                                                                                                                              ttGATAACGCCTCTTTTTGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGAc                  <  1:420723/150‑1 (MQ=255)

                                                                                                                               tgaaagcGCCTCGGTTGGTGCGGG‑‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGa                   <  1:2321/142‑1 (MQ=21)

                                                                                                                               tGATATCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACtt                >  4:1578745/1‑151 (MQ=255)

                                                                                                                                gATATCACCGCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTg               <  2:260114/151‑1 (MQ=255)

                                                                                                                                   atcacttgtgttGGTGCGGG‑‑GG‑G‑GGGGGACTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAc            <  2:1302932/141‑1 (MQ=21)

                                                                                                                                    tCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGaa             <  2:133528/149‑1 (MQ=255)

                                                                                                                                    tCGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAc            <  4:1407533/150‑1 (MQ=255)

                                                                                                                                     caccTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAccc          <  4:110263/149‑1 (MQ=255)

                                                                                                                                     cGCCTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAccc          >  4:728706/1‑151 (MQ=255)

                                                                                                                                        cTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCACTATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAccccc        >  4:73928/1‑150 (MQ=255)

                                                                                                                                        cTCGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAccccc        >  3:38637/1‑150 (MQ=255)

                                                                                                                                          cGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCa       <  4:949570/149‑1 (MQ=255)

                                                                                                                                          cGTTTGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc      <  4:1055933/150‑1 (MQ=255)

                                                                                                                                            tttGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc   >  4:1321271/1‑151 (MQ=255)

                                                                                                                                              ttgTGC‑GAG‑GG‑G‑GGGGGACTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc  <  2:1007769/148‑1 (MQ=255)

                                                                                                                                              tGGTGC‑GAG‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGTAATATCAACGAGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc   >  2:584503/1‑149 (MQ=255)

AGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GG‑G‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCC  >  NC_000913/4606059‑4606348

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