BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	3,815,859	Δ82 bp	100%		[rph]–[rph]	[rph], [rph]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3815859	3815940	82	171 [0]	[1] 169	[rph]–[rph]	[rph], [rph]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3815858	0 (0.000)	136 (1.190)	79/158	NT	100%	pseudogene (24/48 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
*	?	NC_000913	3815941 =	0 (0.000)	136 (1.190)	79/158	NT	100%	pseudogene (609/669 nt)	rph	ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Only 100 of 169 total aligned reads displayed.

CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGC  >  NC_000913/3815842‑3816039

tggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaagatcagtttgccagcgccgccttctgcc                                                                                                    >  1:17146‑M2/100‑100 (MQ=255)

tggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcc                                                                                                    <  1:3127500‑M2/1‑1 (MQ=255)

 ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                   <  1:2657353‑M2/2‑1 (MQ=255)

 ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                   <  1:2573092‑M2/2‑1 (MQ=255)

 ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                   <  1:2260522‑M2/2‑1 (MQ=255)

 ggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCg                                                                                                   <  2:1179054‑M2/2‑1 (MQ=255)

     tcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTccc                                                                                               >  1:130645‑M2/95‑100 (MQ=255)

      catgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTcccc                                                                                              <  1:2743105‑M2/7‑1 (MQ=255)

       atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCt                                                                                             <  2:140400‑M2/8‑1 (MQ=255)

          ccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCa                                                                                          <  2:1480742‑M2/11‑1 (MQ=255)

           cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                         >  2:649417‑M2/89‑100 (MQ=255)

           cttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCAc                                                                                         <  1:459708‑M2/12‑1 (MQ=255)

                 tcctcatcttacttttctacagacaaaaaaaaggcgactcataagtcgcttcaaaaatcagtttgcctgcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                   >  2:2735018‑M2/83‑100 (MQ=255)

                  cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                  >  1:1837066‑M2/82‑100 (MQ=255)

                  cctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                  >  1:2861596‑M2/82‑100 (MQ=255)

                  cctcatcttacttttcaacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                    <  2:510551‑M2/17‑1 (MQ=255)

                   ctcatcttactttaaaaaaaaaaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                   <  2:3324207‑M2/18‑1 (MQ=255)

                        gatctttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                            <  1:1226461‑M2/25‑1 (MQ=255)

                        cttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                            >  1:2415354‑M2/76‑100 (MQ=255)

                         ttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                            <  2:2476430‑M2/25‑1 (MQ=255)

                          tacttttctacagccaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                          >  1:2446644‑M2/74‑100 (MQ=255)

                          tacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGccc                                                                          >  2:1506315‑M2/74‑100 (MQ=255)

                            cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                        >  1:1440666‑M2/72‑100 (MQ=255)

                            cttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagattgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGt                                                                        >  2:1035074‑M2/72‑100 (MQ=255)

                             ttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatgatgcgcagatc                                                                       >  2:1226453‑M2/71‑95 (MQ=255)

                                tctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                    <  1:1142223‑M2/33‑1 (MQ=255)

                                  tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                  >  2:2109429‑M2/66‑100 (MQ=255)

                                  tacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                  <  2:3077922‑M2/35‑1 (MQ=255)

                                   acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                 <  2:2964841‑M2/36‑1 (MQ=255)

                                   acagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                 >  2:2876509‑M2/65‑100 (MQ=255)

                                    cagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                <  1:3059312‑M2/37‑1 (MQ=255)

                                    cagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                                                <  2:389048‑M2/37‑1 (MQ=255)

                                     agacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                               >  2:2435278‑M2/63‑100 (MQ=255)

                                      gacaaaaaaaaggcgactcatcagtcgccttaaacatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                              >  1:1869445‑M2/62‑100 (MQ=255)

                                      gacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                              >  1:2627978‑M2/62‑100 (MQ=255)

                                       acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                             >  2:778526‑M2/61‑100 (MQ=255)

                                       acaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                             <  1:2925426‑M2/40‑1 (MQ=255)

                                        gatttttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                                            >  2:2161696‑M2/60‑100 (MQ=255)

                                        caacaaaaaggcgactcatcattcgccttatagatcagttcgccagcgccgccttctgcCGGCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                              >  2:2925165‑M2/60‑98 (MQ=255)

                                          aaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                          >  1:2820612‑M2/58‑100 (MQ=255)

                                            tttaaggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAc                                                        >  2:3103079‑M2/56‑100 (MQ=255)

                                            aaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAc                                                        <  1:3151046‑M2/45‑1 (MQ=255)

                                             aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                       >  2:3236157‑M2/55‑100 (MQ=255)

                                               aaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                     <  2:3289159‑M2/48‑1 (MQ=255)

                                                aggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                    >  2:760992‑M2/52‑100 (MQ=255)

                                                 ggcgtctcatcagtcgccttaaaactcagtttgccagcgccgccttctgcCTTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                   <  1:2579817‑M2/50‑1 (MQ=255)

                                                 ggcgactgatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                   >  2:1571490‑M2/51‑100 (MQ=255)

                                                   cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                 <  2:484748‑M2/52‑1 (MQ=255)

                                                     actcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                               >  1:2957875‑M2/47‑100 (MQ=255)

                                                     acccatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                               <  1:970380‑M2/54‑1 (MQ=255)

                                                      ctcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                              <  2:2000771‑M2/55‑1 (MQ=255)

                                                       tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                             <  2:3502158‑M2/56‑1 (MQ=255)

                                                       tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                             <  1:3204368‑M2/56‑1 (MQ=255)

                                                        gatgagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                            >  2:976736‑M2/44‑100 (MQ=255)

                                                        catcagtcgccttaaagatcagtttgccagggccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTGTTCGGGCATCACTGCGTTCATGTCGGt                                            >  2:3266051‑M2/44‑100 (MQ=255)

                                                        catcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                            >  2:1862104‑M2/44‑100 (MQ=255)

                                                        catcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                            >  2:2837797‑M2/44‑100 (MQ=255)

                                                        catcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                            <  2:1428231‑M2/57‑1 (MQ=255)

                                                         atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                           >  1:1744496‑M2/43‑100 (MQ=255)

                                                         atcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                           <  1:1488868‑M2/58‑1 (MQ=255)

                                                            agtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGTTctct                                        >  2:222904‑M2/40‑100 (MQ=255)

                                                             gtcgccttaaaaatcagtttgccagcgccgactcctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                       >  2:2386355‑M2/39‑100 (MQ=255)

                                                              tcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                      <  1:760996‑M2/63‑1 (MQ=255)

                                                                gccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                    >  1:666147‑M2/36‑100 (MQ=255)

                                                                 ccttaaaaatcagtttgccggcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                   <  2:2370748‑M2/66‑1 (MQ=255)

                                                                 ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                   <  1:1049156‑M2/66‑1 (MQ=255)

                                                                 ccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                   <  1:1425499‑M2/66‑1 (MQ=255)

                                                                 ccttaaaaatcagtgtgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc                                   <  2:817479‑M2/66‑1 (MQ=255)

                                                                  cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                  >  1:1829245‑M2/34‑100 (MQ=255)

                                                                  cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                  >  1:3442159‑M2/34‑100 (MQ=255)

                                                                  cttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa                                  >  2:1148332‑M2/34‑100 (MQ=255)

                                                                    taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                >  2:1963431‑M2/32‑100 (MQ=255)

                                                                    taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                <  2:3093963‑M2/69‑1 (MQ=255)

                                                                    taaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                >  1:2572632‑M2/32‑100 (MQ=255)

                                                                            agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCaa                        >  1:1367594‑M2/24‑100 (MQ=255)

                                                                             gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                       >  2:1466779‑M2/23‑100 (MQ=255)

                                                                             gtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                       >  2:3395566‑M2/23‑100 (MQ=255)

                                                                              tttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACg                      <  2:63675‑M2/79‑1 (MQ=255)

                                                                              tttcccagcgccgccttctgcCGTCCCCCGCACTTCAATGATGCGCCCGTATTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTAAACg                      <  2:1784273‑M2/79‑1 (MQ=255)

                                                                                 gccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAt                   >  1:2225513‑M2/19‑100 (MQ=255)

                                                                                  ccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                  >  2:489613‑M2/18‑100 (MQ=255)

                                                                                   cagcgccgccttatgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                 >  1:2393208‑M2/17‑100 (MQ=255)

                                                                                    agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                <  1:3410250‑M2/85‑1 (MQ=255)

                                                                                    agcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                <  1:3315073‑M2/85‑1 (MQ=255)

                                                                                     gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCa               >  2:177948‑M2/15‑100 (MQ=255)

                                                                                       gccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGa             <  2:1972020‑M2/88‑1 (MQ=255)

                                                                                         cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc           >  1:2117282‑M2/11‑100 (MQ=255)

                                                                                         cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc           <  1:1718549‑M2/90‑1 (MQ=255)

                                                                                          gccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg          <  2:2632088‑M2/91‑1 (MQ=255)

                                                                                          gccttctgcCGTCCCCTGCACTTCAACGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg          <  1:3550832‑M2/91‑1 (MQ=255)

                                                                                          gccttctgcCGTCCCATGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg          <  1:3582425‑M2/91‑1 (MQ=255)

                                                                                           ccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGc         >  2:3477293‑M2/9‑100 (MQ=255)

                                                                                            cTTCTGCCGTTTCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTCATGTCGGTCTCTTCGGCACAGTCTTCACCCTATTCCACATCGCa        >  2:2979016/1‑100 (MQ=255)

                                                                                               ctgcCGTCCCCTGCGCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc     <  2:2479699‑M2/96‑1 (MQ=255)

                                                                                               ctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc     <  2:639889‑M2/96‑1 (MQ=255)

                                                                                               ctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAc     >  1:1581802‑M2/5‑100 (MQ=255)

                                                                                                tgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc    <  1:2876914‑M2/97‑1 (MQ=255)

                                                                                                 gcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg   >  1:2029133‑M2/3‑100 (MQ=255)

                                                                                                  cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTGCAGATCGCAAACCGc  >  1:3579875‑M2/2‑100 (MQ=255)

                                                                                                  cCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGc  >  2:1342156‑M2/2‑100 (MQ=255)

CCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGC  >  NC_000913/3815842‑3816039

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