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%	510.7 / NA	156	L504R (CTC→CGC)	recQ	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (90/66); total (90/66)

Only 100 of 176 total aligned reads displayed.

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCT‑CGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060

cGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCt                                                                                        >  2:3188606/1‑93 (MQ=255)

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

cGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCt                                                                                        >  2:447547/1‑93 (MQ=255)

  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                      <  1:1004241/93‑3 (MQ=255)

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

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

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    >  1:2826268/1‑89 (MQ=255)

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    >  1:3491097/1‑89 (MQ=255)

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    >  2:2331345/1‑89 (MQ=255)

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    >  2:3430404/1‑89 (MQ=255)

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    <  1:2035040/93‑5 (MQ=255)

    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    <  1:174424/93‑5 (MQ=255)

      ccGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑g                                                                                >  2:3958047/1‑93 (MQ=255)

      ccGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑g                                                                                <  2:2799924/93‑1 (MQ=255)

      ccGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑g                                                                                <  2:2821726/93‑1 (MQ=255)

       cGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gc                                                                               >  2:4622535/1‑92 (MQ=255)

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

          cAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcgaa                                                                            <  1:4447923/93‑5 (MQ=255)

          cAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcgaa                                                                            >  1:1834801/1‑89 (MQ=255)

                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCctct                                                                      <  1:1575982/93‑1 (MQ=255)

                aTTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCctct                                                                      <  1:897559/93‑1 (MQ=255)

                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     >  1:4643914/1‑93 (MQ=255)

                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     >  1:4643913/1‑93 (MQ=255)

                 ttCACCGCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     <  1:2756953/93‑1 (MQ=255)

                  tCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCttt                                                                    <  1:4169062/93‑1 (MQ=255)

                    aCCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGc                                                                  >  1:3617719/1‑93 (MQ=255)

                       gcGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAAc                                                               <  1:3205432/92‑1 (MQ=255)

                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            >  2:793369/1‑93 (MQ=255)

                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            >  2:872523/1‑93 (MQ=255)

                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  1:4694887/1‑93 (MQ=255)

                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        <  1:3536607/93‑1 (MQ=255)

                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        >  2:1656146/1‑93 (MQ=255)

                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgt                                                        <  1:434253/93‑1 (MQ=255)

                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGCGCTGCGCGGAG‑‑AATCCTCTTTGCAACTtgacgt                                                        >  1:3118719/1‑93 (MQ=255)

                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     <  2:1632100/93‑1 (MQ=255)

                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     >  1:947442/1‑93 (MQ=255)

                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     <  1:2701011/93‑1 (MQ=255)

                                     cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 <  2:4558063/93‑1 (MQ=255)

                                     cAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGt                                                 <  2:4558065/93‑1 (MQ=255)

                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  2:2232532/1‑93 (MQ=255)

                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  1:227477/1‑93 (MQ=255)

                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  2:1429339/1‑93 (MQ=255)

                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  1:406158/1‑93 (MQ=255)

                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               >  1:3815913/1‑93 (MQ=255)

                                        aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  1:1964142/93‑1 (MQ=255)

                                        aaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATACTCTTTGCAACTTGCCGTGCCGCGTAtg                                              >  2:3655552/1‑92 (MQ=255)

                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            <  1:766746/93‑1 (MQ=255)

                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            <  1:3045094/93‑1 (MQ=255)

                                           atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           <  1:3780355/93‑1 (MQ=255)

                                           atTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTg                                           >  1:3393073/1‑93 (MQ=255)

                                              gCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcg                                        >  2:674348/1‑93 (MQ=255)

                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgggc                                       >  2:983136/1‑93 (MQ=255)

                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       <  1:1849687/93‑1 (MQ=255)

                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  2:4260555/1‑93 (MQ=255)

                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  2:2296554/1‑93 (MQ=255)

                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:2769977/1‑93 (MQ=255)

                                                 cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  2:585380/1‑93 (MQ=255)

                                                 cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     <  1:2170521/93‑1 (MQ=255)

                                                  aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  1:2111139/1‑93 (MQ=255)

                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  2:2902141/1‑93 (MQ=255)

                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 <  1:1130636/93‑1 (MQ=255)

                                                      ttCTGCCCTACAACTGACAGCGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCTGTTCCGCCTGTCGTGGCGCTCAAAcc                                >  2:3426619/1‑93 (MQ=255)

                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                >  1:1965153/1‑93 (MQ=255)

                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                >  1:364451/1‑93 (MQ=255)

                                                       tCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCg                               >  1:2534286/1‑93 (MQ=255)

                                                         tGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaa                             >  2:2419955/1‑93 (MQ=255)

                                                          gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            <  2:4736812/93‑1 (MQ=255)

                                                          gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            >  2:1699019/1‑93 (MQ=255)

                                                          gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            <  2:1263154/93‑1 (MQ=255)

                                                           cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:243534/93‑1 (MQ=255)

                                                             cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCAGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  1:2642990/1‑93 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:2070370/1‑93 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:4012016/1‑93 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCCTGGCGCTCAAACCGAAAGCGa                        >  2:2667310/1‑93 (MQ=255)

                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCG‑‑GCCAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  2:2370215/1‑93 (MQ=255)

                                                              tACAACTGACAGAGGCCGAGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:299403/1‑93 (MQ=255)

                                                               aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       >  2:2895819/1‑93 (MQ=255)

                                                               aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       <  1:2482346/93‑1 (MQ=255)

                                                                cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATg                      >  2:1803249/1‑93 (MQ=255)

                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     >  1:503906/1‑93 (MQ=255)

                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGc                     <  1:3280406/93‑1 (MQ=255)

                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTTCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  1:3908958/93‑1 (MQ=255)

                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:3675695/1‑93 (MQ=255)

                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:920671/1‑93 (MQ=255)

                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              <  1:460841/93‑1 (MQ=255)

                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              >  1:3947004/1‑93 (MQ=255)

                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              <  1:4149643/93‑1 (MQ=255)

                                                                         gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             <  2:1352407/93‑1 (MQ=255)

                                                                         gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  1:2324589/1‑93 (MQ=255)

                                                                           ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:1574440/93‑1 (MQ=255)

                                                                            gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          >  2:364317/1‑93 (MQ=255)

                                                                             ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  1:3922352/93‑1 (MQ=255)

                                                                             ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  1:81906/93‑1 (MQ=255)

                                                                             ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  1:1311917/93‑1 (MQ=255)

                                                                              cgcgcgCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcgg        >  1:1008248/1‑93 (MQ=255)

                                                                               gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  2:1968897/93‑1 (MQ=255)

                                                                               gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       <  1:4710589/93‑1 (MQ=255)

                                                                                  cgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCgccggc    >  1:1680234/1‑93 (MQ=255)

                                                                                   gcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCa   >  2:4086669/1‑93 (MQ=255)

                                                                                    cccGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  2:1180888/1‑93 (MQ=255)

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCT‑CGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  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` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`