BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,496,430	T→C	I322I (ATT→ATC)	cysG →	uroporphyrinogen III methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,496,430	0	T	C	100.0%	609.1 / NA	189	I322I (ATT→ATC)	cysG	uroporphyrinogen III methylase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (92/97); total (92/97)

Only 100 of 198 total aligned reads displayed.

AGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAG  >  NC_000913/3496340‑3496519

aGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATccc                                                                                         <  1:2005841/93‑1 (MQ=255)

aGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATccc                                                                                         >  1:2313257/1‑93 (MQ=255)

 ggCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCg                                                                                        <  2:2743951/93‑1 (MQ=255)

 ggCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCg                                                                                        <  1:945109/93‑1 (MQ=255)

  gCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGt                                                                                       >  2:2993519/1‑93 (MQ=255)

   cAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGtt                                                                                      >  1:1452400/1‑93 (MQ=255)

    aaaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtc                                                                                     >  2:2792412/1‑93 (MQ=255)

     aaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTtct                                                                                    <  1:2862103/93‑1 (MQ=255)

       cgcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCg                                                                                  >  1:1806829/1‑93 (MQ=255)

        gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCgg                                                                                 <  1:1019143/93‑1 (MQ=255)

          gtggtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGCGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCggtg                                                                               >  1:3602438/1‑93 (MQ=255)

          gtggtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATACCGTTCTCgggg                                                                               >  2:819957/1‑91 (MQ=255)

           tggtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCggtgg                                                                              >  2:2942468/1‑93 (MQ=255)

           tggtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCggtgg                                                                              <  2:1650453/93‑1 (MQ=255)

            ggtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCggtggt                                                                             <  1:651349/93‑1 (MQ=255)

             gtgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGtt                                                                            >  1:3929314/1‑93 (MQ=255)

              tgCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTc                                                                           <  1:4234447/93‑1 (MQ=255)

               gCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTcc                                                                          <  2:2079981/93‑1 (MQ=255)

               gCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTcc                                                                          <  2:1526943/93‑1 (MQ=255)

                  gCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCggg                                                                       >  2:1215915/1‑93 (MQ=255)

                  gCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCggg                                                                       >  1:3247832/1‑93 (MQ=255)

                    tGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCgggca                                                                     <  1:4539437/93‑3 (MQ=255)

                      aaaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAtt                                                                   <  1:1005462/93‑1 (MQ=255)

                      aaaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCAtt                                                                   <  2:17618/93‑1 (MQ=255)

                       aaGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTa                                                                  >  2:4286723/1‑93 (MQ=255)

                        aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                 <  1:3134636/93‑1 (MQ=255)

                        aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAc                                                                 >  1:2216373/1‑93 (MQ=255)

                         ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAcc                                                                >  1:2147401/1‑93 (MQ=255)

                         ggTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTAcc                                                                >  2:3654348/1‑93 (MQ=255)

                             gCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAg                                                            <  1:2103969/93‑1 (MQ=255)

                              cGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGc                                                           <  1:951492/93‑1 (MQ=255)

                               gATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCt                                                          >  1:4641547/1‑93 (MQ=255)

                                aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCtt                                                         >  1:815283/1‑93 (MQ=255)

                                aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCtt                                                         >  2:4370232/1‑93 (MQ=255)

                                aTCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCtt                                                         <  1:1674749/93‑1 (MQ=255)

                                 tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                        <  2:3251754/93‑1 (MQ=255)

                                 tCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                        >  1:2474196/1‑93 (MQ=255)

                                 tCCGTATATTTTTGGCCGTGTTCGCGAAGAGCTGAAAACACTGTGCAACGCGTGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTc                                                        <  1:4424613/93‑1 (MQ=255)

                                   cGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTg                                                      >  2:3388748/1‑93 (MQ=255)

                                     tttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGt                                                    >  1:1224254/1‑93 (MQ=255)

                                      ttATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGATTTACGGCAGTTTCTgttg                                                   >  2:3065940/1‑90 (MQ=255)

                                        aTTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                 <  2:4019923/93‑1 (MQ=255)

                                        aTTTTTGGCCGTGGTGGCGAAGAGCTGCAAACACTGTGCAACGCGTGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGc                                                 <  1:3160662/93‑1 (MQ=255)

                                         tttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGct                                                >  2:551900/1‑93 (MQ=255)

                                         tttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGct                                                >  2:63889/1‑93 (MQ=255)

                                          ttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGctc                                               >  1:3920525/1‑93 (MQ=255)

                                            ttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTg                                             <  1:3578922/93‑1 (MQ=255)

                                            ttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTg                                             >  2:2336597/1‑93 (MQ=255)

                                              ggCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGcc                                           >  2:3181599/1‑93 (MQ=255)

                                               gCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCt                                          >  1:2580671/1‑93 (MQ=255)

                                                ccGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTa                                         >  2:2795907/1‑93 (MQ=255)

                                                ccGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCAGGGCATTACCGCAGCTTCTGGTTGCTCTGCCtt                                         >  1:1343413/1‑92 (MQ=255)

                                                ccGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCATTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTa                                         <  2:1850018/93‑1 (MQ=255)

                                                 cGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAt                                        >  1:3218795/1‑93 (MQ=255)

                                                 cGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAt                                        <  2:4496733/93‑1 (MQ=255)

                                                  gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTAtt                                       >  1:3383694/1‑93 (MQ=255)

                                                   tggtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                      <  2:2323110/93‑1 (MQ=255)

                                                     gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCgg                                    <  1:4049922/93‑1 (MQ=255)

                                                     gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCgg                                    <  2:2867467/93‑1 (MQ=255)

                                                      tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                   >  2:196290/1‑93 (MQ=255)

                                                      tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                   >  2:2431757/1‑93 (MQ=255)

                                                      tggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                   <  2:3562983/93‑1 (MQ=255)

                                                       ggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGt                                  >  2:580826/1‑93 (MQ=255)

                                                       ggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGt                                  >  2:1906256/1‑93 (MQ=255)

                                                         cGAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAt                                <  1:1370609/93‑1 (MQ=255)

                                                          gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                               >  2:4073957/1‑93 (MQ=255)

                                                          gAAGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                               >  2:3459712/1‑93 (MQ=255)

                                                           aaGAGCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTc                              <  2:4633897/93‑1 (MQ=255)

                                                            agagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTcc                             <  1:4449796/93‑1 (MQ=255)

                                                             gagCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCa                            >  2:4421276/1‑93 (MQ=255)

                                                              agCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                           <  1:2037562/93‑1 (MQ=255)

                                                              agCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                           <  1:980979/93‑1 (MQ=255)

                                                              agCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCAc                           <  2:1214856/93‑1 (MQ=255)

                                                               gCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACt                          >  2:2490550/1‑93 (MQ=255)

                                                               gCTGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACt                          >  1:3277047/1‑93 (MQ=255)

                                                                 tGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCa                        >  2:141027/1‑93 (MQ=255)

                                                                 tGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCa                        <  2:746143/93‑1 (MQ=255)

                                                                 tGGAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCa                        <  1:292914/93‑1 (MQ=255)

                                                                  ggAAACACTGTGCAACGCGGGCATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGCTGCTCTGCCTATTAGGGTAGTCCACTCAc                       >  1:1596906/1‑93 (MQ=255)

                                                                   gAAACACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTACACTCACg                      >  2:4247074/1‑93 (MQ=255)

                                                                    aaaCACTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGc                     <  1:4607072/93‑1 (MQ=255)

                                                                       cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                  <  2:1467242/93‑1 (MQ=255)

                                                                       cacTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATc                  >  2:108409/1‑93 (MQ=255)

                                                                         cTGTGCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgc                >  1:3742507/1‑93 (MQ=255)

                                                                            tgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAt             >  1:4172031/1‑93 (MQ=255)

                                                                            tgCAACGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAt             <  2:3751976/93‑1 (MQ=255)

                                                                               aaCGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTAt          <  2:144046/93‑1 (MQ=255)

                                                                               aaCGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTAt          <  2:254401/93‑1 (MQ=255)

                                                                                aCGCGGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATg         >  2:434048/1‑93 (MQ=255)

                                                                                 cgcgGGTGTCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc        <  2:1684115/93‑1 (MQ=255)

                                                                                 cgcgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc        <  1:3321369/93‑1 (MQ=255)

                                                                                 cgcgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGc        <  1:2660211/93‑1 (MQ=255)

                                                                                   cgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGccc      <  2:2847734/93‑1 (MQ=255)

                                                                                   cgGGTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGccc      >  1:1037014/1‑93 (MQ=255)

                                                                                    gggTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCa     >  2:2592096/1‑93 (MQ=255)

                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   <  1:4486168/93‑1 (MQ=255)

                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   <  2:3950375/93‑1 (MQ=255)

                                                                                      gTATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCaga   <  1:270687/93‑1 (MQ=255)

                                                                                       tATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag  >  2:687692/1‑93 (MQ=255)

                                                                                       tATCCCGTTCTCGGTGGTTCCGGGCATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCagag  <  1:2021869/93‑1 (MQ=255)

AGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAG  >  NC_000913/3496340‑3496519

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