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	309 [0]	[0] 328	[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)	218 (1.090) +GCATGCGGATCCTCTAGAGTCGACATATG	100/126	0.0	100%	coding (283/825 nt)	rhaD	rhamnulose‑phosphate aldolase
*	?	NC_000913	4095033 =	0 (0.000)	218 (1.090) +GCATGCGGATCCTCTAGAGTCGACATATG	100/126	0.0	100%	intergenic (‑4/‑284)	rhaB/rhaS	rhamnulokinase/RhaS

Only 100 of 328 total aligned reads displayed.

AGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094941‑4095125

agtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgg                                                                                              >  1:927519‑M2/93‑93 (MQ=255)

agtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgg                                                                                              <  2:2381219‑M2/1‑1 (MQ=255)

agtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgg                                                                                              >  2:4524396‑M2/93‑93 (MQ=255)

 gtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                                             >  2:4882254‑M2/92‑93 (MQ=255)

  tggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTg                                                                                            >  2:4271328‑M2/91‑93 (MQ=255)

   ggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGa                                                                                           <  1:6557558‑M2/4‑1 (MQ=255)

   ggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGa                                                                                           <  2:3371486‑M2/4‑1 (MQ=255)

   ggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGa                                                                                           >  1:2448317‑M2/90‑93 (MQ=255)

    gggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAt                                                                                          >  2:4302023‑M2/89‑93 (MQ=255)

     ggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATc                                                                                         <  2:530252‑M2/6‑1 (MQ=255)

      gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  2:2481845‑M2/7‑1 (MQ=255)

      gacggcttcgttagttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        <  1:1017358‑M2/7‑1 (MQ=255)

       acggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                                       <  1:7315651‑M2/8‑1 (MQ=255)

          gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    >  1:465372‑M2/83‑93 (MQ=255)

          gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    >  1:2052855‑M2/83‑93 (MQ=255)

            ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                  >  2:6875132‑M2/81‑93 (MQ=255)

               gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               >  2:2576939‑M2/78‑93 (MQ=255)

                ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              <  1:1721369‑M2/17‑1 (MQ=255)

                 tggttaacccccaaagagtgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTc                                                                             <  1:6137649‑M2/18‑1 (MQ=255)

                  ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            >  2:2848820‑M2/75‑93 (MQ=255)

                  ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  1:4938383‑M2/19‑1 (MQ=255)

                  ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  2:6003758‑M2/19‑1 (MQ=255)

                  ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                            <  1:6929304‑M2/19‑1 (MQ=255)

                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        <  2:7258981‑M2/23‑1 (MQ=255)

                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        >  1:738214‑M2/71‑93 (MQ=255)

                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        <  1:5280173‑M2/23‑1 (MQ=255)

                       acccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                       >  2:768880‑M2/70‑93 (MQ=255)

                        cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      >  2:5029941‑M2/69‑93 (MQ=255)

                         ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                     >  2:2676932‑M2/68‑93 (MQ=255)

                          cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAcc                                                                    >  2:6145448‑M2/67‑93 (MQ=255)

                            caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  <  2:2157208‑M2/29‑1 (MQ=255)

                            caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  >  1:1430461‑M2/65‑93 (MQ=255)

                            caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  <  2:6567198‑M2/29‑1 (MQ=255)

                             aagaatgtggtagccccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGt                                                                 <  1:1197557‑M2/30‑1 (MQ=255)

                              aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                >  2:3857882‑M2/63‑93 (MQ=255)

                              aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTc                                                                <  1:6189044‑M2/31‑1 (MQ=255)

                               agaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                               <  1:2393032‑M2/32‑1 (MQ=255)

                                   tgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaa                                                           <  1:6481207‑M2/36‑1 (MQ=255)

                                    gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaaaa                                                          <  2:4074198‑M2/37‑1 (MQ=255)

                                       gtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                                       <  1:6301583‑M2/40‑1 (MQ=255)

                                          gcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt                                                    >  1:3243328‑M2/51‑93 (MQ=255)

                                            ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACAAGTCTAAAAAGCGCCTGa                                                  >  1:415355‑M2/49‑93 (MQ=255)

                                             cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGaa                                                 <  1:4110547‑M2/46‑1 (MQ=255)

                                               cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                               >  1:1551176‑M2/46‑93 (MQ=255)

                                               cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                               >  2:7356682‑M2/46‑93 (MQ=255)

                                               cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                               <  2:83373‑M2/48‑1 (MQ=255)

                                                gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              >  2:885961‑M2/45‑93 (MQ=255)

                                                gccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTc                                              <  1:6616730‑M2/49‑1 (MQ=255)

                                                 ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                             >  1:1861816‑M2/44‑93 (MQ=255)

                                                 ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                             <  1:5115100‑M2/50‑1 (MQ=255)

                                                    tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                                          >  2:1029082‑M2/41‑93 (MQ=255)

                                                    tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                                          <  2:3126180‑M2/53‑1 (MQ=255)

                                                      gctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAcc                                        >  2:2194599‑M2/39‑93 (MQ=255)

                                                         gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     <  1:1043561‑M2/58‑1 (MQ=255)

                                                          tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGCATTCGCGACCTtct                                    <  2:5191688‑M2/59‑1 (MQ=255)

                                                          tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    >  2:1064306‑M2/35‑93 (MQ=255)

                                                          tcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtct                                    <  1:5411826‑M2/59‑1 (MQ=255)

                                                           cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                   <  1:4468627‑M2/60‑1 (MQ=255)

                                                           cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                   >  2:5874410‑M2/34‑93 (MQ=255)

                                                            gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                  >  1:1214853‑M2/33‑93 (MQ=255)

                                                             acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                                 <  2:1925942‑M2/62‑1 (MQ=255)

                                                             acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                                 <  2:7099478‑M2/62‑1 (MQ=255)

                                                                catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              <  1:2002289‑M2/65‑1 (MQ=255)

                                                                 atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACt                             <  2:3065894‑M2/66‑1 (MQ=255)

                                                                 atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACt                             <  2:1390526‑M2/66‑1 (MQ=255)

                                                                   gcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                           <  2:1113911‑M2/68‑1 (MQ=255)

                                                                     ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         >  1:3689904‑M2/24‑93 (MQ=255)

                                                                     ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         >  1:5592115‑M2/24‑93 (MQ=255)

                                                                      gatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTGAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAg                        >  2:6907709‑M2/23‑93 (MQ=255)

                                                                        tcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                      >  1:4006883‑M2/21‑93 (MQ=255)

                                                                          ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  1:2681028‑M2/19‑93 (MQ=255)

                                                                           tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   >  1:1816764‑M2/18‑93 (MQ=255)

                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  >  1:4324036‑M2/17‑93 (MQ=255)

                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  <  1:5965573‑M2/77‑1 (MQ=255)

                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  <  1:3113342‑M2/77‑1 (MQ=255)

                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 <  2:3926913‑M2/78‑1 (MQ=255)

                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 <  1:382801‑M2/78‑1 (MQ=255)

                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 <  1:3483633‑M2/78‑1 (MQ=255)

                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 >  2:5484993‑M2/16‑93 (MQ=255)

                                                                                 gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGCCAGGAAAATGGGcc             <  2:6445744‑M2/82‑1 (MQ=255)

                                                                                 gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGcc             <  2:7423022‑M2/82‑1 (MQ=255)

                                                                                 gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGcc             >  1:5800122‑M2/12‑93 (MQ=255)

                                                                                 gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGcc             >  1:3425451‑M2/12‑93 (MQ=255)

                                                                                  tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCa            >  1:7028500‑M2/11‑93 (MQ=255)

                                                                                   cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAt           <  1:4738582‑M2/84‑1 (MQ=255)

                                                                                   cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAt           <  2:6624864‑M2/84‑1 (MQ=255)

                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAtt          >  2:3724704‑M2/9‑93 (MQ=255)

                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAtt          >  2:6742566‑M2/9‑93 (MQ=255)

                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAtt          <  2:1690466‑M2/85‑1 (MQ=255)

                                                                                    gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTAACTGACAGGAAAATGGGCCAtt          <  1:6033633‑M2/85‑1 (MQ=255)

                                                                                     acatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTg         <  2:5606865‑M2/86‑1 (MQ=255)

                                                                                      catatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTgg        <  1:6685650‑M2/87‑1 (MQ=255)

                                                                                       atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGc       <  2:1212415‑M2/88‑1 (MQ=255)

                                                                                        tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCa      >  1:5723379‑M2/5‑93 (MQ=255)

                                                                                          tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAc    >  2:4726745‑M2/3‑93 (MQ=255)

                                                                                          tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAc    >  1:7140851‑M2/3‑93 (MQ=255)

                                                                                           gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAcc   <  1:2797356‑M2/92‑1 (MQ=255)

                                                                                           gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAcc   <  1:4695775‑M2/92‑1 (MQ=255)

                                                                                           gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAcc   >  2:2070154‑M2/2‑93 (MQ=255)

                                                                                            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCa  <  2:6365388/93‑1 (MQ=255)

AGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094941‑4095125

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