BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,644,603	(G)_6→5	coding (743/2313 nt)	pbpC ←	putative peptidoglycan enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,644,598	0	G	.	100.0%	728.4 / NA	160	coding (748/2313 nt)	pbpC	putative peptidoglycan enzyme
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (81/79); total (81/79)

Only 100 of 185 total aligned reads displayed.

AGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTA  >  NC_000913/2644507‑2644686

aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  1:1009371/93‑1 (MQ=255)

aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  2:2495370/93‑1 (MQ=255)

aGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgg                                                                                         <  1:3986115/93‑1 (MQ=255)

 gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                        >  1:1017338/1‑93 (MQ=255)

 gACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggg                                                                                        >  1:445269/1‑93 (MQ=255)

  aCCGGCATCCAGCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggg                                                                                       <  1:898932/93‑1 (MQ=255)

   ccGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggggg                                                                                      >  2:22485/1‑93 (MQ=255)

   ccGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggggg                                                                                      <  1:159533/93‑1 (MQ=255)

   ccGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCggggg                                                                                      <  1:2701320/93‑1 (MQ=255)

     ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    >  2:928871/1‑91 (MQ=255)

     ggCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcc                                                                                    <  1:22275/93‑3 (MQ=255)

      gCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCgggggcca                                                                                   <  2:4172375/93‑4 (MQ=255)

       cATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGAGCCAGTTGCGGCATTTGTCgggggccag                                                                                  >  2:173690/1‑89 (MQ=255)

        aTCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGgccagc                                                                                >  1:39024/1‑93 (MQ=255)

        aTCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGgccagc                                                                                >  2:374900/1‑93 (MQ=255)

         tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGgccagcc                                                                               >  2:1082475/1‑93 (MQ=255)

         tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTATC‑GGGGgccagcc                                                                               <  2:4502561/93‑1 (MQ=255)

            aaCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGa                                                                            >  2:1474695/1‑93 (MQ=255)

            aaCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGa                                                                            <  2:1639308/93‑1 (MQ=255)

             aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGAt                                                                           <  1:2638524/93‑1 (MQ=255)

              cGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATg                                                                          >  1:1967609/1‑93 (MQ=255)

               gtagtaGTTATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATgg                                                                         <  1:3591577/93‑1 (MQ=255)

               gtagtaGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATgg                                                                         <  2:11600/93‑1 (MQ=255)

                tagtagTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATggg                                                                        <  1:1392677/93‑1 (MQ=255)

                 agtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGt                                                                       <  2:2957634/93‑1 (MQ=255)

                 agtagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGt                                                                       <  1:2368159/93‑1 (MQ=255)

                   tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttc                                                                     <  1:811173/93‑1 (MQ=255)

                   tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttc                                                                     >  1:3956856/1‑93 (MQ=255)

                    agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttct                                                                    <  2:4602904/93‑1 (MQ=255)

                    agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttct                                                                    <  2:1053823/93‑1 (MQ=255)

                    agtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttct                                                                    >  1:4293557/1‑93 (MQ=255)

                     gtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttctt                                                                   <  2:853640/93‑1 (MQ=255)

                      tGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttcttc                                                                  >  1:1221410/1‑93 (MQ=255)

                      tGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttcttc                                                                  <  2:591974/93‑1 (MQ=255)

                      tGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGttcttc                                                                  >  1:2052312/1‑93 (MQ=255)

                       gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTcc                                                                 >  1:458390/1‑93 (MQ=255)

                       gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTcc                                                                 >  1:2580177/1‑93 (MQ=255)

                       gATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTcc                                                                 >  2:4646365/1‑93 (MQ=255)

                         ttttGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCt                                                               >  2:747972/1‑93 (MQ=255)

                            tGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGa                                                            >  2:2497534/1‑93 (MQ=255)

                             gTCGCTTGTGCTTTTACCGAGCAGCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGAc                                                           <  1:2128241/93‑1 (MQ=255)

                              tCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGAct                                                          <  2:88328/93‑1 (MQ=255)

                               cGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGActc                                                         <  1:69328/93‑1 (MQ=255)

                               cGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGActc                                                         <  2:4480860/93‑1 (MQ=255)

                                gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGActct                                                        >  2:4095716/1‑93 (MQ=255)

                                gcttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGActat                                                        >  1:3142994/1‑91 (MQ=255)

                                 cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                       >  1:3882750/1‑93 (MQ=255)

                                 cttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtt                                                       <  2:4463149/93‑1 (MQ=255)

                                  ttttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCttt                                                      >  2:3175796/1‑93 (MQ=255)

                                   tttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCtttt                                                     <  1:508183/93‑1 (MQ=255)

                                   tttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGGTCTTCCCTTGACTCtttt                                                     >  1:828533/1‑93 (MQ=255)

                                    ttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTa                                                    >  2:4116553/1‑93 (MQ=255)

                                    ttgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTa                                                    <  1:685865/93‑1 (MQ=255)

                                     tgcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc                                                   <  2:1371849/93‑1 (MQ=255)

                                     tccttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc                                                   <  2:1992032/91‑1 (MQ=255)

                                      gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                  <  2:1274540/93‑1 (MQ=255)

                                      gcttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAcc                                                  <  2:2788697/93‑1 (MQ=255)

                                        ttttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTg                                                >  1:4495579/1‑93 (MQ=255)

                                         tttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGc                                               >  1:2555970/1‑93 (MQ=255)

                                         tttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGc                                               <  1:1518379/93‑1 (MQ=255)

                                         tttACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGc                                               >  2:121135/1‑93 (MQ=255)

                                            aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            <  1:1638339/93‑1 (MQ=255)

                                            aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            <  2:3528378/93‑1 (MQ=255)

                                            aCCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCa                                            <  2:3253065/93‑1 (MQ=255)

                                                   atcatGCGCGAAAACACCGGTGCCCGTTGCGGCACTTGTC‑GGGGGCCAGCAAGATGGGTTCTTACCTTGACTCTTTTATCTGCTCACGGGccc                                     >  2:2018183/1‑93 (MQ=255)

                                                     catGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcac                                   <  2:3143628/93‑1 (MQ=255)

                                                       tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                 <  2:2013101/93‑1 (MQ=255)

                                                       tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                 >  2:501065/1‑93 (MQ=255)

                                                       tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacac                                 >  1:2773139/1‑93 (MQ=255)

                                                        gcgcGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacaca                                <  2:2583905/93‑1 (MQ=255)

                                                        gcgcGAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcacaca                                <  1:2979359/93‑1 (MQ=255)

                                                          gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              >  2:433311/1‑93 (MQ=255)

                                                          gcgAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACAcc                              >  1:3712055/1‑93 (MQ=255)

                                                            gAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                            >  2:3584333/1‑93 (MQ=255)

                                                            gAAAACAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCtt                            <  2:3925780/93‑1 (MQ=255)

                                                              aaaCAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGc                          <  1:3035251/93‑1 (MQ=255)

                                                              aaaCAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGc                          <  2:1714262/93‑1 (MQ=255)

                                                               aaCAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCa                         >  1:3078387/1‑93 (MQ=255)

                                                               aaCAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCa                         <  2:342483/93‑1 (MQ=255)

                                                                aCAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCAc                        >  2:3252563/1‑93 (MQ=255)

                                                                 cAGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGGTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACg                       >  2:526054/1‑93 (MQ=255)

                                                                  aGCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACgg                      >  2:4275969/1‑93 (MQ=255)

                                                                   gCGGTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGGTATTCCCTTGACTCTTTTACCTGCTCACGGCACCACACACCTTGCACGGc                     >  2:363047/1‑93 (MQ=255)

                                                                      gTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCAt                  >  2:320727/1‑93 (MQ=255)

                                                                      gTGCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCAt                  <  1:1667660/93‑1 (MQ=255)

                                                                        gCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATcc                <  1:3374957/93‑1 (MQ=255)

                                                                        gCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATcc                <  1:215172/93‑1 (MQ=255)

                                                                        gCCAGTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATcc                >  1:4258061/1‑93 (MQ=255)

                                                                            gTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc            <  1:3493607/93‑1 (MQ=255)

                                                                            gTTGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTc            >  1:2529963/1‑93 (MQ=255)

                                                                             ttGCGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCg           >  2:3060000/1‑93 (MQ=255)

                                                                                gggCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGc        <  2:387689/92‑1 (MQ=255)

                                                                                cGGCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGc        <  1:1208305/93‑1 (MQ=255)

                                                                                 ggCATTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCa       <  1:2453676/93‑1 (MQ=255)

                                                                                  gCATTTGTC‑GGGCGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTTCCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCAc      <  2:3377836/93‑1 (MQ=255)

                                                                                    aTTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACGACACACCTTGCACGGCAATCCGTTCGAGCACtt    >  1:3258881/1‑93 (MQ=255)

                                                                                    aTTTGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGTACGGCCATCCGTTCGAGCACtt    >  2:2366037/1‑93 (MQ=255)

                                                                                     tttGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACttt   <  1:4145886/93‑1 (MQ=255)

                                                                                     tttGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACttt   <  1:749826/93‑1 (MQ=255)

                                                                                      ttGTC‑GGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTtta  <  1:172540/93‑1 (MQ=255)

AGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTA  >  NC_000913/2644507‑2644686

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`