BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,004,812	T→C	R448R (CGT→CGC)	recQ →	ATP‑dependent DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,004,812	0	T	C	100.0%	727.5 / NA	225	R448R (CGT→CGC)	recQ	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (120/105); total (120/105)

Only 100 of 238 total aligned reads displayed.

TCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCT  >  NC_000913/4004722‑4004901

tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCAt                                                                                         >  2:2455483/1‑93 (MQ=255)

 cAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATc                                                                                        >  2:5405339/1‑93 (MQ=255)

  aGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATcc                                                                                       <  1:2723180/93‑1 (MQ=255)

  aGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATcc                                                                                       <  2:2019761/93‑1 (MQ=255)

  aGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATcc                                                                                       <  2:4764614/93‑1 (MQ=255)

   gATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcg                                                                                      >  1:1158306/1‑93 (MQ=255)

    aTTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcgc                                                                                     <  1:4851430/93‑1 (MQ=255)

    aTTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcgc                                                                                     >  1:6152815/1‑93 (MQ=255)

       gCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGAc                                                                                  <  1:3608948/93‑1 (MQ=255)

        cccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACt                                                                                 <  2:3289222/93‑1 (MQ=255)

        cccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACt                                                                                 >  1:714992/1‑93 (MQ=255)

         ccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTa                                                                                <  1:5043372/93‑1 (MQ=255)

         ccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTa                                                                                >  1:3581186/1‑93 (MQ=255)

           tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATg                                                                              <  2:2162704/93‑1 (MQ=255)

           tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATg                                                                              <  2:3979424/93‑1 (MQ=255)

              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAGGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           <  2:563618/93‑1 (MQ=255)

              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           <  2:1158338/93‑1 (MQ=255)

              ccaccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTc                                                                           >  1:36872/1‑93 (MQ=255)

               caccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAACAGCGCATCCGCGACTATGGTCa                                                                          >  1:1037906/1‑93 (MQ=255)

                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         >  2:2996924/1‑93 (MQ=255)

                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         >  1:3373116/1‑93 (MQ=255)

                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         >  1:64123/1‑93 (MQ=255)

                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTCATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         <  2:3558823/93‑1 (MQ=255)

                 ccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATg                                                                        >  2:5867541/1‑93 (MQ=255)

                 ccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATg                                                                        >  2:4061614/1‑93 (MQ=255)

                 ccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATg                                                                        >  2:5998673/1‑93 (MQ=255)

                 ccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATg                                                                        >  1:4939648/1‑93 (MQ=255)

                  caTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGa                                                                       <  2:5777621/93‑1 (MQ=255)

                   aTTGGTCGTGTGAATCCGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      >  1:2647097/1‑93 (MQ=255)

                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      >  2:3707897/1‑93 (MQ=255)

                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      >  2:4429470/1‑93 (MQ=255)

                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      >  2:5949510/1‑93 (MQ=255)

                     tGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACaa                                                                    <  1:4160580/93‑1 (MQ=255)

                       gTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAAc                                                                  <  2:2563581/93‑1 (MQ=255)

                         cGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTg                                                                <  1:5439022/93‑1 (MQ=255)

                         cGTGTGAATCAGCGGTTTGGGATGGGGTATGTTGTGGAAGTGATTCGTGGTGCTAATAACCAGCTCATCCGCGACTATGGTCATGACAAACTg                                                                >  1:4058595/1‑93 (MQ=255)

                          gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGa                                                               >  1:2035470/1‑93 (MQ=255)

                           tgtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaa                                                              >  1:2928773/1‑93 (MQ=255)

                            gtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaaa                                                             <  2:2668682/93‑1 (MQ=255)

                            gtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaaa                                                             <  1:3086964/93‑1 (MQ=255)

                              gAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGt                                                           <  2:3072788/93‑1 (MQ=255)

                               aaTCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTc                                                          >  1:3007427/1‑93 (MQ=255)

                                  cAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATACCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTAt                                                       <  1:3865948/93‑1 (MQ=255)

                                  cAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTAt                                                       <  1:4330752/93‑1 (MQ=255)

                                   aGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATg                                                      >  1:635899/1‑93 (MQ=255)

                                    gCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATgg                                                     >  2:6526660/1‑93 (MQ=255)

                                      ggTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCa                                                   >  2:6602468/1‑93 (MQ=255)

                                       gTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCAt                                                  >  1:4105532/1‑93 (MQ=255)

                                        tttGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATg                                                 <  1:4605221/93‑1 (MQ=255)

                                         ttGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATgg                                                >  1:4393380/1‑93 (MQ=255)

                                         ttGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATgg                                                >  1:1606191/1‑93 (MQ=255)

                                           gggATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGc                                              >  1:6719273/1‑93 (MQ=255)

                                              aTGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGt                                           <  2:153033/93‑1 (MQ=255)

                                               tGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTg                                          >  1:5950079/1‑93 (MQ=255)

                                                 ggTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGAt                                        <  1:6216553/93‑1 (MQ=255)

                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       >  2:628604/1‑93 (MQ=255)

                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       <  1:4995433/93‑1 (MQ=255)

                                                    tATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaa                                     <  2:3784864/93‑1 (MQ=255)

                                                     aTGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaaa                                    <  2:126445/93‑1 (MQ=255)

                                                     aTGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaaa                                    <  1:915526/93‑1 (MQ=255)

                                                         ggtggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCa                                <  1:5578668/93‑1 (MQ=255)

                                                         ggtggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCa                                <  1:3713126/93‑1 (MQ=255)

                                                          gtggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCAt                               <  1:1361809/93‑1 (MQ=255)

                                                            ggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGa                             <  1:992062/93‑1 (MQ=255)

                                                             gAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGaa                            <  1:2499309/93‑1 (MQ=255)

                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         >  2:2221272/1‑93 (MQ=255)

                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         >  2:4900213/1‑93 (MQ=255)

                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         >  2:2855004/1‑93 (MQ=255)

                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       >  2:574212/1‑93 (MQ=255)

                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       <  2:3328464/93‑1 (MQ=255)

                                                                   aTTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTgg                      <  1:4115749/93‑1 (MQ=255)

                                                                   aTTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTgg                      >  1:2330605/1‑93 (MQ=255)

                                                                   aTTCGTGGTGCTAATAACCAGCGCACCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATCGGCCGTGCTAAAAGCCATGAACATTgg                      <  1:1076486/93‑1 (MQ=255)

                                                                      cGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTg                   <  2:2585744/93‑1 (MQ=255)

                                                                      cGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTg                   <  2:4832098/93‑1 (MQ=255)

                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTgg                  >  2:5504579/1‑92 (MQ=255)

                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa                  <  2:1614035/93‑1 (MQ=255)

                                                                        tggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAg                 <  2:741358/93‑1 (MQ=255)

                                                                         ggtgCTAATAACCAGCGCAACCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGc                >  1:6539485/1‑93 (MQ=255)

                                                                          gtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCg               <  2:448028/93‑1 (MQ=255)

                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCgg              >  1:1808550/1‑92 (MQ=255)

                                                                            gCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTg             <  1:2504939/93‑1 (MQ=255)

                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            <  2:1247695/93‑1 (MQ=255)

                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            >  2:3858905/1‑93 (MQ=255)

                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            >  1:2320370/1‑93 (MQ=255)

                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            >  1:3339762/1‑93 (MQ=255)

                                                                               aataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATc          >  2:4179104/1‑93 (MQ=255)

                                                                               aataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATc          >  2:3466061/1‑93 (MQ=255)

                                                                                ataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATcc         <  1:325585/93‑1 (MQ=255)

                                                                                 taaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCg        >  2:847248/1‑93 (MQ=255)

                                                                                  aaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGc       >  1:249418/1‑93 (MQ=255)

                                                                                  aaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCAATGAACATTGGGTGAGCGTGATCCGc       >  2:5095469/1‑93 (MQ=255)

                                                                                   aCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGcc      <  2:5515733/93‑1 (MQ=255)

                                                                                    ccAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCa     >  1:1739026/1‑93 (MQ=255)

                                                                                      aGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGc   >  1:2176114/1‑93 (MQ=255)

                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  >  1:3515545/1‑93 (MQ=255)

                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  >  1:3472880/1‑93 (MQ=255)

                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  >  2:1116114/1‑93 (MQ=255)

                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  >  2:4600393/1‑93 (MQ=255)

                                                                                       gcgcATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCt  >  2:3416742/1‑93 (MQ=255)

TCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCT  >  NC_000913/4004722‑4004901

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