BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	66,533	3540 bp→CCCGAGACCCGCACGCGTCGACCTGCA		[araD]–araB	[araD], araA, araB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	66533	70072	3540	157 [0]	[0] 175	[araD]–araB	[araD],araA,araB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 66532	0 (0.000)	116 (0.900) +CCCGAGACCCGCACGCGTCGACCTGCA	70/130	0.3	100%	coding (19/696 nt)	araD	L‑ribulose‑5‑phosphate 4‑epimerase
*	?	NC_000913	70073 =	0 (0.000)	116 (0.900) +CCCGAGACCCGCACGCGTCGACCTGCA	70/130	0.3	100%	intergenic (‑25/‑314)	araB/araC	L‑ribulokinase/transcriptional regulator for ara operon

Only 100 of 175 total aligned reads displayed.

GTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGA  >  NC_000913/69981‑70165

ttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcag                                                                                              >  2:1707653‑M2/93‑93 (MQ=255)

ttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcag                                                                                              >  2:3526886‑M2/93‑93 (MQ=255)

 tgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcagg                                                                                             <  1:2413367‑M2/2‑1 (MQ=255)

 tgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcagg                                                                                             <  2:3958182‑M2/2‑1 (MQ=255)

  gccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGt                                                                                            >  1:3833476‑M2/91‑93 (MQ=255)

  gccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGt                                                                                            <  2:3790694‑M2/3‑1 (MQ=255)

   ccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTa                                                                                           <  1:2018120‑M2/4‑1 (MQ=255)

   ccccatgtgaacgtgaccaggttgtcttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTa                                                                                           <  1:4253928‑M2/4‑1 (MQ=255)

    cccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTAt                                                                                          <  1:3917354‑M2/5‑1 (MQ=255)

     ccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATg                                                                                         <  2:3863219‑M2/6‑1 (MQ=255)

     ccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATg                                                                                         >  2:1283849‑M2/88‑93 (MQ=255)

     ccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATg                                                                                         >  1:2927755‑M2/88‑93 (MQ=255)

      catgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATgg                                                                                        >  1:2848265‑M2/87‑93 (MQ=255)

        tgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGgag                                                                                      >  2:4207929‑M2/85‑93 (MQ=255)

         gtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGgaga                                                                                     <  1:279154‑M2/10‑1 (MQ=255)

          tgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGaa                                                                                    >  1:981084‑M2/83‑93 (MQ=255)

            agcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAAc                                                                                  <  1:1900626‑M2/13‑1 (MQ=255)

             gcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACa                                                                                 <  1:2665450‑M2/14‑1 (MQ=255)

             gcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACa                                                                                 <  2:3416580‑M2/14‑1 (MQ=255)

                tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgctcgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                              <  1:3295325‑M2/17‑1 (MQ=255)

                tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                              <  1:3689912‑M2/17‑1 (MQ=255)

                tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                              <  1:1140386‑M2/17‑1 (MQ=255)

                 gaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTag                                                                             <  2:2814140‑M2/18‑1 (MQ=255)

                  accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            <  1:2706959‑M2/19‑1 (MQ=255)

                  accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            <  2:4545447‑M2/19‑1 (MQ=255)

                   ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                           >  2:247435‑M2/74‑93 (MQ=255)

                    caggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagaga                                                                          >  1:3204834‑M2/73‑93 (MQ=255)

                    caggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagaga                                                                          >  2:65413‑M2/73‑93 (MQ=255)

                    caggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATAGAGAAACAGTagaga                                                                          <  2:143914‑M2/21‑1 (MQ=255)

                     aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagagag                                                                         >  1:2328814‑M2/72‑93 (MQ=255)

                     aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagagag                                                                         >  1:3389003‑M2/72‑93 (MQ=255)

                     aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagagag                                                                         >  1:4679147‑M2/72‑93 (MQ=255)

                     aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagagag                                                                         >  2:483507‑M2/72‑93 (MQ=255)

                      ggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGt                                                                        >  1:955026‑M2/71‑93 (MQ=255)

                            gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGAATGGAGAAACAGTAGAGAGGTGAGAt                                                                  >  1:194869‑M2/65‑93 (MQ=255)

                             gtttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATa                                                                 <  1:1103275‑M2/30‑1 (MQ=255)

                              tttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaa                                                                >  1:245574‑M2/63‑93 (MQ=255)

                               ttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaa                                                               >  2:1754990‑M2/62‑93 (MQ=255)

                                tggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaa                                                              <  1:3642414‑M2/33‑1 (MQ=255)

                                tggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGGTGCGATaaaa                                                              >  2:1755175‑M2/61‑93 (MQ=255)

                                  gcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAg                                                            <  1:2094040‑M2/35‑1 (MQ=255)

                                     gcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGt                                                         <  2:1226166‑M2/38‑1 (MQ=255)

                                     gcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGt                                                         <  1:697499‑M2/38‑1 (MQ=255)

                                       gccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCa                                                       <  2:2635557‑M2/40‑1 (MQ=255)

                                       gccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCa                                                       <  2:4459811‑M2/40‑1 (MQ=255)

                                       gccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCa                                                       <  1:1316842‑M2/40‑1 (MQ=255)

                                         caggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAgg                                                     >  2:1605597‑M2/52‑93 (MQ=255)

                                         caggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAgg                                                     >  1:2251003‑M2/52‑93 (MQ=255)

                                          aggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                    >  2:3073151‑M2/51‑93 (MQ=255)

                                          aggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                    <  2:1915925‑M2/43‑1 (MQ=255)

                                          aggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                    <  1:3751940‑M2/43‑1 (MQ=255)

                                          agcttgtcttctaatacccggcgcccgagacccgcacgcgtccacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                    <  1:1214339‑M2/43‑1 (MQ=255)

                                           ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGTGATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   <  1:1853832‑M2/44‑1 (MQ=255)

                                            gttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAg                                                  >  2:76146‑M2/49‑93 (MQ=255)

                                            gttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAg                                                  >  1:4286813‑M2/49‑93 (MQ=255)

                                            gttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAg                                                  <  1:1579702‑M2/45‑1 (MQ=255)

                                             ttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAgg                                                 >  1:3010940‑M2/48‑93 (MQ=255)

                                             ttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAgg                                                 >  2:2867987‑M2/48‑93 (MQ=255)

                                               ggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGGATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGAt                                               >  2:2706464‑M2/46‑93 (MQ=255)

                                                gcttctaatacccggcgaccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATc                                              <  2:3714403‑M2/49‑1 (MQ=255)

                                                 cttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATcc                                             <  1:3320991‑M2/50‑1 (MQ=255)

                                                   tctaatacccggcgcccgagccccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGc                                           <  2:3391449‑M2/52‑1 (MQ=255)

                                                     taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTa                                         <  2:487144‑M2/54‑1 (MQ=255)

                                                       atacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAAt                                       <  1:4056396‑M2/56‑1 (MQ=255)

                                                        tacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                      <  2:4413978‑M2/57‑1 (MQ=255)

                                                        tacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                      <  2:2295284‑M2/57‑1 (MQ=255)

                                                        tacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                      >  1:3739115‑M2/37‑93 (MQ=255)

                                                            cggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTAt                                  <  2:3122730‑M2/61‑1 (MQ=255)

                                                            cggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTAt                                  <  1:820634‑M2/61‑1 (MQ=255)

                                                                gcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGAt                              <  1:2055430‑M2/65‑1 (MQ=255)

                                                                 cccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATa                             >  1:3547248‑M2/28‑93 (MQ=255)

                                                                  ccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaa                            >  2:918535‑M2/27‑93 (MQ=255)

                                                                   cgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaa                           <  2:1223126‑M2/68‑1 (MQ=255)

                                                                    gagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaa                          >  1:57671‑M2/25‑93 (MQ=255)

                                                                     agacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                         <  2:631114‑M2/70‑1 (MQ=255)

                                                                     agacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                         <  2:1391816‑M2/70‑1 (MQ=255)

                                                                     agacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                         <  2:2064031‑M2/70‑1 (MQ=255)

                                                                      gacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAAt                        <  2:3291783‑M2/71‑1 (MQ=255)

                                                                      gacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAAt                        <  2:759377‑M2/71‑1 (MQ=255)

                                                                      gacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGGAGGATCCGCTAATCTTATGGATAAAAAt                        >  2:829546‑M2/23‑93 (MQ=255)

                                                                          cgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTa                    <  1:393111‑M2/75‑1 (MQ=255)

                                                                           gcacgtgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTAt                   <  1:4307985‑M2/76‑1 (MQ=255)

                                                                           gcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTAt                   >  2:1692177‑M2/18‑93 (MQ=255)

                                                                            cacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATg                  <  1:867539‑M2/77‑1 (MQ=255)

                                                                             acgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATgg                 >  2:3421679‑M2/16‑93 (MQ=255)

                                                                             acgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATgg                 >  1:4506548‑M2/16‑93 (MQ=255)

                                                                              cgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGc                >  2:1556533‑M2/15‑93 (MQ=255)

                                                                              cgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGc                >  2:1093926‑M2/15‑93 (MQ=255)

                                                                              cgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGc                <  1:74935‑M2/79‑1 (MQ=255)

                                                                                  tcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAg            >  1:3438608‑M2/11‑93 (MQ=255)

                                                                                     acctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaa         >  2:1425552‑M2/8‑93 (MQ=255)

                                                                                      cctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaaa        >  2:1489669‑M2/7‑93 (MQ=255)

                                                                                        tgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgt      <  2:1808207‑M2/89‑1 (MQ=255)

                                                                                         gcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgtg     <  1:3884842‑M2/90‑1 (MQ=255)

                                                                                           aGGTATGGAGAAACAGTAGATAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGtgtg   <  2:4581477‑M2/92‑1 (MQ=255)

                                                                                           aGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGtgtg   >  1:3612813‑M2/2‑93 (MQ=255)

                                                                                            ggTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGa  <  1:4449444/93‑1 (MQ=255)

                                                                                            ggTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGa  >  1:165726/1‑93 (MQ=255)

                                                                                            ggTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGa  <  1:178471/93‑1 (MQ=255)

                                                                                            ggTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGa  <  2:2517200/93‑1 (MQ=255)

GTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGA  >  NC_000913/69981‑70165

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