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	244 [0]	[0] 262	[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)	168 (0.990) +GCATGCGGATCCTCTAGAGTCGACATATG	86/126	0.3	100%	coding (283/825 nt)	rhaD	rhamnulose‑phosphate aldolase
*	?	NC_000913	4095033 =	0 (0.000)	168 (0.990) +GCATGCGGATCCTCTAGAGTCGACATATG	86/126	0.3	100%	intergenic (‑4/‑284)	rhaB/rhaS	rhamnulokinase/RhaS

Only 100 of 262 total aligned reads displayed.

AGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094941‑4095125

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

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

 gtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                                             <  2:526601‑M2/2‑1 (MQ=255)

     ggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagaggcgacatatgGTGATc                                                                                         >  1:6461687‑M2/88‑93 (MQ=255)

      gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  2:1105848‑M2/87‑93 (MQ=255)

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

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

        cggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctg                                                                                      >  2:4655900‑M2/85‑93 (MQ=255)

         ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     >  1:225666‑M2/84‑93 (MQ=255)

          gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    >  2:3387585‑M2/83‑93 (MQ=255)

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

          gcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                    <  2:3224414‑M2/11‑1 (MQ=255)

            ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                  >  1:2835229‑M2/81‑93 (MQ=255)

            ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                  >  1:3364048‑M2/81‑93 (MQ=255)

            ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                  >  1:2336658‑M2/81‑93 (MQ=255)

             tcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 <  1:3859261‑M2/14‑1 (MQ=255)

              cgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAt                                                                                <  1:2399655‑M2/15‑1 (MQ=255)

              cgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAt                                                                                <  1:4563748‑M2/15‑1 (MQ=255)

               gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               <  1:4802057‑M2/16‑1 (MQ=255)

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

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

               gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                               >  1:3549169‑M2/78‑93 (MQ=255)

                ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              >  2:633793‑M2/77‑93 (MQ=255)

                   gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           >  1:3621456‑M2/74‑93 (MQ=255)

                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        >  2:6005434‑M2/71‑93 (MQ=255)

                      aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                        >  2:4593416‑M2/71‑93 (MQ=255)

                        cccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGa                                                                      >  1:3508643‑M2/69‑93 (MQ=255)

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

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

                          cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAcc                                                                    <  1:1595479‑M2/27‑1 (MQ=255)

                             aaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCagc                                                                 >  1:341261‑M2/64‑92 (MQ=255)

                             aaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGt                                                                 <  1:3427604‑M2/30‑1 (MQ=255)

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

                               agaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                               >  2:5117252‑M2/62‑93 (MQ=255)

                               agaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                               >  1:1034957‑M2/62‑93 (MQ=255)

                                 aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaa                                                             >  2:687508‑M2/60‑93 (MQ=255)

                                  atgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaaa                                                            <  2:3610240‑M2/35‑1 (MQ=255)

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

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

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

                                         agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc                                                     >  2:118804‑M2/52‑93 (MQ=255)

                                         agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGcc                                                     <  2:490258‑M2/42‑1 (MQ=255)

                                          gcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt                                                    <  1:179493‑M2/43‑1 (MQ=255)

                                          gcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCt                                                    >  2:5749649‑M2/51‑93 (MQ=255)

                                            ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                                  <  1:212978‑M2/45‑1 (MQ=255)

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

                                              gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                <  1:5491015‑M2/47‑1 (MQ=255)

                                              gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                >  1:2066134‑M2/47‑93 (MQ=255)

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

                                                   gtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg                                           >  2:1869506‑M2/42‑93 (MQ=255)

                                                       ctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCt                                       <  1:5866189‑M2/56‑1 (MQ=255)

                                                       ctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCt                                       <  1:3146818‑M2/56‑1 (MQ=255)

                                                       ctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCt                                       >  2:721893‑M2/38‑93 (MQ=255)

                                                        tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                      >  1:5541802‑M2/37‑93 (MQ=255)

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

                                                            gacgcatgcggatcctctagagtcgacatatcGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                  <  2:2498154‑M2/61‑1 (MQ=255)

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

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

                                                             acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                                 >  1:1539071‑M2/32‑93 (MQ=255)

                                                              cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                >  2:1757482‑M2/31‑93 (MQ=255)

                                                              cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                >  1:323769‑M2/31‑93 (MQ=255)

                                                              cgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                >  2:6372756‑M2/31‑93 (MQ=255)

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

                                                                catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                              >  1:5215993‑M2/29‑93 (MQ=255)

                                                                  tgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTg                            >  2:2940342‑M2/27‑93 (MQ=255)

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

                                                                     ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         >  2:5424958‑M2/24‑93 (MQ=255)

                                                                     ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                         <  1:4998686‑M2/70‑1 (MQ=255)

                                                                      gatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAg                        >  2:3639172‑M2/23‑93 (MQ=255)

                                                                       atcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAgg                       >  1:43123‑M2/22‑93 (MQ=255)

                                                                       atcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAgg                       >  1:534742‑M2/22‑93 (MQ=255)

                                                                        tcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTATTGACAGGa                      <  1:6263539‑M2/73‑1 (MQ=255)

                                                                        tcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                      <  1:6355880‑M2/73‑1 (MQ=255)

                                                                           tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   >  2:611019‑M2/18‑93 (MQ=255)

                                                                           tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   >  2:3078259‑M2/18‑93 (MQ=255)

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

                                                                           tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                   <  2:2352744‑M2/76‑1 (MQ=255)

                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  >  2:5801067‑M2/17‑93 (MQ=255)

                                                                            ctagagtcgacatatgGTGATCCTGCTGAATTTCATCACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                  <  2:2539190‑M2/77‑1 (MQ=255)

                                                                             tagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAACAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                 >  1:3672172‑M2/16‑93 (MQ=255)

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

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

                                                                              agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATgg                >  2:2430201‑M2/15‑93 (MQ=255)

                                                                              agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATgg                >  1:5905191‑M2/15‑93 (MQ=255)

                                                                               gagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATggg               >  1:1756309‑M2/14‑93 (MQ=255)

                                                                                agtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGc              <  1:4495032‑M2/81‑1 (MQ=255)

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

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

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

                                                                                  tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCa            <  1:3095014‑M2/83‑1 (MQ=255)

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

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

                                                                                         atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCaa     <  1:2742033‑M2/90‑1 (MQ=255)

                                                                                          tgGTGATCCTGCTGAATTACATTACGACCAGTCTAAAAACCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAc    <  1:97839‑M2/91‑1 (MQ=255)

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

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

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

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

                                                                                            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCa  >  2:4107177/1‑93 (MQ=255)

                                                                                            gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCa  >  2:5443941/1‑93 (MQ=255)

AGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094941‑4095125

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