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 236 total aligned reads displayed.

GGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094942‑4095125

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

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

     gacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATcc                                                                                        >  1:92619‑M2/87‑93 (MQ=255)

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

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

        ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     >  2:3141262‑M2/84‑93 (MQ=255)

        ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     <  2:684028‑M2/10‑1 (MQ=255)

        ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     >  2:3397979‑M2/84‑93 (MQ=255)

        ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgc                                                                                     <  1:2060320‑M2/10‑1 (MQ=255)

        ggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagccgacatatgGTGATCctgc                                                                                     >  1:279753‑M2/84‑93 (MQ=255)

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

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

            tcgttggttaccccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 <  2:2143384‑M2/14‑1 (MQ=255)

            tcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 >  1:4255011‑M2/80‑93 (MQ=255)

            tcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGaa                                                                                 >  1:3738523‑M2/80‑93 (MQ=255)

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

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

              gttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctattgtcgacatatgGTGATCCTGCTGAAtt                                                                               >  1:4683414‑M2/78‑93 (MQ=255)

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

               ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              >  1:2137536‑M2/77‑93 (MQ=255)

               ttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAttt                                                                              >  1:2885792‑M2/77‑93 (MQ=255)

                tggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTc                                                                             >  2:4254815‑M2/76‑93 (MQ=255)

                tggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTc                                                                             <  2:1840251‑M2/18‑1 (MQ=255)

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

                  gttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAt                                                                           >  2:1578273‑M2/74‑93 (MQ=255)

                   ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                                          <  2:3841059‑M2/21‑1 (MQ=255)

                   ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCAtt                                                                          >  2:377278‑M2/73‑93 (MQ=255)

                   ttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCCGAATTTCAtt                                                                          >  2:2976406‑M2/73‑93 (MQ=255)

                    taacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTa                                                                         <  1:4456769‑M2/22‑1 (MQ=255)

                    taacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTa                                                                         <  1:740326‑M2/22‑1 (MQ=255)

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

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

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

                          ccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggctcctctagaatcgacatatgGTGATCCTGGTGAATTTCATTACGACCa                                                                   >  2:2045505‑M2/66‑93 (MQ=255)

                          ccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                                   <  2:2717552‑M2/28‑1 (MQ=255)

                          ccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCa                                                                   >  2:2577150‑M2/66‑93 (MQ=255)

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

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

                           caaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                  >  2:1419709‑M2/65‑93 (MQ=255)

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

                              agaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGCATTTCATTACGACCAGTCt                                                               >  1:2968082‑M2/62‑93 (MQ=255)

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

                               gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                                              >  2:479615‑M2/61‑93 (MQ=255)

                               gaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTa                                                              <  1:857253‑M2/33‑1 (MQ=255)

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

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

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

                                    tggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAg                                                         >  1:3248706‑M2/56‑93 (MQ=255)

                                    tggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAg                                                         >  2:2871926‑M2/56‑93 (MQ=255)

                                     ggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAgc                                                        >  2:1734705‑M2/55‑93 (MQ=255)

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

                                          cccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTg                                                   >  1:3592557‑M2/50‑93 (MQ=255)

                                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                                  >  2:1509816‑M2/49‑93 (MQ=255)

                                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                                  >  1:4104495‑M2/49‑93 (MQ=255)

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

                                           ccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGa                                                  >  1:1090527‑M2/49‑93 (MQ=255)

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

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

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

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

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

                                                  gtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTcgcg                                           >  1:532203‑M2/42‑93 (MQ=255)

                                                    cgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAc                                         <  2:2041631‑M2/54‑1 (MQ=255)

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

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

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

                                                        gtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTtc                                     >  1:745938‑M2/36‑93 (MQ=255)

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

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

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

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

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

                                                              gcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTa                               <  2:675728‑M2/64‑1 (MQ=255)

                                                              gcatgcgcatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTa                               >  1:190547‑M2/30‑93 (MQ=255)

                                                                atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACt                             >  2:903992‑M2/28‑93 (MQ=255)

                                                                atgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACt                             >  1:4165529‑M2/28‑93 (MQ=255)

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

                                                                  gcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGa                           >  1:3211643‑M2/26‑93 (MQ=255)

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

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

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

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

                                                                     gatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAg                        <  1:3050926‑M2/71‑1 (MQ=255)

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

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

                                                                         ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  2:794209‑M2/19‑93 (MQ=255)

                                                                         ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  2:1023560‑M2/19‑93 (MQ=255)

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

                                                                         ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaa                    >  2:3523646‑M2/19‑93 (MQ=255)

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

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

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

                                                                             agagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAAtga                >  2:3940972‑M2/15‑92 (MQ=255)

                                                                                 tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa            >  1:3948985‑M2/11‑93 (MQ=255)

                                                                                   gacatatgGTGATGCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt          >  2:4502351‑M2/9‑93 (MQ=255)

                                                                                      atatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGc       >  1:222791‑M2/6‑93 (MQ=255)

                                                                                      atatgGAGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGTGCCTGAATTCGCGACCTTCTCGTTACTGGCAGGAAAATGAGCCATTGGc       >  2:734235‑M2/6‑93 (MQ=255)

                                                                                       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa      >  2:3609400‑M2/5‑93 (MQ=255)

                                                                                          gGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAcc   <  1:2899489‑M2/92‑1 (MQ=255)

                                                                                           gtgATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCa  >  2:1616416/1‑93 (MQ=255)

GGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCA  >  NC_000913/4094942‑4095125

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