BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,847,298	A→G	K146K (AAA→AAG)	sppA →	protease IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,847,298	0	A	G	100.0%	573.6 / NA	184	K146K (AAA→AAG)	sppA	protease IV
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (98/86); total (98/86)

Only 100 of 191 total aligned reads displayed.

TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387

tgtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGcc                                                                                         >  2:6976457/1‑93 (MQ=255)

 gtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCg                                                                                        >  2:3323596/1‑93 (MQ=255)

  tgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCgg                                                                                       >  2:4553316/1‑93 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  1:2263577/93‑1 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     >  1:7481347/1‑93 (MQ=255)

      gtaNCTGAAAAACTTCGCAGGCGCNNCCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTGTCGTGACAGCGGGAAGCCGGNTTa                                                                                   <  2:7547633/91‑1 (MQ=255)

      ggATCTGAAAAGCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  1:6572965/93‑1 (MQ=255)

      ggATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  1:3639649/93‑1 (MQ=255)

      ggATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  1:7176076/93‑1 (MQ=255)

      ggATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   >  1:844519/1‑93 (MQ=255)

       gATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTAt                                                                                  <  1:5915862/93‑1 (MQ=255)

         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                >  1:7830429/1‑93 (MQ=255)

           tGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCg                                                                              <  1:4364892/93‑1 (MQ=255)

            gAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGt                                                                             <  2:221745/93‑1 (MQ=255)

               aaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTgg                                                                          >  2:1660307/1‑93 (MQ=255)

                 aCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCg                                                                        >  2:7287230/1‑93 (MQ=255)

                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGGGACAGCGGGAAGCCGGTTTATGCCGTTCGCga                                                                       >  1:1254239/1‑93 (MQ=255)

                    tCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCgaga                                                                     >  2:3653274/1‑93 (MQ=255)

                     cGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGaa                                                                    <  2:823349/93‑1 (MQ=255)

                     cGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGaa                                                                    <  2:5404245/93‑1 (MQ=255)

                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAACAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   >  2:7226868/1‑93 (MQ=255)

                        aGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTa                                                                 <  2:4448027/93‑1 (MQ=255)

                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                >  2:939273/1‑93 (MQ=255)

                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                >  1:5730843/1‑93 (MQ=255)

                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                <  1:3164033/93‑1 (MQ=255)

                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  2:4721023/93‑1 (MQ=255)

                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  2:3518350/93‑1 (MQ=255)

                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  1:2379030/93‑1 (MQ=255)

                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGCACCCGGTTTATGCCGTTGGCGAGAACTACAGcc                                                            >  1:3130445/1‑93 (MQ=255)

                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGcc                                                            >  1:2063556/1‑93 (MQ=255)

                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTATGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGcc                                                            >  2:7798687/1‑93 (MQ=255)

                              cgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCa                                                           <  2:1011749/93‑1 (MQ=255)

                               gaccgaccGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAg                                                          >  1:873248/1‑93 (MQ=255)

                               gACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAg                                                          >  1:1108982/1‑93 (MQ=255)

                               gACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAg                                                          >  1:2444874/1‑93 (MQ=255)

                               gACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAg                                                          <  1:6725716/93‑1 (MQ=255)

                                aCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                         >  1:6756817/1‑93 (MQ=255)

                                  cAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgggg                                                       >  1:1980952/1‑93 (MQ=255)

                                   aaCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                      <  1:1437451/93‑1 (MQ=255)

                                    aCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCa                                                     >  1:3234813/1‑93 (MQ=255)

                                     ccGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCaa                                                    >  1:3776899/1‑93 (MQ=255)

                                      cGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAat                                                   <  2:4088693/93‑1 (MQ=255)

                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  >  1:7474143/1‑93 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:5205293/1‑93 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:341677/1‑93 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:6180688/1‑93 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:2778324/1‑93 (MQ=255)

                                         cTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatt                                                <  1:4282167/93‑1 (MQ=255)

                                          tATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatta                                               >  1:7679013/1‑93 (MQ=255)

                                          tATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatta                                               <  2:2719005/93‑1 (MQ=255)

                                          tATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatta                                               <  1:6118802/93‑1 (MQ=255)

                                           aTGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtattat                                              <  2:1364945/93‑1 (MQ=255)

                                           aTGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtattat                                              <  2:1364944/93‑1 (MQ=255)

                                            tGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtc                                             <  2:3252005/93‑1 (MQ=255)

                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            <  2:6833583/93‑1 (MQ=255)

                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            <  1:5845685/93‑1 (MQ=255)

                                              cAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtctc                                           <  1:7510588/93‑1 (MQ=255)

                                              cAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtctc                                           <  1:4000243/93‑1 (MQ=255)

                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         >  1:2525062/1‑93 (MQ=255)

                                                 tACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGcc                                        <  1:3860747/93‑1 (MQ=255)

                                                  gcATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCa                                       <  2:1920283/92‑1 (MQ=255)

                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAggt                                    >  2:919730/1‑91 (MQ=255)

                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGGTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    >  1:5353818/1‑93 (MQ=255)

                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  >  2:1274932/1‑93 (MQ=255)

                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  <  2:7497932/93‑1 (MQ=255)

                                                          aaaGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGcc                               >  2:2906886/1‑93 (MQ=255)

                                                          aaaGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGcc                               >  2:6255073/1‑93 (MQ=255)

                                                              ctctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATa                           <  2:5080790/93‑1 (MQ=255)

                                                               tctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGATAACTAGAGCCAGGGGCAGTGTTATCTCGCCGGTTTCGCGGATaa                          >  2:2437705/1‑93 (MQ=255)

                                                               tctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaa                          >  2:3300490/1‑93 (MQ=255)

                                                               tctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaa                          >  1:7704663/1‑93 (MQ=255)

                                                                ctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaaa                         <  2:7714227/93‑1 (MQ=255)

                                                                 tGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaaaa                        <  2:5731381/93‑1 (MQ=255)

                                                                 tGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaaaa                        >  2:4868564/1‑93 (MQ=255)

                                                                  gAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAt                       <  2:3360886/93‑1 (MQ=255)

                                                                   aaaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAtt                      <  1:2121371/93‑1 (MQ=255)

                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     <  1:5129243/93‑1 (MQ=255)

                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     >  1:432323/1‑93 (MQ=255)

                                                                      gagTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTgg                   <  2:4650479/93‑1 (MQ=255)

                                                                      gagTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTgg                   >  1:4459226/1‑93 (MQ=255)

                                                                      gagTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTgg                   <  2:7432263/93‑1 (MQ=255)

                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 <  1:3894052/93‑1 (MQ=255)

                                                                          ttCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGt               >  1:5832084/1‑93 (MQ=255)

                                                                           tCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtc              <  1:290675/93‑1 (MQ=255)

                                                                           tCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtc              <  2:6508529/93‑1 (MQ=255)

                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            >  2:3161638/1‑93 (MQ=255)

                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            >  2:2931624/1‑93 (MQ=255)

                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            <  1:1646186/93‑1 (MQ=255)

                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           <  2:3873978/93‑1 (MQ=255)

                                                                               gACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCg          <  2:5963451/93‑1 (MQ=255)

                                                                                  aGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCGATAAAATTTGGCTGTCTCCGCaa       >  2:4299182/1‑93 (MQ=255)

                                                                                  aGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCaa       <  1:6462089/93‑1 (MQ=255)

                                                                                   gCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAg      <  1:3023584/93‑1 (MQ=255)

                                                                                   gCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAg      >  1:776889/1‑93 (MQ=255)

                                                                                    cGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg     >  2:2118613/1‑93 (MQ=255)

                                                                                     gggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGc    >  1:475532/1‑93 (MQ=255)

                                                                                     gggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGc    <  1:6347954/93‑1 (MQ=255)

                                                                                      ggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCg   >  1:4186133/1‑93 (MQ=255)

                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  <  2:6071576/93‑1 (MQ=255)

                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  >  2:3547357/1‑93 (MQ=255)

TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387

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