BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,050	G→A	Q446* (CAA→TAA)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,050	0	G	A	100.0%	601.7 / NA	171	Q446* (CAA→TAA)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (94/77); total (94/77)

Only 100 of 172 total aligned reads displayed.

TTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGC  >  NC_000913/4582960‑4583139

tttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACg                                                                                         <  1:2187432/93‑1 (MQ=255)

  tCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                                       <  1:1112558/93‑1 (MQ=255)

  tCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                                       <  1:806838/93‑1 (MQ=255)

  tCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                                       <  1:445320/93‑1 (MQ=255)

    aaCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTc                                                                                     <  2:3171603/93‑1 (MQ=255)

    aaCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTc                                                                                     <  2:4513265/93‑1 (MQ=255)

      ccGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                                   >  1:33212/1‑93 (MQ=255)

       cgcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                                  <  1:4389442/93‑1 (MQ=255)

        gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTg                                                                                 >  2:2442380/1‑93 (MQ=255)

         cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                <  1:3436218/93‑1 (MQ=255)

         cgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                <  1:2499610/93‑1 (MQ=255)

          gcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                               >  2:1615317/1‑93 (MQ=255)

             gACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAg                                                                            >  2:4151495/1‑93 (MQ=255)

             gACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAg                                                                            >  1:434994/1‑93 (MQ=255)

              aCGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGc                                                                           >  1:1293637/1‑93 (MQ=255)

               cGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCa                                                                          >  2:384952/1‑93 (MQ=255)

                 gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAtt                                                                        >  1:2410847/1‑93 (MQ=255)

                   aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTc                                                                      >  2:3649350/1‑93 (MQ=255)

                   aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTc                                                                      >  2:3649351/1‑93 (MQ=255)

                    tCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                                     >  2:1072043/1‑93 (MQ=255)

                    tCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                                     <  1:2969689/93‑1 (MQ=255)

                     cTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                                    <  2:2578677/93‑1 (MQ=255)

                     cTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCca                                                                    <  2:4467645/93‑1 (MQ=255)

                      ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcag                                                                   >  2:1118911/1‑93 (MQ=255)

                       tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagc                                                                  <  2:3800844/93‑1 (MQ=255)

                        ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagca                                                                 <  2:4128713/93‑1 (MQ=255)

                         cTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                                >  1:4100631/1‑93 (MQ=255)

                         cTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                                <  1:384086/93‑1 (MQ=255)

                         cTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                                <  1:556744/93‑1 (MQ=255)

                          tGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                               <  2:722742/93‑1 (MQ=255)

                          tGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                               >  1:1210119/1‑93 (MQ=255)

                           gATAATAATGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                              <  1:4177370/93‑1 (MQ=255)

                           gATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                              >  2:3062174/1‑93 (MQ=255)

                           gATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                              >  2:805492/1‑93 (MQ=255)

                            ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                                             <  2:1970142/93‑1 (MQ=255)

                            ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                                             <  2:2345576/93‑1 (MQ=255)

                             taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctct                                                            >  2:4713578/1‑93 (MQ=255)

                             taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctct                                                            <  2:1497525/93‑1 (MQ=255)

                             taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctct                                                            <  1:478335/93‑1 (MQ=255)

                             taataaCGCAGGCCCAGCTCGCTCATGCCAGGGGTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCCGctct                                                            >  2:403151/1‑93 (MQ=255)

                                taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTcc                                                         >  1:1386272/1‑93 (MQ=255)

                                   cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCggg                                                      <  2:1986425/93‑1 (MQ=255)

                                    gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                                     >  1:4148465/1‑93 (MQ=255)

                                    gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                                     >  2:2294025/1‑93 (MQ=255)

                                        gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTg                                                 <  1:4027788/93‑1 (MQ=255)

                                            aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCAc                                             >  2:279048/1‑93 (MQ=255)

                                            aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCAc                                             >  2:128072/1‑93 (MQ=255)

                                            aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCAc                                             >  1:221708/1‑93 (MQ=255)

                                            aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCAc                                             <  1:1592706/93‑1 (MQ=255)

                                             gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACt                                            >  1:3219258/1‑93 (MQ=255)

                                               tcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCg                                          >  2:1402762/1‑93 (MQ=255)

                                                cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACAGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTccg                                         >  1:2552406/1‑91 (MQ=255)

                                                  ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCgggg                                       >  1:300018/1‑92 (MQ=255)

                                                  ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggga                                       >  2:4107993/1‑92 (MQ=255)

                                                  ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGt                                       <  2:1556123/93‑1 (MQ=255)

                                                   tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                                      >  2:126296/1‑93 (MQ=255)

                                                    cATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                                     >  1:587000/1‑93 (MQ=255)

                                                    cATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                                     <  2:1972345/93‑1 (MQ=255)

                                                     aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                                    >  1:1574373/1‑93 (MQ=255)

                                                      tGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAgg                                   <  1:242849/93‑1 (MQ=255)

                                                       gCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGGAAGGc                                  >  1:1066806/1‑93 (MQ=255)

                                                        ccAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCt                                 >  1:2025973/1‑93 (MQ=255)

                                                         cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                <  2:2861775/93‑1 (MQ=255)

                                                         cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                <  2:2209466/93‑1 (MQ=255)

                                                          aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                               >  2:1762077/1‑93 (MQ=255)

                                                          aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                               >  2:559634/1‑93 (MQ=255)

                                                           gggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                              <  2:3891359/93‑1 (MQ=255)

                                                           gggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc                              >  1:2762917/1‑93 (MQ=255)

                                                              ttATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc                           >  1:1702991/1‑93 (MQ=255)

                                                               tATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACa                          <  1:4600209/93‑1 (MQ=255)

                                                                aTCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACat                         >  1:1418851/1‑93 (MQ=255)

                                                                 tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACata                        <  1:4173729/93‑1 (MQ=255)

                                                                 tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACata                        >  2:3042568/1‑93 (MQ=255)

                                                                   ggCAAACCACTGTTTTTTTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt                      <  1:4022134/93‑1 (MQ=255)

                                                                   ggCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt                      >  1:4440003/1‑93 (MQ=255)

                                                                   cgCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAtt                      <  1:1348358/92‑1 (MQ=255)

                                                                    gCAAACCACTGCTTCTGTTTTTACGGTTCGCTGCCGAGCCTTTCCAGCAGCTCTTCCGGCGGGTGCCACCCGGGTGAGGCTTTCGACATATTg                     >  1:45901/1‑93 (MQ=255)

                                                                     cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc                    >  2:2494603/1‑93 (MQ=255)

                                                                     cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc                    <  1:1749222/93‑1 (MQ=255)

                                                                      aacccACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCg                   <  2:4210513/90‑1 (MQ=255)

                                                                       aaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCgg                  >  2:1787510/1‑93 (MQ=255)

                                                                       aaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCgg                  >  1:2101015/1‑93 (MQ=255)

                                                                         ccACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGt                <  2:448958/93‑1 (MQ=255)

                                                                          cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTa               <  2:734098/93‑1 (MQ=255)

                                                                          cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTa               <  1:3858071/93‑1 (MQ=255)

                                                                           aCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTAt              >  2:3525376/1‑93 (MQ=255)

                                                                           aCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTAt              >  1:3963072/1‑93 (MQ=255)

                                                                             tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATca            >  1:3570712/1‑93 (MQ=255)

                                                                             tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATca            <  1:4267654/93‑1 (MQ=255)

                                                                             tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATca            <  2:227786/93‑1 (MQ=255)

                                                                              gATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcac           <  1:603416/93‑1 (MQ=255)

                                                                               aTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacg          <  1:3052793/93‑1 (MQ=255)

                                                                                ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgc         >  1:376566/1‑93 (MQ=255)

                                                                                ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgc         >  1:3044544/1‑93 (MQ=255)

                                                                                ttCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgc         >  2:2820791/1‑93 (MQ=255)

                                                                                 tCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgca        >  1:1729365/1‑93 (MQ=255)

                                                                                   tGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgcacg      >  1:2280719/1‑93 (MQ=255)

                                                                                   tGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgcacg      <  2:2718860/93‑1 (MQ=255)

                                                                                       tttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCAGGTcgc  >  1:2747123/1‑93 (MQ=255)

                                                                                       tttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcgc  >  2:2121620/1‑93 (MQ=255)

TTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGC  >  NC_000913/4582960‑4583139

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