BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,072	T→C	D73G (GAC→GGC)	b3460 ←	LivJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,072	0	T	C	100.0%	1152.1 / NA	349	D73G (GAC→GGC)	b3460	LivJ
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (179/170); total (179/170)

Only 100 of 351 total aligned reads displayed.

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

gaagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCa                                                                                         <  2:3091975/93‑1 (MQ=255)

 aagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCat                                                                                        <  2:948291/93‑1 (MQ=255)

 aagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCat                                                                                        <  2:2666162/93‑1 (MQ=255)

 aagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCat                                                                                        <  2:381525/93‑1 (MQ=255)

  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGGATCGCCata                                                                                       >  1:6223307/1‑93 (MQ=255)

  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       >  2:5237805/1‑93 (MQ=255)

  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       >  2:6415579/1‑93 (MQ=255)

  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       >  1:5440491/1‑93 (MQ=255)

  agaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCCCCTGTTTCGGGTCACAGGCATCGCCata                                                                                       <  2:2034173/93‑1 (MQ=255)

    aaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                     <  2:2360883/93‑1 (MQ=255)

     aCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttt                                                                                    >  1:322428/1‑93 (MQ=255)

       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  2:4645527/1‑93 (MQ=255)

       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  <  2:2597946/93‑1 (MQ=255)

        gagGTGACCATTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 <  1:1761069/93‑1 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:498081/1‑93 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:3916894/1‑93 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  1:5000757/1‑93 (MQ=255)

          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               <  2:6135028/93‑1 (MQ=255)

          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               >  1:6033802/1‑93 (MQ=255)

            tGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAtt                                                                             <  1:1072607/93‑1 (MQ=255)

              aCCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTg                                                                           >  2:1914070/1‑93 (MQ=255)

               ccAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGc                                                                          >  1:1905176/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  1:3811673/1‑93 (MQ=255)

                   tCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtt                                                                      <  2:2846576/93‑1 (MQ=255)

                   tCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtt                                                                      >  1:4505433/1‑93 (MQ=255)

                    cacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGttt                                                                     <  1:748138/93‑1 (MQ=255)

                    cacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGttt                                                                     <  2:5916201/93‑1 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtttt                                                                    >  1:3073896/1‑92 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    >  2:2779491/1‑93 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    <  2:627946/93‑1 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    >  2:4853589/1‑93 (MQ=255)

                      caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   <  1:4136289/93‑1 (MQ=255)

                      caTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgt                                                                   <  1:161155/93‑1 (MQ=255)

                        tgtttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCGATATTTTACGATTTGCAGTttgttg                                                                 >  2:6481565/3‑93 (MQ=255)

                         atTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGc                                                                >  1:4420579/1‑93 (MQ=255)

                           ttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCt                                                              <  2:1748523/93‑1 (MQ=255)

                            tAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCtt                                                             <  2:1869206/93‑1 (MQ=255)

                             aaTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCttt                                                            <  1:6277422/93‑1 (MQ=255)

                              aTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTa                                                           >  1:5979861/1‑93 (MQ=255)

                              aTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGACTTGCAGTTTGGTGCCTTTa                                                           <  1:6569290/93‑1 (MQ=255)

                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:2744922/93‑1 (MQ=255)

                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          >  1:3596050/1‑93 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         <  1:5851035/93‑1 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         <  2:1419278/93‑1 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         >  1:5265916/1‑93 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         <  2:4114601/93‑1 (MQ=255)

                                 ccGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATg                                                        <  1:5230805/93‑1 (MQ=255)

                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       <  2:667728/93‑1 (MQ=255)

                                   gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      >  1:4204195/1‑93 (MQ=255)

                                   gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      <  1:5405578/93‑1 (MQ=255)

                                      gtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgcc                                                   >  1:2398007/1‑93 (MQ=255)

                                       ttAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCt                                                  <  1:2337643/93‑1 (MQ=255)

                                        tAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                 <  2:6560477/93‑1 (MQ=255)

                                         aaCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCttt                                                <  2:4564337/93‑1 (MQ=255)

                                          aCGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTa                                               >  2:1500571/1‑93 (MQ=255)

                                            gACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGc                                             <  2:5746892/93‑1 (MQ=255)

                                              cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           >  2:5071374/1‑93 (MQ=255)

                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  1:1205117/1‑93 (MQ=255)

                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  1:2258374/1‑93 (MQ=255)

                                                 tGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGa                                        >  1:3448405/1‑93 (MQ=255)

                                                 tGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGa                                        <  2:1328768/93‑1 (MQ=255)

                                                 tGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGa                                        <  2:4304769/93‑1 (MQ=255)

                                                  gTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGat                                       <  2:5326408/93‑1 (MQ=255)

                                                   ttCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGata                                      >  1:4880738/1‑93 (MQ=255)

                                                     cGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATc                                    <  1:3003742/93‑1 (MQ=255)

                                                     cGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATc                                    >  2:3375173/1‑93 (MQ=255)

                                                     cGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATc                                    <  1:2199599/93‑1 (MQ=255)

                                                      gCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   >  2:689700/1‑93 (MQ=255)

                                                      gCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATATGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   >  2:5941992/1‑93 (MQ=255)

                                                      cccACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   <  1:2577069/92‑1 (MQ=255)

                                                        cACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGc                                 <  2:121244/93‑1 (MQ=255)

                                                         aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  1:946676/1‑93 (MQ=255)

                                                         aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  2:5573100/1‑93 (MQ=255)

                                                           cGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAc                              <  2:2110982/93‑1 (MQ=255)

                                                            gCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAcc                             >  1:4702106/1‑93 (MQ=255)

                                                             cAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCg                            >  2:1783625/1‑93 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          <  1:3164765/93‑1 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  2:4396231/1‑93 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  1:2080849/1‑93 (MQ=255)

                                                                 cGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTg                        <  2:2828213/93‑1 (MQ=255)

                                                                 cGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTg                        <  1:4359047/93‑1 (MQ=255)

                                                                 cGCCTGTTTCGGGTCACAGGCATCGCCATATGTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTg                        <  2:2489468/93‑1 (MQ=255)

                                                                  cccTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGc                       <  2:3976671/92‑1 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     >  1:2823678/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     <  2:443511/93‑1 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTTATATCCGCAACCGCCTGctct                    <  2:3161571/93‑1 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    <  2:5286116/93‑1 (MQ=255)

                                                                       tttCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgc                  <  1:398037/93‑1 (MQ=255)

                                                                         tCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgcgc                >  2:2188112/1‑93 (MQ=255)

                                                                         tCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgcgc                >  2:2146566/1‑93 (MQ=255)

                                                                          cGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGcc               <  2:1963162/93‑1 (MQ=255)

                                                                             gTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGt            <  2:3458902/93‑1 (MQ=255)

                                                                               cacaGGAATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaa          >  1:1404692/1‑93 (MQ=255)

                                                                                acaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaaa         >  2:5610995/1‑93 (MQ=255)

                                                                                acaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaaa         >  2:4923361/1‑93 (MQ=255)

                                                                                acaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaaa         >  2:3007287/1‑93 (MQ=255)

                                                                                acaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaaa         >  1:661148/1‑93 (MQ=255)

                                                                                    gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  2:803118/93‑1 (MQ=255)

                                                                                    gCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTcc     <  1:1002883/93‑1 (MQ=255)

                                                                                     cATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCt    <  1:5717582/93‑1 (MQ=255)

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

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