BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,474,849	20 bp→GTATTTCCCGGCATCTTTATAGCGATAGCA	intergenic (+24/+2)	yjgL → / ← argI	hypothetical protein/ornithine carbamoyltransferase 1

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4474849	4474868	20	193 [0]	[0] 215	yjgL/argI	hypothetical protein/ornithine carbamoyltransferase 1

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4474848	0 (0.000)	131 (1.080) +GTATTTCCCGGCATCTTTATAGCGATAGCA	75/124	0.1	100%	intergenic (+23/+22)	yjgL/argI	hypothetical protein/ornithine carbamoyltransferase 1
*	?	NC_000913	4474869 =	0 (0.000)	131 (1.080) +GTATTTCCCGGCATCTTTATAGCGATAGCA	75/124	0.1	100%	intergenic (+44/+1)	yjgL/argI	hypothetical protein/ornithine carbamoyltransferase 1

Only 100 of 193 total aligned reads displayed.

AATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAA  >  NC_000913/4474757‑4474940

aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGttt                                                                                              >  1:789358/1‑93 (MQ=255)

 aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtttt                                                                                             >  2:3953431/1‑92 (MQ=255)

   ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                                                           <  2:2087834‑M1/93‑4 (MQ=255)

     aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                                                         >  2:1789479‑M1/1‑88 (MQ=255)

     aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                                                         <  1:368232‑M1/93‑6 (MQ=255)

     aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                                                         <  2:3374637‑M1/93‑6 (MQ=255)

     aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                                                         >  1:2780744‑M1/1‑88 (MQ=255)

      ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtattt                                                                                        >  1:1680397‑M1/1‑87 (MQ=255)

        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      <  2:1324248‑M1/93‑9 (MQ=255)

        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      >  2:3723324‑M1/1‑85 (MQ=255)

        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      >  1:2510700‑M1/1‑85 (MQ=255)

        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      >  1:4223780‑M1/1‑85 (MQ=255)

        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      >  1:2232630‑M1/1‑85 (MQ=255)

         gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                                                     <  2:2415231‑M1/93‑10 (MQ=255)

             aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                                                 <  1:4407574‑M1/93‑14 (MQ=255)

             aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                                                 >  2:4159878‑M1/1‑80 (MQ=255)

              tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                                                >  2:4609846‑M1/1‑79 (MQ=255)

                 aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                                                             <  1:4661119‑M1/93‑18 (MQ=255)

                 aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                                                             <  2:583285‑M1/93‑18 (MQ=255)

                 aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                                                             >  2:2528334‑M1/1‑76 (MQ=255)

                  cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                                                            >  2:3073767‑M1/1‑75 (MQ=255)

                  cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                                                            >  1:2914544‑M1/1‑75 (MQ=255)

                   gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                                                           <  1:604422‑M1/93‑20 (MQ=255)

                   gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                                                           <  1:2412716‑M1/93‑20 (MQ=255)

                   gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                                                           >  2:1276962‑M1/1‑74 (MQ=255)

                   gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                                                           >  1:1257783‑M1/1‑74 (MQ=255)

                    tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                                                          >  1:2132486‑M1/1‑73 (MQ=255)

                     ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                                                         <  1:3018678‑M1/93‑22 (MQ=255)

                     ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                                                         >  1:3505764‑M1/1‑72 (MQ=255)

                       aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                                                       >  2:2605331‑M1/1‑70 (MQ=255)

                        gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                                                      <  1:2748881‑M1/93‑25 (MQ=255)

                           aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                                                   >  1:3299056‑M1/1‑66 (MQ=255)

                             aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                                                 >  1:3271599‑M1/1‑64 (MQ=255)

                              gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccctgcatctttatagcgatagca                                                                >  2:3317154‑M1/1‑63 (MQ=255)

                              gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                                                <  1:4065657‑M1/93‑31 (MQ=255)

                                acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                                              <  1:340888‑M1/93‑33 (MQ=255)

                                acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                                              >  2:2944466‑M1/1‑61 (MQ=255)

                                    tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                                                          <  1:4039905‑M1/93‑37 (MQ=255)

                                     ttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt                                                         <  1:1407119‑M1/93‑38 (MQ=255)

                                       gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                                                       <  2:1125997‑M1/93‑40 (MQ=255)

                                       gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                                                       >  2:1178017‑M1/1‑54 (MQ=255)

                                        ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                                                      <  1:378060‑M1/93‑41 (MQ=255)

                                        ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                                                      <  2:2322350‑M1/93‑41 (MQ=255)

                                         gTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                                                     <  1:3184425‑M1/93‑42 (MQ=255)

                                           ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                                                   >  2:2510399‑M1/1‑50 (MQ=255)

                                            gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                                                  >  2:1274176‑M1/1‑49 (MQ=255)

                                             aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                                                 <  2:3767591‑M1/93‑46 (MQ=255)

                                             aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                                                 >  2:1884684‑M1/1‑48 (MQ=255)

                                              agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                                                >  1:3517189‑M1/1‑47 (MQ=255)

                                               gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                                               <  1:1377050‑M1/93‑48 (MQ=255)

                                                   aGACCTATCACTGCATAAAATAATGTGATGCCGGGATTGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga                                           <  1:1621536‑M1/93‑52 (MQ=255)

                                                    gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac                                          <  1:4290575‑M1/93‑53 (MQ=255)

                                                    gACCTATCACAGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac                                          <  1:175000‑M1/93‑53 (MQ=255)

                                                      ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca                                        <  1:3627434‑M1/93‑55 (MQ=255)

                                                      ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca                                        >  1:1140704‑M1/1‑39 (MQ=255)

                                                        tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc                                      >  1:1900283‑M1/1‑37 (MQ=255)

                                                          tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac                                    <  2:4050593‑M1/93‑59 (MQ=255)

                                                          tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac                                    >  1:1175883‑M1/1‑35 (MQ=255)

                                                          tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac                                    <  1:3821052‑M1/93‑59 (MQ=255)

                                                           cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc                                   >  1:1011401‑M1/1‑34 (MQ=255)

                                                            aCTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccg                                  >  1:1487183‑M1/1‑33 (MQ=255)

                                                            aCTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccg                                  >  2:3277186‑M1/1‑33 (MQ=255)

                                                              tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct                                <  1:2395477‑M1/93‑63 (MQ=255)

                                                              tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct                                >  1:1119332‑M1/1‑31 (MQ=255)

                                                                cataaaataatgtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgcttt                              <  1:1031394‑M1/89‑65 (MQ=255)

                                                                cataaaataatgtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgcttt                              >  1:4747234‑M1/5‑29 (MQ=255)

                                                                  taaaataatgTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttga                            <  1:2318563‑M1/91‑67 (MQ=255)

                                                                  taaaataatgTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttga                            <  1:2669800‑M1/91‑67 (MQ=255)

                                                                   aaaataatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgat                           <  2:2478406‑M1/92‑68 (MQ=255)

                                                                    aaataaTGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgata                          <  2:2162209‑M1/93‑69 (MQ=255)

                                                                    aaataaTGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgata                          <  2:3220595‑M1/93‑69 (MQ=255)

                                                                     aataatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatag                         <  1:3237595‑M1/93‑70 (MQ=255)

                                                                     aataatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatag                         >  2:886105‑M1/1‑24 (MQ=255)

                                                                     aataatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatag                         <  1:1584345‑M1/93‑70 (MQ=255)

                                                                      ataatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagt                        <  2:4539194‑M1/93‑71 (MQ=255)

                                                                       taatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagta                       >  2:671875‑M1/1‑22 (MQ=255)

                                                                       taatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagta                       <  2:3614466‑M1/93‑72 (MQ=255)

                                                                        aatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtat                      >  1:3474333‑M1/1‑21 (MQ=255)

                                                                         atGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatg                     <  2:1471860‑M1/93‑74 (MQ=255)

                                                                          tgtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgc                    <  1:3129244‑M1/93‑75 (MQ=255)

                                                                           gtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgca                   <  2:2137770‑M1/93‑76 (MQ=255)

                                                                             gATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcata                 >  2:4298564‑M1/1‑16 (MQ=255)

                                                                             gATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcaca                 >  1:1231192‑M1/1‑16 (MQ=255)

                                                                              aTGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatac                <  2:1339055‑M1/93‑79 (MQ=255)

                                                                              aTGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatac                <  2:2133254‑M1/93‑79 (MQ=255)

                                                                               tGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacg               >  2:1463159‑M1/1‑14 (MQ=255)

                                                                                gCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacgg              <  2:3644368‑M1/93‑81 (MQ=255)

                                                                                  cGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggtt            <  1:471118‑M1/93‑83 (MQ=255)

                                                                                  cGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacgggt            >  1:895925‑M1/1‑11 (MQ=255)

                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccggtttgatagtatgcatacgggtt           >  2:3439581‑M1/1‑10 (MQ=255)

                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttt           >  2:1432983‑M1/1‑10 (MQ=255)

                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttt           <  1:687769‑M1/93‑84 (MQ=255)

                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttt           >  1:1713973‑M1/1‑10 (MQ=255)

                                                                                    ggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggtttt          >  2:161749‑M1/1‑9 (MQ=255)

                                                                                    ggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggtttt          >  1:4324531‑M1/1‑9 (MQ=255)

                                                                                     gATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttttc         <  2:2323806‑M1/93‑86 (MQ=255)

                                                                                      aTGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttttcc        <  1:20134‑M1/93‑87 (MQ=255)

                                                                                       tGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttttccg       >  1:173496‑M1/1‑6 (MQ=255)

                                                                                        ggTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttttccgc      >  2:853090‑M1/1‑5 (MQ=255)

                                                                                            tgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgctttgatagtatgcatacggttttccgcctga  <  2:4681936‑M1/93‑93 (MQ=255)

AATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAA  >  NC_000913/4474757‑4474940

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