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	289 [0]	[0] 354	[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)	229 (1.110) +CCCGAGACCCGCACGCGTCGACCTGCA	96/130	0.2	100%	coding (19/696 nt)	araD	L‑ribulose‑5‑phosphate 4‑epimerase
*	?	NC_000913	70073 =	0 (0.000)	229 (1.110) +CCCGAGACCCGCACGCGTCGACCTGCA	96/130	0.2	100%	intergenic (‑25/‑314)	araB/araC	L‑ribulokinase/transcriptional regulator for ara operon

Only 100 of 354 total aligned reads displayed.

GTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGA  >  NC_000913/69981‑70165

ttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcag                                                                                              >  1:827566‑M2/93‑93 (MQ=255)

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

 tgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcagg                                                                                             >  1:708024‑M2/92‑93 (MQ=255)

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

   ccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTa                                                                                           >  1:2122394‑M2/90‑93 (MQ=255)

   ccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTa                                                                                           >  1:4057260‑M2/90‑93 (MQ=255)

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

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

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

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

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

       atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGga                                                                                       >  2:6102‑M2/86‑93 (MQ=255)

       atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGga                                                                                       >  2:4914134‑M2/86‑93 (MQ=255)

        tgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGgag                                                                                      >  1:6532532‑M2/85‑93 (MQ=255)

         gtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccgtcgcccgagacccgcacgcgtcgacctgcaGGTATGgaga                                                                                     <  1:6226083‑M2/10‑1 (MQ=255)

           gagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGaaa                                                                                   <  2:3719036‑M2/12‑1 (MQ=255)

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

              cgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                <  2:1899451‑M2/15‑1 (MQ=255)

               gtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGt                                                                               >  1:7212230‑M2/78‑93 (MQ=255)

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

                  accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            >  2:4790940‑M2/75‑93 (MQ=255)

                  accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            >  1:187457‑M2/75‑93 (MQ=255)

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

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

                      ggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGt                                                                        <  2:1487668‑M2/23‑1 (MQ=255)

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

                       gttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGtt                                                                       >  2:7416239‑M2/70‑93 (MQ=255)

                       gttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGtt                                                                       >  1:3915664‑M2/70‑93 (MQ=255)

                        ttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTg                                                                      >  2:2196361‑M2/69‑93 (MQ=255)

                         tgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGc                                                                     <  2:1679419‑M2/26‑1 (MQ=255)

                          gtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCg                                                                    >  1:460110‑M2/67‑93 (MQ=255)

                          gtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCg                                                                    >  1:7463005‑M2/67‑93 (MQ=255)

                           tgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGa                                                                   >  2:3115655‑M2/66‑93 (MQ=255)

                           tgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGa                                                                   >  2:1873441‑M2/66‑93 (MQ=255)

                            gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGAt                                                                  >  1:4876148‑M2/65‑93 (MQ=255)

                             ttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATa                                                                 >  2:2586784‑M2/64‑93 (MQ=255)

                               ttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaa                                                               <  2:4469446‑M2/32‑1 (MQ=255)

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

                                 ggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaaa                                                             <  2:4083813‑M2/34‑1 (MQ=255)

                                 ggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaaa                                                             <  2:2807715‑M2/34‑1 (MQ=255)

                                   cagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGc                                                           <  1:5912380‑M2/36‑1 (MQ=255)

                                   cagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGc                                                           <  1:3154364‑M2/36‑1 (MQ=255)

                                    agcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                          >  1:5381537‑M2/57‑93 (MQ=255)

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

                                         caggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCCgg                                                     >  1:4917972‑M2/52‑93 (MQ=255)

                                         caggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAgg                                                     <  1:573549‑M2/42‑1 (MQ=255)

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

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

                                          aggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                    >  1:4538593‑M2/51‑93 (MQ=255)

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

                                             ttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAgg                                                 <  2:432746‑M2/46‑1 (MQ=255)

                                              tggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGa                                                >  2:5739785‑M2/47‑93 (MQ=255)

                                                gcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAggann                                              >  1:2479630‑M2/45‑91 (MQ=255)

                                                gcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATc                                              <  1:2511316‑M2/49‑1 (MQ=255)

                                                 cttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATcc                                             >  1:2668226‑M2/44‑93 (MQ=255)

                                                   tctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGc                                           <  1:4419935‑M2/52‑1 (MQ=255)

                                                     taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTa                                         >  1:3682921‑M2/40‑93 (MQ=255)

                                                     taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTa                                         >  2:6158447‑M2/40‑93 (MQ=255)

                                                      aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                        >  2:2439998‑M2/39‑93 (MQ=255)

                                                      aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                        <  2:3320643‑M2/55‑1 (MQ=255)

                                                      aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                        <  2:4964241‑M2/55‑1 (MQ=255)

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

                                                           ccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTa                                   >  2:170970‑M2/34‑93 (MQ=255)

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

                                                             ggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATg                                 <  1:2798409‑M2/62‑1 (MQ=255)

                                                             ggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATg                                 <  2:2654041‑M2/62‑1 (MQ=255)

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

                                                                  ccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGGGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaa                            >  2:6012940‑M2/27‑93 (MQ=255)

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

                                                                   cgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaa                           >  1:425078‑M2/26‑93 (MQ=255)

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

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

                                                                         ccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCt                     >  2:73615‑M2/20‑93 (MQ=255)

                                                                         ccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCt                     <  1:6757025‑M2/74‑1 (MQ=255)

                                                                          cgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTa                    >  1:4710244‑M2/19‑93 (MQ=255)

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

                                                                           gcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTAt                   >  1:6038965‑M2/18‑93 (MQ=255)

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

                                                                             acgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATgg                 <  2:6832033‑M2/78‑1 (MQ=255)

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

                                                                             acgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATgg                 <  2:4829245‑M2/78‑1 (MQ=255)

                                                                               gcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCa               <  2:939060‑M2/80‑1 (MQ=255)

                                                                                cgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCAt              <  1:3275704‑M2/81‑1 (MQ=255)

                                                                                cgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCAt              >  2:6319361‑M2/13‑93 (MQ=255)

                                                                                 gtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATa             <  2:5352615‑M2/82‑1 (MQ=255)

                                                                                 gtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATa             >  2:7018375‑M2/12‑93 (MQ=255)

                                                                                 gtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATa             <  2:6212638‑M2/82‑1 (MQ=255)

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

                                                                                  tcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAg            <  1:5299994‑M2/83‑1 (MQ=255)

                                                                                    gacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCa          >  2:2860075‑M2/9‑93 (MQ=255)

                                                                                    gacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCa          >  1:3432229‑M2/9‑93 (MQ=255)

                                                                                    gacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCa          >  1:4147752‑M2/9‑93 (MQ=255)

                                                                                     acctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaa         >  1:3020722‑M2/8‑93 (MQ=255)

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

                                                                                       ctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAg       >  2:7162829‑M2/6‑93 (MQ=255)

                                                                                        tgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgt      >  2:679601‑M2/5‑93 (MQ=255)

                                                                                          caGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgtgt    >  1:736073‑M2/3‑93 (MQ=255)

                                                                                          caGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgtgt    <  1:1755842‑M2/91‑1 (MQ=255)

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

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

GTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGA  >  NC_000913/69981‑70165

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