BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	708,076	(G)_5→6	coding (520/1407 nt)	ybfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	708,071	1	.	G	100.0%	287.0 / NA	91	Q172R (CAG→CGG)	ybfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (34/57); total (34/57)

Only 100 of 113 total aligned reads displayed.

ATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑GGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGC  >  NC_000913/707979‑708157

aTACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCa                                                                                         <  1:664819/93‑1 (MQ=255)

aTACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCa                                                                                         >  1:3852083/1‑93 (MQ=255)

aTACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAAAGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCa                                                                                         <  1:2515130/93‑1 (MQ=255)

 tACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑g                                                                                       >  2:442398/1‑93 (MQ=255)

 tACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑g                                                                                       >  1:1360557/1‑93 (MQ=255)

 tACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑g                                                                                       >  2:1352855/1‑93 (MQ=255)

 tACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑g                                                                                       <  1:3127536/93‑1 (MQ=255)

 tACGGTCCGGCTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑g                                                                                       >  2:2568981/1‑93 (MQ=255)

  aCGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑gg                                                                                      <  1:821830/93‑1 (MQ=255)

  aCGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑gg                                                                                      <  2:547543/93‑1 (MQ=255)

   cGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑ggg                                                                                     <  1:3536444/93‑1 (MQ=255)

   cGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑ggg                                                                                     <  2:1173429/93‑1 (MQ=255)

    ggTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑gggg                                                                                    <  2:3578161/93‑1 (MQ=255)

    ggTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑gggg                                                                                    <  2:1004832/93‑1 (MQ=255)

    ggTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑gggg                                                                                    <  2:825571/93‑1 (MQ=255)

     gTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑ggggg                                                                                   >  2:2625619/1‑93 (MQ=255)

        cGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑ggggggcg                                                                                >  1:947358/1‑90 (MQ=255)

           tttGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGaa                                                                              <  1:2610261/93‑1 (MQ=255)

           tttGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGaa                                                                              >  1:3147764/1‑93 (MQ=255)

            ttGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAg                                                                             <  2:824085/93‑1 (MQ=255)

             tGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGc                                                                            <  2:3766187/93‑1 (MQ=255)

             tGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGc                                                                            <  2:1852941/93‑1 (MQ=255)

             tGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGc                                                                            <  2:2082901/93‑1 (MQ=255)

               ggCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCg                                                                          >  1:4599600/1‑93 (MQ=255)

               ggCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCg                                                                          <  2:124247/93‑1 (MQ=255)

                gCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCTCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCgg                                                                         <  2:2930757/93‑1 (MQ=255)

                gCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCgg                                                                         <  2:3174890/93‑1 (MQ=255)

                gCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCgg                                                                         <  2:653069/93‑1 (MQ=255)

                gCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCgg                                                                         <  2:4445497/93‑1 (MQ=255)

                  gAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCgggg                                                                       >  1:4049416/1‑93 (MQ=255)

                  gAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCgggg                                                                       <  2:3845517/93‑1 (MQ=255)

                        aGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAAtt                                                                 <  2:3395405/93‑1 (MQ=255)

                          gTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAAtttt                                                               >  2:833863/1‑93 (MQ=255)

                          gTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAAtttt                                                               <  1:2725500/93‑1 (MQ=255)

                           ttACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTg                                                              <  1:5243/93‑1 (MQ=255)

                             aCATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGAt                                                            <  1:1143968/93‑1 (MQ=255)

                             aCATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGAt                                                            >  2:849633/1‑93 (MQ=255)

                              cATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGAtt                                                           >  1:22413/1‑93 (MQ=255)

                               aTTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTa                                                          <  2:235316/93‑1 (MQ=255)

                               aTTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTa                                                          <  1:1670199/93‑1 (MQ=255)

                                   aGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGc                                                      >  1:3594042/1‑93 (MQ=255)

                                   aGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGc                                                      <  2:4300643/93‑1 (MQ=255)

                                   aGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGc                                                      <  2:3978490/93‑1 (MQ=255)

                                    gCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCg                                                     >  2:2016060/1‑93 (MQ=255)

                                     ccAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGa                                                    <  1:4496960/93‑1 (MQ=255)

                                     ccAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGa                                                    <  1:4576903/93‑1 (MQ=255)

                                      cAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGAt                                                   >  2:2977580/1‑93 (MQ=255)

                                      cAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGAt                                                   <  2:3513443/93‑1 (MQ=255)

                                      cAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGAt                                                   >  1:4302967/1‑93 (MQ=255)

                                       aaCTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATg                                                  <  1:4578265/93‑1 (MQ=255)

                                       aaCTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATg                                                  >  1:4250757/1‑93 (MQ=255)

                                        aCTGGTCAAACTCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGc                                                 <  2:21263/93‑1 (MQ=255)

                                        aCTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGc                                                 >  1:21770/1‑93 (MQ=255)

                                          tGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTg                                               >  2:3465384/1‑93 (MQ=255)

                                           ggTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTgg                                              <  1:2882160/93‑1 (MQ=255)

                                            gTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGt                                             <  1:1646948/93‑1 (MQ=255)

                                            gTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGt                                             >  2:1852509/1‑93 (MQ=255)

                                              cAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGGATTTTGATTACGGCGATGCTGGTGc                                           >  2:3537622/1‑93 (MQ=255)

                                                aaCGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGt                                         >  1:3616595/1‑93 (MQ=255)

                                                 aCGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGtt                                        >  1:2436078/1‑93 (MQ=255)

                                                 aCGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGtt                                        <  2:426219/93‑1 (MQ=255)

                                                   gCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGa                                      <  1:2714210/93‑1 (MQ=255)

                                                     tGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGt                                    >  1:1266913/1‑93 (MQ=255)

                                                     tGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGt                                    >  1:4477933/1‑93 (MQ=255)

                                                      gTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGtt                                   <  1:2295408/93‑1 (MQ=255)

                                                      gTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGtt                                   >  1:1497422/1‑93 (MQ=255)

                                                      gTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGtt                                   <  2:4033661/93‑1 (MQ=255)

                                                       ttAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGttt                                  <  1:2356227/93‑1 (MQ=255)

                                                        tAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTtc                                 <  1:3580108/93‑1 (MQ=255)

                                                        tAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTtc                                 >  2:2744677/1‑93 (MQ=255)

                                                        tAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTtc                                 <  2:3540469/93‑1 (MQ=255)

                                                             ccgcACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTa                            <  1:4144487/93‑1 (MQ=255)

                                                             ccgcACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTa                            <  1:1950159/93‑1 (MQ=255)

                                                               gcACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACa                          >  1:583437/1‑93 (MQ=255)

                                                               gcACTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACa                          <  1:3081728/93‑1 (MQ=255)

                                                                 aCTGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACAtg                        <  1:3213680/93‑1 (MQ=255)

                                                                   tGGAGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACAtgtg                      >  2:885523/1‑93 (MQ=255)

                                                                      aGCTTTATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGAc                   >  2:4569065/1‑93 (MQ=255)

                                                                         tttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCaa                <  1:1265816/93‑1 (MQ=255)

                                                                         tttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCaa                <  2:3358040/93‑1 (MQ=255)

                                                                         tttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCaa                >  2:2310979/1‑93 (MQ=255)

                                                                         tttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCaa                >  2:815369/1‑93 (MQ=255)

                                                                          ttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAAc               <  2:2803529/93‑1 (MQ=255)

                                                                          ttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAAc               <  2:2531922/93‑1 (MQ=255)

                                                                          ttATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAAc               <  1:3867025/93‑1 (MQ=255)

                                                                           tATGCCAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACg              >  1:3795937/1‑93 (MQ=255)

                                                                               ccAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATa          <  2:1183874/93‑1 (MQ=255)

                                                                                cAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATAc         >  2:3469128/1‑93 (MQ=255)

                                                                                cAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATAc         <  1:3505328/93‑1 (MQ=255)

                                                                                cAGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATAc         <  2:4182518/93‑1 (MQ=255)

                                                                                 aGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACa        >  1:3641889/1‑93 (MQ=255)

                                                                                 aGGTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACa        >  2:3529561/1‑93 (MQ=255)

                                                                                  ggTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACaa       <  1:4094990/93‑1 (MQ=255)

                                                                                   gTACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACaaa      <  2:4033156/93‑1 (MQ=255)

                                                                                    tACTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAg     <  1:1308393/93‑1 (MQ=255)

                                                                                     aCTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAgc    <  1:2601469/93‑1 (MQ=255)

                                                                                      cTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAgcg   <  2:422696/93‑1 (MQ=255)

                                                                                      cTTATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAgcg   <  1:2100521/93‑1 (MQ=255)

                                                                                       ttATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAgcgc  >  2:1934423/1‑93 (MQ=255)

                                                                                       ttATCAGGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAgcgc  >  2:4628520/1‑93 (MQ=255)

ATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCA‑GGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGC  >  NC_000913/707979‑708157

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