BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,091,572	3461 bp→GCATGCGGATCCTCTAGAGTCGACATATG		[rhaD]–rhaB	[rhaD], rhaA, rhaB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4091572	4095032	3461	232 [0]	[0] 236	[rhaD]–rhaB	[rhaD],rhaA,rhaB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4091571	0 (0.000)	164 (1.320) +GCATGCGGATCCTCTAGAGTCGACATATG	86/126	0.1	100%	coding (283/825 nt)	rhaD	rhamnulose‑phosphate aldolase
*	?	NC_000913	4095033 =	0 (0.000)	164 (1.320) +GCATGCGGATCCTCTAGAGTCGACATATG	86/126	0.1	100%	intergenic (‑4/‑284)	rhaB/rhaS	rhamnulokinase/RhaS

Only 100 of 232 total aligned reads displayed.

GTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCA  >  NC_000913/4091479‑4091660

  gggAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgc                                                                                         <  2:1608657‑M1/93‑3 (MQ=255)

      aaGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                                                                                     <  2:4355263‑M1/93‑7 (MQ=255)

       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                    >  2:4078339‑M1/1‑86 (MQ=255)

       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                    >  1:2232504‑M1/1‑86 (MQ=255)

       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                    <  2:1178832‑M1/93‑8 (MQ=255)

       aGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcg                                                                                    >  2:4560262‑M1/1‑86 (MQ=255)

         gAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcgga                                                                                  >  2:2299961‑M1/1‑84 (MQ=255)

         gAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcgga                                                                                  <  2:2170360‑M1/93‑10 (MQ=255)

            gTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcc                                                                               >  1:1771276‑M1/1‑81 (MQ=255)

              gAGCGGGAAGTTCGGAAGTGGGGACGTCTTCGTTAGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctc                                                                             <  2:3539330‑M1/93‑15 (MQ=255)

               aGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctct                                                                            >  1:380513‑M1/1‑78 (MQ=255)

                  cGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctaga                                                                         <  2:2884305‑M1/93‑19 (MQ=255)

                  cGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctaga                                                                         >  1:4553978‑M1/1‑75 (MQ=255)

                   ggAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagag                                                                        >  1:4217285‑M1/1‑74 (MQ=255)

                   ggAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagag                                                                        <  1:4477083‑M1/93‑20 (MQ=255)

                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAACGAATGTGGTAGCCCGCGCAGTCGCTGTCGACgcatgcggatcctctagagt                                                                       >  1:2918876‑M1/1‑73 (MQ=255)

                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                       >  1:3220861‑M1/1‑73 (MQ=255)

                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                       <  2:935021‑M1/93‑21 (MQ=255)

                    gAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagt                                                                       >  2:3818586‑M1/1‑73 (MQ=255)

                      aGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcg                                                                     >  2:556705‑M1/1‑71 (MQ=255)

                       gTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcga                                                                    <  1:1599200‑M1/93‑24 (MQ=255)

                        ttCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgac                                                                   <  2:4532314‑M1/93‑25 (MQ=255)

                         tCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgaca                                                                  <  1:4249312‑M1/93‑26 (MQ=255)

                          cGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacat                                                                 <  2:3169855‑M1/93‑27 (MQ=255)

                          cGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacat                                                                 <  2:4464028‑M1/93‑27 (MQ=255)

                           ggAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacata                                                                <  2:2071036‑M1/93‑28 (MQ=255)

                           ggAAGTGGGCACGGCTTCGTTGGTTACCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacata                                                                <  1:2302681‑M1/93‑28 (MQ=255)

                            gAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatat                                                               <  2:2526160‑M1/93‑29 (MQ=255)

                             aaGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatg                                                              >  1:2631841‑M1/1‑64 (MQ=255)

                              aGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatgg                                                             <  1:936855‑M1/93‑31 (MQ=255)

                               gTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggt                                                            <  1:2085722‑M1/93‑32 (MQ=255)

                                   ggACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagaggcgacatatggggatc                                                        >  1:3259445‑M1/1‑58 (MQ=255)

                                    gACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcc                                                       >  2:472445‑M1/1‑57 (MQ=255)

                                     aCGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcct                                                      >  1:754789‑M1/1‑56 (MQ=255)

                                      cGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctg                                                     >  1:537804‑M1/1‑55 (MQ=255)

                                       ggCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgc                                                    <  2:684028‑M1/93‑40 (MQ=255)

                                       ggCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgc                                                    <  1:2060320‑M1/93‑40 (MQ=255)

                                       ggCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgc                                                    >  2:3397979‑M1/1‑54 (MQ=255)

                                       ggCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagccgacatatggtgatcctgc                                                    >  1:279753‑M1/1‑54 (MQ=255)

                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagcgacgacctacggtgatcctgct                                                   >  1:2246501‑M1/1‑53 (MQ=255)

                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                   <  1:116021‑M1/93‑41 (MQ=255)

                                        gCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                   >  1:3038736‑M1/1‑53 (MQ=255)

                                         cTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctg                                                  >  1:1334794‑M1/1‑52 (MQ=255)

                                          ttCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctga                                                 >  2:3682931‑M1/1‑51 (MQ=255)

                                           tCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaa                                                >  1:3738523‑M1/1‑50 (MQ=255)

                                            cGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaat                                               <  2:595464‑M1/93‑45 (MQ=255)

                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctattgtcgacatatggtgatcctgctgaatt                                              >  1:4683414‑M1/1‑48 (MQ=255)

                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                              >  1:1543900‑M1/1‑48 (MQ=255)

                                             gTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                              >  1:2784401‑M1/1‑48 (MQ=255)

                                              ttGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaattt                                             >  1:2885792‑M1/1‑47 (MQ=255)

                                              ttGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaattt                                             >  1:2137536‑M1/1‑47 (MQ=255)

                                               tGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttc                                            >  1:2923760‑M1/1‑46 (MQ=255)

                                               tGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttc                                            >  2:4254815‑M1/1‑46 (MQ=255)

                                                 gTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcat                                          <  2:1022283‑M1/93‑50 (MQ=255)

                                                  ttAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatt                                         <  2:3841059‑M1/93‑51 (MQ=255)

                                                  ttAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgccgaatttcatt                                         >  2:2976406‑M1/1‑43 (MQ=255)

                                                   tAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatta                                        <  2:2791891‑M1/93‑52 (MQ=255)

                                                   tAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatta                                        <  2:2023994‑M1/93‑52 (MQ=255)

                                                   tAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcatta                                        <  1:4456769‑M1/93‑52 (MQ=255)

                                                     aCCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacg                                      >  2:1374702‑M1/1‑40 (MQ=255)

                                                      cccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacga                                     >  2:2091705‑M1/1‑39 (MQ=255)

                                                       ccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgac                                    <  2:3292129‑M1/93‑56 (MQ=255)

                                                       ccccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgac                                    <  1:1364682‑M1/93‑56 (MQ=255)

                                                        cccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacc                                   <  2:170723‑M1/93‑57 (MQ=255)

                                                         ccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggctcctctagaatcgacatatggtgatcctggtgaatttcattacgacca                                  >  2:2045505‑M1/1‑36 (MQ=255)

                                                         ccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacca                                  <  2:2717552‑M1/93‑58 (MQ=255)

                                                         ccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgacca                                  <  2:3766252‑M1/93‑58 (MQ=255)

                                                          cAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccag                                 >  2:1419709‑M1/1‑35 (MQ=255)

                                                          cAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccag                                 <  2:2266976‑M1/93‑59 (MQ=255)

                                                             aGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtct                              >  1:901162‑M1/1‑32 (MQ=255)

                                                              gAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtcta                             <  1:857253‑M1/93‑63 (MQ=255)

                                                                aTGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaa                           <  1:3423484‑M1/93‑65 (MQ=255)

                                                                aTGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaa                           <  2:3989786‑M1/93‑65 (MQ=255)

                                                                 tgtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaa                          <  1:957139‑M1/93‑66 (MQ=255)

                                                                 tgtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaa                          >  1:1887692‑M1/1‑28 (MQ=255)

                                                                 tgtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaa                          >  2:1145484‑M1/1‑28 (MQ=255)

                                                                  gtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaa                         <  1:668097‑M1/93‑67 (MQ=255)

                                                                  gtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaa                         <  2:80432‑M1/93‑67 (MQ=255)

                                                                  gtgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaa                         <  2:4163773‑M1/93‑67 (MQ=255)

                                                                   tgGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaag                        >  2:2871926‑M1/1‑26 (MQ=255)

                                                                      tAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagagc                     >  1:258087‑M1/1‑23 (MQ=255)

                                                                       aGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcc                    >  1:2913753‑M1/1‑22 (MQ=255)

                                                                         cccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctg                  >  2:1406763‑M1/1‑20 (MQ=255)

                                                                          ccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctga                 >  1:1483677‑M1/1‑19 (MQ=255)

                                                                          ccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctga                 >  1:4104495‑M1/1‑19 (MQ=255)

                                                                          ccGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctga                 <  1:23644‑M1/93‑75 (MQ=255)

                                                                           cgcgCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaa                <  1:3625477‑M1/93‑76 (MQ=255)

                                                                              gcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc             <  1:1054084‑M1/93‑79 (MQ=255)

                                                                              gcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc             <  2:1934514‑M1/93‑79 (MQ=255)

                                                                              gcCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattc             >  1:4508942‑M1/1‑15 (MQ=255)

                                                                               ccGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcg            <  2:1790839‑M1/93‑80 (MQ=255)

                                                                               ccGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcg            >  1:4359898‑M1/1‑14 (MQ=255)

                                                                                  tCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcga         <  2:3941850‑M1/93‑83 (MQ=255)

                                                                                    gCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacc       >  2:1212392‑M1/1‑9 (MQ=255)

                                                                                     cTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacct      >  2:3417623‑M1/1‑8 (MQ=255)

                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcgtgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt     >  1:2512873‑M1/1‑7 (MQ=255)

                                                                                      tGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgacctt     >  2:877636‑M1/1‑7 (MQ=255)

                                                                                        tCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttct   >  1:3927332‑M1/1‑5 (MQ=255)

                                                                                         cGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctc  <  2:1561996‑M1/93‑90 (MQ=255)

                                                                                         cGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctc  <  2:1289877‑M1/93‑90 (MQ=255)

GTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCA  >  NC_000913/4091479‑4091660

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