BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,004,979	+G	coding (1511/1833 nt)	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,979	1	.	G	100.0%	948.8 / NA	255	L504R (CTC→CGC)	recQ	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (134/121); total (134/121)

Only 100 of 273 total aligned reads displayed.

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCT‑CGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060

cGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCt                                                                                      <  1:5144650/93‑1 (MQ=255)

  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                    >  2:2866758/1‑91 (MQ=255)

   gATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcg                                                                                   <  2:1126584/93‑4 (MQ=255)

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGACCGGTGCtgcgc                                                                                  >  2:404837/1‑89 (MQ=255)

       cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                                             <  1:3868945/93‑2 (MQ=255)

       cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGCGCTGCGCg‑‑gc                                                                             >  1:6637189/1‑92 (MQ=255)

        gCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcg                                                                            >  2:4204738/1‑91 (MQ=255)

         ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                           >  1:1061408/1‑90 (MQ=255)

         ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                           >  1:1547140/1‑90 (MQ=255)

                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      >  1:425766/1‑93 (MQ=255)

                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      >  1:4409878/1‑93 (MQ=255)

                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      >  1:5035237/1‑93 (MQ=255)

                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCctct                                                                      >  1:881846/1‑93 (MQ=255)

                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCtt                                                                     <  2:5462687/93‑1 (MQ=255)

                  tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                                    <  2:6285110/93‑1 (MQ=255)

                  tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCttt                                                                    <  2:3949102/93‑1 (MQ=255)

                   cACCTCGGCCTCGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTg                                                                   <  2:1434071/93‑1 (MQ=255)

                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGc                                                                  >  2:2351807/1‑93 (MQ=255)

                      cTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCaa                                                                >  2:5528345/1‑93 (MQ=255)

                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAAc                                                               <  2:3678340/93‑1 (MQ=255)

                        cGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACt                                                              >  2:3450051/1‑93 (MQ=255)

                         ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                                             >  1:3440197/1‑93 (MQ=255)

                         ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                                             <  2:2248876/93‑1 (MQ=255)

                         ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACtt                                                             >  2:3827808/1‑93 (MQ=255)

                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtg                                                            >  2:89049/1‑93 (MQ=255)

                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTACTTGCAACttc                                                            >  2:5797874/1‑92 (MQ=255)

                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  2:2445936/1‑93 (MQ=255)

                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgcc                                                          >  1:3532876/1‑93 (MQ=255)

                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        >  1:1970437/1‑93 (MQ=255)

                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgt                                                        <  1:942699/93‑1 (MQ=255)

                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtg                                                       <  1:1137982/93‑1 (MQ=255)

                                ttaCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                      >  1:4206525/3‑93 (MQ=255)

                                tGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                      >  1:2798286/1‑93 (MQ=255)

                                tGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgc                                                      >  2:5533596/1‑93 (MQ=255)

                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     <  1:4692807/93‑1 (MQ=255)

                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     >  2:4860133/1‑93 (MQ=255)

                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgcc                                                     <  1:6175279/93‑1 (MQ=255)

                                  aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgccg                                                    <  1:4996926/93‑1 (MQ=255)

                                  aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgccg                                                    >  1:4476028/1‑93 (MQ=255)

                                  aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTtgccgtgccg                                                    >  2:6409229/1‑93 (MQ=255)

                                   cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgc                                                   >  1:3569266/1‑93 (MQ=255)

                                   cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgc                                                   >  2:4986383/1‑93 (MQ=255)

                                    gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  2:5408098/1‑93 (MQ=255)

                                    gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  >  2:3084029/1‑93 (MQ=255)

                                    gCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCcgcg                                                  <  1:1663870/93‑1 (MQ=255)

                                     cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 >  2:4553041/1‑93 (MQ=255)

                                     cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 >  1:5321941/1‑93 (MQ=255)

                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  1:3796306/1‑93 (MQ=255)

                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                <  2:2737790/93‑1 (MQ=255)

                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                <  1:3035372/93‑1 (MQ=255)

                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               <  1:1182913/93‑1 (MQ=255)

                                        aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  2:6091057/93‑1 (MQ=255)

                                         atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             >  2:4739655/1‑93 (MQ=255)

                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            <  1:3732760/93‑1 (MQ=255)

                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            >  1:5074929/1‑93 (MQ=255)

                                            ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          <  1:2013969/93‑1 (MQ=255)

                                             tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:4398854/93‑1 (MQ=255)

                                             tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  1:2285299/93‑1 (MQ=255)

                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:4065274/1‑93 (MQ=255)

                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  2:2109755/1‑93 (MQ=255)

                                                    cATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaaa                                  <  1:2202581/93‑1 (MQ=255)

                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  2:1673591/93‑1 (MQ=255)

                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  1:5803758/1‑93 (MQ=255)

                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  2:603041/93‑1 (MQ=255)

                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                >  1:6239538/1‑93 (MQ=255)

                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                <  2:1225383/93‑1 (MQ=255)

                                                        cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              <  1:542731/93‑1 (MQ=255)

                                                        cTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGa                              <  2:5313187/93‑1 (MQ=255)

                                                         tGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             <  1:5566516/93‑1 (MQ=255)

                                                          gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            >  1:5236077/1‑93 (MQ=255)

                                                          gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            <  2:669212/93‑1 (MQ=255)

                                                           cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:2855157/93‑1 (MQ=255)

                                                             cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  2:834880/1‑93 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  2:5393313/93‑1 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  1:700555/93‑1 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  2:1501340/93‑1 (MQ=255)

                                                                cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGTAACTTGTCGTGCCGCGTATCGTGTTGCTCAAACCGAAAGCGATg                      >  2:6648105/1‑93 (MQ=255)

                                                                cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      <  1:2913218/93‑1 (MQ=255)

                                                                cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      >  1:6453401/1‑93 (MQ=255)

                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  1:2659361/1‑93 (MQ=255)

                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  2:5376001/93‑1 (MQ=255)

                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  2:1390165/93‑1 (MQ=255)

                                                                    tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGa                  >  1:5755617/1‑93 (MQ=255)

                                                                     gACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaa                 >  2:3462609/1‑93 (MQ=255)

                                                                      aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  1:4759993/93‑1 (MQ=255)

                                                                      aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  1:3240661/93‑1 (MQ=255)

                                                                      aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                >  1:5284617/1‑93 (MQ=255)

                                                                       cAGAGGCCGCGCGCCCGTTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  2:2037664/93‑1 (MQ=255)

                                                                       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:3180982/93‑1 (MQ=255)

                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              >  2:1952046/1‑93 (MQ=255)

                                                                         gagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  1:4388157/1‑93 (MQ=255)

                                                                         gagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  2:3234725/1‑93 (MQ=255)

                                                                           ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           >  1:6353530/1‑93 (MQ=255)

                                                                              cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        <  2:5546961/93‑1 (MQ=255)

                                                                              cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        >  2:5153682/1‑93 (MQ=255)

                                                                               gcgcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  1:3944096/1‑93 (MQ=255)

                                                                                 gcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcgg     <  1:4536870/93‑1 (MQ=255)

                                                                                  cgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGCTCggcggc    >  2:1042163/1‑93 (MQ=255)

                                                                                    cccGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  1:245047/1‑93 (MQ=255)

                                                                                    cccGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  <  2:2097224/93‑1 (MQ=255)

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCT‑CGCG‑‑GAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060

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: `‑`