BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,600,078	C→T	W135* (TGG→TAG)	yjjP ←	putative structural protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,600,078	0	C	T	100.0%	1076.6 / NA	274	W135* (TGG→TAG)	yjjP	putative structural protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (157/117); total (157/117)

Only 100 of 276 total aligned reads displayed.

GCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATG  >  NC_000913/4599988‑4600167

 cAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATc                                                                                        <  2:1013624/93‑1 (MQ=255)

  aaaGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCt                                                                                       >  2:3139538/1‑93 (MQ=255)

  aaaGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCt                                                                                       >  2:6632172/1‑93 (MQ=255)

  aaaGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCt                                                                                       >  1:2243514/1‑93 (MQ=255)

   aagaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCtt                                                                                      >  2:807969/1‑93 (MQ=255)

     gaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTgg                                                                                    >  2:6265098/1‑93 (MQ=255)

       agGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGt                                                                                  <  1:4373590/93‑1 (MQ=255)

       agGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGt                                                                                  >  1:46297/1‑93 (MQ=255)

        ggTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTa                                                                                 <  1:391549/93‑1 (MQ=255)

         gTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTaa                                                                                <  1:3762909/93‑1 (MQ=255)

               aCGGCACCATCCCAGCCACCGTTATTGAGTTTAGAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                          >  2:4905453/1‑93 (MQ=255)

               aCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                          <  1:4932454/93‑1 (MQ=255)

                cGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATg                                                                         <  2:2776605/93‑1 (MQ=255)

                 ggCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATgg                                                                        >  1:5325493/1‑93 (MQ=255)

                 ggCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATgg                                                                        <  1:1911307/93‑1 (MQ=255)

                  gCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGc                                                                       >  1:1576659/1‑93 (MQ=255)

                  gCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGc                                                                       >  1:3385295/1‑93 (MQ=255)

                  gCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGc                                                                       <  2:291825/93‑1 (MQ=255)

                  gCACCATCCCAGACACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGc                                                                       >  1:1606214/1‑93 (MQ=255)

                   cACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCt                                                                      <  1:6353404/93‑1 (MQ=255)

                   cACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCt                                                                      >  2:2549081/1‑93 (MQ=255)

                    aCCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGTGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTg                                                                     >  2:4155001/1‑93 (MQ=255)

                     ccATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGa                                                                    >  2:5443212/1‑93 (MQ=255)

                      cATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAGCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGaa                                                                   <  2:3529396/93‑1 (MQ=255)

                      cATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGaa                                                                   <  2:6375798/93‑1 (MQ=255)

                      cATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGaa                                                                   >  2:5751517/1‑93 (MQ=255)

                      cATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGaa                                                                   >  2:3287840/1‑93 (MQ=255)

                       aTCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAt                                                                  >  2:4442008/1‑93 (MQ=255)

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 <  1:4224394/93‑1 (MQ=255)

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 <  2:5917574/93‑1 (MQ=255)

                         cccAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAttt                                                                <  2:6043785/93‑1 (MQ=255)

                         cccAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAttt                                                                >  1:1962289/1‑93 (MQ=255)

                          ccAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTg                                                               >  2:227648/1‑93 (MQ=255)

                            aGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGc                                                             >  1:6448367/1‑93 (MQ=255)

                            aGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGc                                                             >  2:3161758/1‑93 (MQ=255)

                             gCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCt                                                            >  2:3352737/1‑93 (MQ=255)

                             gCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCt                                                            >  2:1070748/1‑93 (MQ=255)

                               cACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTCGGCTaa                                                          >  1:283785/1‑93 (MQ=255)

                                aCCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTaaa                                                         >  2:5512084/1‑93 (MQ=255)

                                aCCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTaaa                                                         >  1:6010982/1‑93 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        <  1:5287859/93‑1 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        <  1:6152776/93‑1 (MQ=255)

                                 ccGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        >  2:5016171/1‑93 (MQ=255)

                                  cGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATc                                                       <  1:3277671/93‑1 (MQ=255)

                                  cGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATc                                                       <  2:6453415/93‑1 (MQ=255)

                                   gTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCg                                                      >  1:5881199/1‑93 (MQ=255)

                                   gTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCg                                                      >  2:2040673/1‑93 (MQ=255)

                                   gTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCg                                                      >  1:120015/1‑93 (MQ=255)

                                    ttATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGt                                                     <  2:1768851/93‑1 (MQ=255)

                                     tATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGtt                                                    <  1:2561188/93‑1 (MQ=255)

                                     tATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGtt                                                    >  1:4713553/1‑93 (MQ=255)

                                     tATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGtt                                                    <  2:4645903/93‑1 (MQ=255)

                                       ttGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtc                                                  <  1:816307/93‑1 (MQ=255)

                                       ttGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtc                                                  <  2:4586178/93‑1 (MQ=255)

                                       ttGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtc                                                  <  2:4770310/93‑1 (MQ=255)

                                         gAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtctc                                                >  2:5266453/1‑93 (MQ=255)

                                          aGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTctct                                               <  2:1316094/93‑1 (MQ=255)

                                          aGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTctct                                               >  1:4783502/1‑93 (MQ=255)

                                            tttACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTAc                                             <  1:2179370/93‑1 (MQ=255)

                                             ttACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACg                                            <  2:2689538/93‑1 (MQ=255)

                                              tACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGc                                           >  1:4386201/1‑93 (MQ=255)

                                               aCAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGcc                                          <  2:3372696/93‑1 (MQ=255)

                                               aCAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGcc                                          >  2:5139981/1‑93 (MQ=255)

                                                cAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCt                                         <  2:2650107/93‑1 (MQ=255)

                                                  gAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCttt                                       >  1:3067817/1‑93 (MQ=255)

                                                  gAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCttt                                       <  1:3594266/93‑1 (MQ=255)

                                                    aaCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGt                                     >  2:301822/1‑93 (MQ=255)

                                                       aGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAt                                  <  1:2939792/93‑1 (MQ=255)

                                                       aGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAt                                  >  2:4961229/1‑93 (MQ=255)

                                                         gcgcAAGAAAGGCCAACCATTAAGGCTACCAGCTTTCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATcc                                >  2:4714772/1‑93 (MQ=255)

                                                           gcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcag                              >  1:2509220/1‑93 (MQ=255)

                                                             aaGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagca                            >  2:2174689/1‑93 (MQ=255)

                                                               gAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAga                          >  1:5823860/1‑93 (MQ=255)

                                                               gAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAga                          >  2:2937489/1‑93 (MQ=255)

                                                               gAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAga                          <  1:4498759/93‑1 (MQ=255)

                                                                aaaGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgat                         >  2:2083703/1‑93 (MQ=255)

                                                                aaaGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgat                         <  2:3904215/93‑1 (MQ=255)

                                                                 aaGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatg                        >  1:5694923/1‑93 (MQ=255)

                                                                  aGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatga                       <  2:3026089/93‑1 (MQ=255)

                                                                  aGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatga                       >  2:223693/1‑93 (MQ=255)

                                                                    gCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGatgatg                     <  1:1050705/93‑1 (MQ=255)

                                                                    gCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGatgatg                     >  1:3986769/1‑93 (MQ=255)

                                                                    gCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGatgatg                     >  2:4335592/1‑93 (MQ=255)

                                                                    gCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGatgatg                     >  2:3765775/1‑93 (MQ=255)

                                                                       aaCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTc                  >  2:2448463/1‑93 (MQ=255)

                                                                        aCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTcc                 >  2:1481257/1‑93 (MQ=255)

                                                                        aCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTcc                 <  2:3046709/93‑1 (MQ=255)

                                                                        aCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTcc                 >  1:1384682/1‑93 (MQ=255)

                                                                          cATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGc               >  1:5063095/1‑93 (MQ=255)

                                                                           aTTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa              >  2:3338344/1‑93 (MQ=255)

                                                                             tAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAg            >  2:672742/1‑93 (MQ=255)

                                                                               aGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGaa          <  1:3857226/93‑1 (MQ=255)

                                                                               aGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGaa          >  1:2115828/1‑93 (MQ=255)

                                                                                ggCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAAt         <  1:1762160/93‑1 (MQ=255)

                                                                                 gCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATc        >  1:1154833/1‑93 (MQ=255)

                                                                                 gCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATc        >  1:5481235/1‑93 (MQ=255)

                                                                                  cTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATca       >  2:198255/1‑93 (MQ=255)

                                                                                     ccAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACaa    <  2:364413/93‑1 (MQ=255)

                                                                                       aGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg  >  2:6509374/1‑93 (MQ=255)

                                                                                       aGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg  <  1:693239/93‑1 (MQ=255)

GCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATG  >  NC_000913/4599988‑4600167

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