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	135 [3]	[1] 133	[rph]–[rph]	[rph], [rph]

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

Only 100 of 135 total aligned reads displayed.

TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCA  >  NC_000913/3815760‑3815956

tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                   >  2:2359886‑M1/1‑99 (MQ=255)

 ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                  <  2:1214459‑M1/100‑3 (MQ=255)

 ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                  <  2:3162079‑M1/100‑3 (MQ=255)

 ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                  <  2:2752905‑M1/100‑3 (MQ=255)

  gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                 >  1:2072365‑M1/1‑97 (MQ=255)

  gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                  <  1:659266‑M1/99‑3 (MQ=255)

   tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                                <  2:2327023‑M1/100‑5 (MQ=255)

    ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                               >  2:1359225‑M1/1‑95 (MQ=255)

    ttCATGCCTTCGCTCCTCATCTTACTTAAAAAAAAAA‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                               <  2:3212779‑M1/100‑6 (MQ=255)

     tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                              <  1:2934276‑M1/100‑7 (MQ=255)

      cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                             <  2:3236448‑M1/100‑8 (MQ=255)

      cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                             <  2:3085416‑M1/100‑8 (MQ=255)

       aTGCCTTCGCTCCTCATCTTACTTTTCTACATAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                                                                                            <  1:2479352/100‑1 (MQ=255)

        tGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                           >  1:1210807‑M1/1‑91 (MQ=255)

         gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                                                          <  2:1091327‑M1/100‑11 (MQ=255)

         gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                                                          <  2:769503‑M1/100‑11 (MQ=255)

          ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                         <  2:2107343‑M1/100‑12 (MQ=255)

          ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                         >  2:247482‑M1/1‑89 (MQ=255)

           cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCAACTCATCAGTCGCCTTAAAAATCAGATTGCCAGCGCCGCCTTCGGCcgtaccctgcac                                                                                        >  2:2396227‑M1/1‑88 (MQ=255)

               gCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                                                                    <  2:909372‑M1/100‑17 (MQ=255)

                ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                   <  1:2262729‑M1/100‑18 (MQ=255)

                ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                                   <  1:3032435/100‑1 (MQ=255)

                 tcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                  >  1:1435187‑M1/1‑82 (MQ=255)

                      aTCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgc                                                                             >  2:2218809‑M1/1‑77 (MQ=255)

                        cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                           >  1:383418‑M1/1‑75 (MQ=255)

                                tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                   <  2:237966‑M1/100‑34 (MQ=255)

                                 cTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                                                                  <  2:400736‑M1/100‑35 (MQ=255)

                                  tccagAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                 >  1:2212685‑M1/3‑65 (MQ=255)

                                  tACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                 >  1:1819933‑M1/1‑65 (MQ=255)

                                   aCAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                 >  1:333014‑M1/1‑65 (MQ=255)

                                   aCAGAA‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                 <  1:2607823‑M1/99‑36 (MQ=255)

                                    cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                               <  1:133506‑M1/100‑38 (MQ=255)

                                    cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                               >  1:1194829‑M1/1‑63 (MQ=255)

                                    cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                               <  2:1868785‑M1/100‑38 (MQ=255)

                                       ctgatttttaaggcgactgaTGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                                             >  1:2019320‑M1/10‑61 (MQ=255)

                                       aC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtcccctgcacttcaatgatgctcccgtcttcggtcatc                                                            <  2:1618175‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  1:142147‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  2:429331‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  1:1035443‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  1:1126853‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  1:2766749‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  2:2807414‑M1/100‑41 (MQ=255)

                                       aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                            <  1:2167250‑M1/100‑41 (MQ=255)

                                        c‑AAAAAAAAGGCGACTCATCAGTCGCGTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcgatgatgcgcccgtcttcggtcatca                                                           >  1:1993578‑M1/1‑59 (MQ=255)

                                        c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                           >  2:2267794‑M1/1‑59 (MQ=255)

                                        c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                           >  1:437316‑M1/1‑59 (MQ=255)

                                         gatttttaaggcgactGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                           >  1:1214406‑M1/8‑59 (MQ=255)

                                         gatttttaaggcgactGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccactgcacttcaatgatgcgcccgtcttcggtcatca                                                           >  2:1236548‑M1/8‑59 (MQ=255)

                                           aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                                         <  2:1641243‑M1/100‑44 (MQ=255)

                                           aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAATTTGCCAGCGCCGCATTCTGCcgtaccctgcacttccatgatgcgcccgtcttcggtcatcact                                                         >  2:2984867‑M1/1‑57 (MQ=255)

                                            ttttaaggcgACTGATGAGTCGCCTTAAAAATCAGTTTGCCGGCGCCGCCTTCTGCcgtccactgcacttcaatgatgcgcccgtcttcggtcatcact                                                         >  2:757737‑M1/5‑56 (MQ=255)

                                            aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                        >  2:368417‑M1/1‑56 (MQ=255)

                                              aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                      >  2:413408‑M1/1‑54 (MQ=255)

                                                aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                     >  1:1422371‑M1/1‑52 (MQ=255)

                                                 aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                   <  2:3138971‑M1/100‑50 (MQ=255)

                                                  ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca                                                  <  2:788985‑M1/100‑51 (MQ=255)

                                                  ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca                                                  <  2:1093290‑M1/100‑51 (MQ=255)

                                                    cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                >  1:966323‑M1/1‑48 (MQ=255)

                                                    cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                <  2:47152‑M1/100‑53 (MQ=255)

                                                     gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                               >  2:2558218‑M1/1‑47 (MQ=255)

                                                      aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                              <  1:2080602‑M1/100‑55 (MQ=255)

                                                       cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                             >  2:556928‑M1/1‑45 (MQ=255)

                                                       cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg                                             <  1:2831068‑M1/100‑56 (MQ=255)

                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                            <  2:1581375‑M1/100‑57 (MQ=255)

                                                        tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                            <  1:919839‑M1/100‑57 (MQ=255)

                                                         catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt                                           >  2:3199937‑M1/1‑43 (MQ=255)

                                                            caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                        >  1:1157030‑M1/1‑40 (MQ=255)

                                                            caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                        >  1:2712001‑M1/1‑40 (MQ=255)

                                                             aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctct                                       >  2:1924454‑M1/1‑39 (MQ=255)

                                                              gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctg                                      >  2:543512‑M1/1‑38 (MQ=255)

                                                              gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctg                                      >  1:630484‑M1/1‑38 (MQ=255)

                                                                 gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgc                                     >  2:2743227‑M1/1‑35 (MQ=255)

                                                                   cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                 >  1:1683833‑M1/1‑33 (MQ=255)

                                                                   cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                 >  1:201910‑M1/1‑33 (MQ=255)

                                                                    ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcag                                >  2:2318380‑M1/1‑32 (MQ=255)

                                                                     tAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                               >  1:3244616‑M1/1‑31 (MQ=255)

                                                                      aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagag                              <  2:1187237‑M1/100‑71 (MQ=255)

                                                                         aaTCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                           >  1:1051222‑M1/1‑27 (MQ=255)

                                                                              gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                      >  2:933129‑M1/1‑22 (MQ=255)

                                                                              gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                      <  1:527645‑M1/100‑79 (MQ=255)

                                                                              gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                      >  1:744073‑M1/1‑22 (MQ=255)

                                                                              gTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac                      <  1:867787‑M1/100‑79 (MQ=255)

                                                                               tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                     >  2:343718‑M1/1‑21 (MQ=255)

                                                                               tttGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg                     <  2:2121057‑M1/100‑80 (MQ=255)

                                                                                  gCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtat                  >  1:1874664‑M1/1‑18 (MQ=255)

                                                                                  gCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtat                  >  1:2667072‑M1/1‑18 (MQ=255)

                                                                                   ccAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtatt                 >  2:2450585‑M1/1‑17 (MQ=255)

                                                                                    cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                >  2:2351506‑M1/1‑16 (MQ=255)

                                                                                    cAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctatgcggcagagtcttcaacgtatt                 >  2:624937‑M1/1‑16 (MQ=255)

                                                                                     aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc               <  2:704255‑M1/100‑86 (MQ=255)

                                                                                     aGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc               <  2:2278720‑M1/100‑86 (MQ=255)

                                                                                       cgccgcCTTCTGCcgtcccctgcacttcgatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgacagagtcttcaacgtattcc               <  2:1993572‑M1/98‑86 (MQ=255)

                                                                                         ccgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagat           >  2:356419‑M1/1‑11 (MQ=255)

                                                                                          cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc          <  2:133581‑M1/100‑91 (MQ=255)

                                                                                           gccTTCTGCcgtcccctgcacttcaattatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcg         <  1:2984884‑M1/100‑92 (MQ=255)

                                                                                            ccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgc        <  2:1164595‑M1/100‑93 (MQ=255)

                                                                                              ttCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaa      >  2:322882‑M1/1‑6 (MQ=255)

                                                                                                cTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacttcgttcatgtcggtctctgcggcagagtcttcaccgtattccagatcgcaaac    >  2:67723‑M1/1‑4 (MQ=255)

                                                                                                  gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccg  >  1:1178108‑M1/1‑2 (MQ=255)

                                                                                                  gCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccg  >  2:1086512‑M1/1‑2 (MQ=255)

TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCA  >  NC_000913/3815760‑3815956

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