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%	1314.4 / NA	367	W135* (TGG→TAG)	yjjP	putative structural protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (186/181); total (186/181)

Only 100 of 372 total aligned reads displayed.

GCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATG  >  NC_000913/4599988‑4600167

gCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTAt                                                                                         >  1:407498/1‑93 (MQ=255)

  aaaGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCt                                                                                       <  1:6090811/93‑1 (MQ=255)

  aaaGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCt                                                                                       <  2:1605701/93‑1 (MQ=255)

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

  aaaGAAGGGGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAGCCATTAAGGCTACCAGCTATCt                                                                                       <  1:8094539/93‑1 (MQ=255)

    agaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTg                                                                                     >  1:5283429/1‑93 (MQ=255)

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

          tgatgaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAAc                                                                               >  2:3562746/1‑93 (MQ=255)

              gaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTaa                                                                           <  1:13943/93‑1 (MQ=255)

              gaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTaa                                                                           >  2:754096/1‑93 (MQ=255)

               aCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                          >  1:4909775/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                    aCCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTg                                                                     >  2:4099386/1‑93 (MQ=255)

                     ccATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTCACGTAATGGCTGa                                                                    >  1:4278290/1‑93 (MQ=255)

                     ccATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGa                                                                    <  2:7184240/93‑1 (MQ=255)

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

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

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 >  1:3925330/1‑93 (MQ=255)

                        tCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                 >  1:1037518/1‑93 (MQ=255)

                         cccAGCCACGGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAttt                                                                <  2:8017383/93‑1 (MQ=255)

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

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

                          ccAGCCTCCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAATGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTg                                                               <  2:1311935/93‑1 (MQ=255)

                           cAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTgg                                                              >  2:8710870/1‑93 (MQ=255)

                            aGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGc                                                             <  2:1086776/93‑1 (MQ=255)

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

                             gCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCt                                                            <  2:1361084/93‑1 (MQ=255)

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

                              ccACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTa                                                           <  2:349365/93‑1 (MQ=255)

                               cACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTaa                                                          >  2:1791016/1‑93 (MQ=255)

                               cACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTaa                                                          >  1:7255839/1‑93 (MQ=255)

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

                                 ccGTTATTGAGTTTACAGAAACAGGCGCTAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAAt                                                        >  2:6430843/1‑93 (MQ=255)

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

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

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

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

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

                                    ttATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGt                                                     >  1:641289/1‑93 (MQ=255)

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

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

                                      aTTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGttt                                                   <  1:985655/93‑1 (MQ=255)

                                        tGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtct                                                 <  1:3729758/93‑1 (MQ=255)

                                         gAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtctc                                                <  2:2076550/93‑1 (MQ=255)

                                           gTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTa                                              >  1:7687544/1‑93 (MQ=255)

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

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

                                              tACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGc                                           <  1:7776478/93‑1 (MQ=255)

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

                                                   aaaCAGGCGCAAGAAAGGCCAACCATTAAGGCGACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTg                                      <  2:3115335/93‑1 (MQ=255)

                                                    aaCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTga                                     >  1:2730426/1‑92 (MQ=255)

                                                     aCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTa                                    <  2:7667084/93‑1 (MQ=255)

                                                     aCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTa                                    >  2:8419383/1‑93 (MQ=255)

                                                      cAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTaa                                   <  1:2128234/93‑1 (MQ=255)

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

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

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                 <  2:2162991/93‑1 (MQ=255)

                                                        ggCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                 <  1:339211/93‑1 (MQ=255)

                                                         gcgcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATcc                                <  2:5043635/93‑1 (MQ=255)

                                                           gcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcag                              <  1:8052204/93‑1 (MQ=255)

                                                           gcAAGAAAGGCCAACCATTAAGGCCACCAGCTATCTTGGGTAACGTAATGGCTGAATTTAGCTAAATCGTTTCTCTACGCCTTTGTAATCcag                              >  1:5890511/1‑93 (MQ=255)

                                                              aGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagcag                           >  2:2289865/1‑93 (MQ=255)

                                                              aGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagcag                           >  2:5419765/1‑93 (MQ=255)

                                                              aGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagcag                           >  1:657852/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                                       aaCCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTc                  <  2:5895060/93‑1 (MQ=255)

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

                                                                        aCCATTAATGCTATCAGCTTTTTTGGGTAACGTTATGGCTGAATTTGGCTAAATTGTTTCTCTATGTCTTTGTTATCCAGCAGATGATGCTcc                 >  2:8077577/1‑93 (MQ=255)

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

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

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

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

                                                                          cATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGc               <  2:1113068/93‑1 (MQ=255)

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

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

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

                                                                              aaGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGa           >  1:7979965/1‑93 (MQ=255)

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

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

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

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

                                                                                    aCCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATcaca     >  2:6092702/1‑93 (MQ=255)

                                                                                    aCCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATcaca     <  2:7778487/93‑1 (MQ=255)

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

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

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

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

GCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATG  >  NC_000913/4599988‑4600167

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