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	585 [0]	[0] 589	[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)	384 (0.960) +GCATGCGGATCCTCTAGAGTCGACATATG	113/132	0.7	100%	coding (283/825 nt)	rhaD	rhamnulose‑phosphate aldolase
*	?	NC_000913	4095033 =	0 (0.000)	384 (0.960) +GCATGCGGATCCTCTAGAGTCGACATATG	113/132	0.7	100%	intergenic (‑4/‑284)	rhaB/rhaS	rhamnulokinase/RhaS

Only 100 of 589 total aligned reads displayed.

TCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGG  >  NC_000913/4094939‑4095128

gaagtgggggcggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                                                >  1:12971089‑M2/95‑96 (MQ=255)

gaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                                                >  1:1843010‑M2/95‑96 (MQ=255)

gaagtggggacggcttcgttggttaacccccaaagaatgtggcagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGt                                                                                                <  2:628978‑M2/2‑1 (MQ=255)

 aagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatggcg                                                                                               <  1:11862681‑M2/3‑3 (MQ=255)

 aagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTg                                                                                               <  2:1263104‑M2/3‑1 (MQ=255)

   gtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAt                                                                                             >  1:8007470‑M2/92‑96 (MQ=255)

      gggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCct                                                                                          <  1:5156698‑M2/8‑1 (MQ=255)

       ggacggcttccttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctg                                                                                         <  2:13177841‑M2/9‑1 (MQ=255)

        gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGAtccttc                                                                                        <  1:10675499‑M2/10‑3 (MQ=255)

         acggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgtcctatgGTGATCctgct                                                                                       >  1:2865071‑M2/86‑96 (MQ=255)

          cggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgctg                                                                                      <  2:12413973‑M2/12‑1 (MQ=255)

           ggcttcgttggttaacccccaaagaatggggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                     >  1:2193867‑M2/84‑96 (MQ=255)

           agcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGa                                                                                     >  1:13840956‑M2/84‑96 (MQ=255)

             cttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagccgacatatgGTGATCCTGCTGAAt                                                                                   >  2:7788552‑M2/82‑96 (MQ=255)

              ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGTAtt                                                                                  >  2:1015959‑M2/81‑96 (MQ=255)

              ttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                                  >  1:5543543‑M2/81‑96 (MQ=255)

                cgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgccatatgGTGATCCTGCTGAATTTc                                                                                <  2:3632234‑M2/18‑1 (MQ=255)

                cgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTc                                                                                <  2:1667869‑M2/18‑1 (MQ=255)

                 gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCa                                                                               >  2:8620097‑M2/78‑96 (MQ=255)

                    ggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTa                                                                            <  2:2729000‑M2/22‑1 (MQ=255)

                      ttaaccccccaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                          <  1:4232836‑M2/24‑1 (MQ=255)

                      ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                          <  2:11162955‑M2/24‑1 (MQ=255)

                      ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACg                                                                          >  2:7549245‑M2/73‑96 (MQ=255)

                        aacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGAc                                                                        >  1:11022859‑M2/71‑96 (MQ=255)

                           ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                     <  1:8548120‑M2/29‑1 (MQ=255)

                           ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                     <  1:2932637‑M2/29‑1 (MQ=255)

                           ccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                     <  2:5251254‑M2/29‑1 (MQ=255)

                              caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                                  <  1:2369637‑M2/32‑1 (MQ=255)

                              caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                                  <  2:11998523‑M2/32‑1 (MQ=255)

                                aagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTaa                                                                <  2:11574904‑M2/34‑1 (MQ=255)

                                     tgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgc                                                           <  2:4008766‑M2/39‑1 (MQ=255)

                                      gtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                                          >  2:6693582‑M2/57‑96 (MQ=255)

                                      gtgggagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgcg                                                          >  1:13987422‑M2/57‑96 (MQ=255)

                                           agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                                     >  2:8762893‑M2/52‑96 (MQ=255)

                                             cccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                   <  1:5909596‑M2/47‑1 (MQ=255)

                                             cccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAt                                                   >  2:7308719‑M2/50‑96 (MQ=255)

                                              ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAAtt                                                  >  1:141864‑M2/49‑96 (MQ=255)

                                                gcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcg                                                >  2:4699885‑M2/47‑96 (MQ=255)

                                                 cgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgc                                               <  2:12185216‑M2/51‑1 (MQ=255)

                                                   ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                                             >  2:1346276‑M2/44‑96 (MQ=255)

                                                   ccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGa                                             >  1:9021480‑M2/44‑96 (MQ=255)

                                                    agtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc                                            >  1:11423547‑M2/43‑96 (MQ=255)

                                                      tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCt                                          >  2:6028553‑M2/41‑96 (MQ=255)

                                                      tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCt                                          <  1:944586‑M2/56‑1 (MQ=255)

                                                      tcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCCGAATTCGCGACCt                                          >  2:4991495‑M2/41‑96 (MQ=255)

                                                       cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCtt                                         >  1:1029855‑M2/40‑96 (MQ=255)

                                                          tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                      <  2:966988‑M2/60‑1 (MQ=255)

                                                          tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                      >  2:9028613‑M2/37‑96 (MQ=255)

                                                          tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                      >  1:13435263‑M2/37‑96 (MQ=255)

                                                          tgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtctc                                      <  1:1626608‑M2/60‑1 (MQ=255)

                                                           gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                     >  2:9920427‑M2/36‑96 (MQ=255)

                                                           gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                     >  1:8225224‑M2/36‑96 (MQ=255)

                                                           gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCg                                     >  2:12927769‑M2/36‑96 (MQ=255)

                                                           gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAACTTCTCg                                     >  1:9847716‑M2/36‑96 (MQ=255)

                                                             cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                   <  1:7611101‑M2/63‑1 (MQ=255)

                                                             cgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGtt                                   <  2:11773799‑M2/63‑1 (MQ=255)

                                                              gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTa                                  <  2:2206104‑M2/64‑1 (MQ=255)

                                                              gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTa                                  >  1:13253341‑M2/33‑96 (MQ=255)

                                                              gacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTa                                  >  1:7137356‑M2/33‑96 (MQ=255)

                                                                  catgcggctcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                              <  1:6251588‑M2/68‑1 (MQ=255)

                                                                  catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                              >  2:2531285‑M2/29‑96 (MQ=255)

                                                                  catgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                              >  1:11566416‑M2/29‑96 (MQ=255)

                                                                   atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGAc                             <  1:2071965‑M2/69‑1 (MQ=255)

                                                                   atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGAc                             <  1:6830202‑M2/69‑1 (MQ=255)

                                                                   atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGAc                             >  1:326292‑M2/28‑96 (MQ=255)

                                                                   atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGAc                             >  1:3543411‑M2/28‑96 (MQ=255)

                                                                   atgcggatcctctagaggcgacatatgGGGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGAc                             >  1:4508074‑M2/28‑96 (MQ=255)

                                                                    tgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACa                            <  2:11443531‑M2/70‑1 (MQ=255)

                                                                     gcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAg                           >  1:1220013‑M2/26‑96 (MQ=255)

                                                                     gcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAg                           >  2:4314629‑M2/26‑96 (MQ=255)

                                                                       ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                         >  2:9242033‑M2/24‑96 (MQ=255)

                                                                       ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                         >  1:1239532‑M2/24‑96 (MQ=255)

                                                                       ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                         >  1:6405141‑M2/24‑96 (MQ=255)

                                                                         atcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                       >  1:7802505‑M2/22‑96 (MQ=255)

                                                                          tcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                      >  2:10292861‑M2/21‑96 (MQ=255)

                                                                           cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                     <  2:14218937‑M2/77‑1 (MQ=255)

                                                                           cctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAt                     <  1:1450902‑M2/77‑1 (MQ=255)

                                                                            ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATg                    >  1:9660595‑M2/19‑96 (MQ=255)

                                                                             tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATgg                   >  1:2556027‑M2/18‑96 (MQ=255)

                                                                             tctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTGCTCGTTACTGACAGGAAaatag                   >  1:4916074‑M2/18‑94 (MQ=255)

                                                                              ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATggg                  >  1:9207971‑M2/17‑96 (MQ=255)

                                                                              ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATggg                  >  1:8730382‑M2/17‑96 (MQ=255)

                                                                                agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGcc                >  2:1276039‑M2/15‑96 (MQ=255)

                                                                                agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGcc                <  2:1874256‑M2/82‑1 (MQ=255)

                                                                                  agtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGggccct              <  1:7548527‑M2/84‑3 (MQ=255)

                                                                                  agtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAt              <  2:4459853‑M2/84‑1 (MQ=255)

                                                                                   gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCAtt             <  2:2946130‑M2/85‑1 (MQ=255)

                                                                                   gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTAACAGGAAAATGGGCCAtt             <  1:12641566‑M2/85‑1 (MQ=255)

                                                                                     cgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTgg           <  1:12520449‑M2/87‑1 (MQ=255)

                                                                                      gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAGTGACAGGAAAATGGGCCATTGGc          >  1:10071516‑M2/9‑96 (MQ=255)

                                                                                      gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGc          >  2:6470475‑M2/9‑96 (MQ=255)

                                                                                       acatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCa         >  2:3254912‑M2/8‑96 (MQ=255)

                                                                                         atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAc       >  2:10713998‑M2/6‑96 (MQ=255)

                                                                                          tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAAcc      <  2:8211694‑M2/92‑1 (MQ=255)

                                                                                           atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCGATGTGCAACCa     >  1:379683‑M2/4‑96 (MQ=255)

                                                                                           atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCa     <  2:1037554‑M2/93‑1 (MQ=255)

                                                                                           atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAATGACAGGAAAATGGGCCATTGGCAACCa     <  2:13917747‑M2/93‑1 (MQ=255)

                                                                                            tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAg    >  1:4136879‑M2/3‑96 (MQ=255)

                                                                                             gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAgg   <  2:7037412‑M2/95‑1 (MQ=255)

                                                                                              gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAggg  >  1:8461306/1‑96 (MQ=255)

TCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGG  >  NC_000913/4094939‑4095128

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