BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,500,613	T→C	5.2%	V29A (GTG→GCG)	mdeA →	methionine gamma‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,500,613	0	T	C	5.2%	810.7 / 16.1	246	V29A (GTG→GCG)	mdeA	methionine gamma‑lyase
Reads supporting (aligned to +/- strand): ref base T (134/99); new base C (4/9); total (138/108)
Fisher's exact test for biased strand distribution p-value = 8.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

ATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAG  >  NC_002947/1500464‑1500761

 ttttGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTg                                                                                                                                                     >  2:2613216/1‑150 (MQ=255)

   ttGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGcc                                                                                                                                                   >  1:2861090/1‑150 (MQ=255)

     ggAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCAc                                                                                                                                                 >  1:1385337/1‑150 (MQ=255)

     ggAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCAc                                                                                                                                                 >  2:2707783/1‑150 (MQ=255)

      gAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCAcc                                                                                                                                                <  2:423945/150‑1 (MQ=255)

       aaaaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCg                                                                                                                                               >  2:500015/1‑150 (MQ=255)

        aaaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCgg                                                                                                                                              >  2:688116/1‑150 (MQ=255)

        aaaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGGTGCCATGCGCGACTCCCACCACACCACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCgg                                                                                                                                              >  2:4143149/1‑150 (MQ=255)

         aaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCATACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCTCTTTCCCATGGTGGCACCAGGGCGCCAccatg                                                                                                                                             <  2:4553759/150‑4 (MQ=255)

         aaaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgt                                                                                                                                             <  1:387173/150‑1 (MQ=255)

          aaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgtg                                                                                                                                            <  2:3450917/150‑1 (MQ=255)

          aaaTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCACGGTGGCACCAGGGCGCCAccatgg                                                                                                                                            <  1:1603613/150‑5 (MQ=255)

            aTCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTa                                                                                                                                          <  2:2178951/150‑1 (MQ=255)

                  cAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGAc                                                                                                                                    <  1:1063066/150‑1 (MQ=255)

                   aGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACg                                                                                                                                   >  2:2469548/1‑150 (MQ=255)

                            gCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTAt                                                                                                                          >  2:5526090/1‑150 (MQ=255)

                            gCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTAt                                                                                                                          <  1:408538/150‑1 (MQ=255)

                             cATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATg                                                                                                                         <  2:602568/150‑1 (MQ=255)

                               tAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGcc                                                                                                                       >  1:4457501/1‑150 (MQ=255)

                                  gCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTc                                                                                                                    <  1:1814932/150‑1 (MQ=255)

                                    ttATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTccc                                                                                                                  >  2:3291609/1‑150 (MQ=255)

                                        aGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACt                                                                                                              >  1:2047426/1‑150 (MQ=255)

                                          cAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGACTTCCCCACTGt                                                                                                            >  2:114181/1‑150 (MQ=255)

                                               cacaAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGt                                                                                                       >  1:510451/1‑150 (MQ=255)

                                               cacaAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGt                                                                                                       >  2:3963259/1‑150 (MQ=255)

                                                    gagGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGt                                                                                                  >  2:2357073/1‑150 (MQ=255)

                                                       gCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGTGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGtgct                                                                                               >  1:5354072/1‑150 (MQ=255)

                                                       gCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGtgct                                                                                               <  2:5347298/150‑1 (MQ=255)

                                                       gCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGtgct                                                                                               >  1:4100881/1‑150 (MQ=255)

                                                             gCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGctgcctgc                                                                                         >  1:4661065/1‑150 (MQ=255)

                                                                  gcgcGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgtg                                                                                                                                            <  2:185145/94‑1 (MQ=255)

                                                                  gcgcGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGgtg                                                                                                                                            >  1:185145/1‑94 (MQ=255)

                                                                  gcgcGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGc                                                                                    >  1:5434091/1‑150 (MQ=255)

                                                                   cgcgACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGcc                                                                                   >  2:1501224/1‑150 (MQ=255)

                                                                    gcgACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCg                                                                                  >  2:2142224/1‑150 (MQ=255)

                                                                    gcgACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCg                                                                                  >  1:4789177/1‑150 (MQ=255)

                                                                      gACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCggg                                                                                >  1:4196488/1‑150 (MQ=255)

                                                                      gACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCgg                                                                                 <  2:3512294/149‑1 (MQ=255)

                                                                      gACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCagg                                                                                >  1:3370980/1‑150 (MQ=255)

                                                                       aCTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCgggg                                                                               >  2:3589280/1‑150 (MQ=255)

                                                                          cccACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggagg                                                                            <  1:2402588/150‑1 (MQ=255)

                                                                          cccACCACAACACCGGTTTTTCCACGCGGGCCATTCATCAGGGCTATGACCCGCTTTCCCCTGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTCCGGTGCTGCCTGCTTCGCCGGggagg                                                                            >  2:2566329/1‑150 (MQ=255)

                                                                           ccaccaCAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGggagga                                                                           <  2:5429285/150‑1 (MQ=255)

                                                                            caccacAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGCGTCCGGTGCTGCCTGCGTCGCCGGGgaggag                                                                          >  1:5502307/1‑150 (MQ=255)

                                                                              ccacAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGcc                                                                        >  2:1589057/1‑150 (MQ=255)

                                                                                cccacacCGGTTTGTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCCTGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTg                                                                      <  1:4211012/147‑1 (MQ=255)

                                                                                 caacaCCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTgg                                                                     <  1:3748274/150‑1 (MQ=255)

                                                                                   acacCGGTGTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGc                                                                   <  2:3201936/150‑1 (MQ=255)

                                                                                      ccGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCAtt                                                                >  2:3391235/1‑150 (MQ=255)

                                                                                      ccGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCAtt                                                                >  1:3688611/1‑150 (MQ=255)

                                                                                       cGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCTGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGATGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCAttt                                                               >  1:1287447/1‑150 (MQ=255)

                                                                                        ggTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTc                                                              >  2:181804/1‑150 (MQ=255)

                                                                                        ggTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTc                                                              <  1:3768707/150‑1 (MQ=255)

                                                                                         gTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCt                                                             <  1:3829718/150‑1 (MQ=255)

                                                                                            tttCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACa                                                          <  1:257617/150‑1 (MQ=255)

                                                                                            tttCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACa                                                          <  1:1505852/150‑1 (MQ=255)

                                                                                              tCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGc                                                        >  1:4584614/1‑150 (MQ=255)

                                                                                                 aCGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGc                                                     >  1:2265509/1‑150 (MQ=255)

                                                                                                 aCGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGc                                                     >  2:3776659/1‑150 (MQ=255)

                                                                                                 aCGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGc                                                     >  1:4423820/1‑150 (MQ=255)

                                                                                                   gcgGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCggg                                                                                <  2:4771060/121‑1 (MQ=255)

                                                                                                   gcgGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCACCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCAt                                                   >  2:2811421/1‑150 (MQ=255)

                                                                                                       gCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTcc                                               <  2:4013832/150‑1 (MQ=255)

                                                                                                         cATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTCCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCaa                                             <  1:2461969/150‑1 (MQ=255)

                                                                                                          aTTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAc                                            <  1:4169420/150‑1 (MQ=255)

                                                                                                          aTTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAc                                            <  1:4739333/150‑1 (MQ=255)

                                                                                                          aTTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAc                                            >  2:1142084/1‑150 (MQ=255)

                                                                                                           ttCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAcc                                           >  1:4940542/1‑150 (MQ=255)

                                                                                                           ttCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAAcc                                           >  2:323182/1‑150 (MQ=255)

                                                                                                              atcatGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCa                                        >  1:738626/1‑150 (MQ=255)

                                                                                                                 atGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGc                                     >  1:832738/1‑150 (MQ=255)

                                                                                                                 atGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGc                                     >  2:2014700/1‑150 (MQ=255)

                                                                                                                     cTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGc                                 <  1:3493001/150‑1 (MQ=255)

                                                                                                                          aCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGc                            <  2:90740/150‑1 (MQ=255)

                                                                                                                          aCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGc                            >  1:2058050/1‑150 (MQ=255)

                                                                                                                          aCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGc                            >  1:3585141/1‑150 (MQ=255)

                                                                                                                          aCCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGc                            <  2:4728693/150‑1 (MQ=255)

                                                                                                                           cccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCt                           >  2:4825994/1‑150 (MQ=255)

                                                                                                                            ccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                          >  2:708331/1‑150 (MQ=255)

                                                                                                                            ccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                          >  1:966246/1‑150 (MQ=255)

                                                                                                                            ccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                          >  1:1335197/1‑150 (MQ=255)

                                                                                                                            ccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                          >  1:45455/1‑150 (MQ=255)

                                                                                                                            ccGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTg                          >  1:3232171/1‑150 (MQ=255)

                                                                                                                              gCTTTCCTATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGa                        >  2:4601397/1‑150 (MQ=255)

                                                                                                                                   cccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGagcagc                   >  1:2346665/1‑150 (MQ=255)

                                                                                                                                   cccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCAGGagcagc                   >  2:1344373/1‑150 (MQ=255)

                                                                                                                                    ccATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAgcg                  >  2:4601907/1‑150 (MQ=255)

                                                                                                                                      aTGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCa                >  2:2997889/1‑150 (MQ=255)

                                                                                                                                      aTGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCa                <  1:2707783/150‑1 (MQ=255)

                                                                                                                                      aTGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCa                <  2:1385337/150‑1 (MQ=255)

                                                                                                                                      aTGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCa                <  2:841016/150‑1 (MQ=255)

                                                                                                                                       cggtggCACCAGGGCGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCAt               >  2:5060551/2‑150 (MQ=255)

                                                                                                                                       cggtggCACCAGGGCGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCAt               >  2:2045855/2‑150 (MQ=255)

                                                                                                                                         gtggCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATgg             <  1:4828089/150‑1 (MQ=255)

                                                                                                                                           ggCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCt           >  1:5309865/1‑150 (MQ=255)

                                                                                                                                             cgcCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTc         <  1:2521722/150‑1 (MQ=255)

                                                                                                                                               cccTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGc       <  2:5215536/150‑1 (MQ=255)

                                                                                                                                               cccTGGTGCCACCGGTGTACCAGACAGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGc       <  1:3066916/150‑1 (MQ=255)

                                                                                                                                                   ggTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGa   <  1:1771828/150‑1 (MQ=255)

                                                                                                                                                    gTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAg  <  1:294179/150‑1 (MQ=255)

ATTTTGGAAAAAATCGGGCAGCTCGGTGGCATAAGCTTATAGCAAACCACAAGAGGCTGTTGCCATGCGCGACTCCCACCACAACACCGGTTTTTCCACGCGGGCCATTCATCATGGCTATGACCCGCTTTCCCATGGTGGCGCCCTGGTGCCACCGGTGTACCAGACGGCTACCTATGCCTTCCCCACTGTCGAGTACGGTGCTGCCTGCTTCGCCGGGGAGGAGCCTGGGCATTTCTACAGCCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAG  >  NC_002947/1500464‑1500761

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